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Preface 



One of the most important things in one person's life is to have good friends, 
whom one can exchange thoughts and ideas with. The present volume is a col- 
lection of contributions by friends of Holger Bech Nielsen for his 60 th birthday. 
The collection of papers covers a broad range of physics, which shows that Hol- 
ger Bech Nielsen has been working and is working on many topics — also be- 
cause different fields of physics are much more connected as one would believe, 
at least when fundamental questions of physics are concerned. The most fasci- 
nating works in physics bring new understanding of known phenomena, since 
they lead to theories beyond the known ones. Holger Bech Nielsen is one of (very 
rare and because of that additionally precious ) scientists who is always willing 
to play with new thoughts and progressive ideas, trying to formulate them into 
formal proofs built on well defined assumptions and into formulas, pointing out 
in such a way that mathematics is a part of Nature. 

It is my pleasure and my privilege to have Holger as a real friend in the last 
five years, since we have started — together with Colin Froggatt — the annual 
workshop at Bled, Slovenia, which is entitled "What comes Beyond the Standard 
Model". During the workshop the ideas connected with the open problems of (el- 
ementary particle) physics and cosmology are very vividly and openly discussed. 
The first volume of this Proceedings is a Festschrift dedicated to Holger. I wish to 
thank all his friends who are contributing to this volume and also to those friends, 
who have sent their contributions too late to be included. 

The second volume is collecting contributions and discussions from the last two 
workshops and will appear with a little delay. 

The editors would like to thank Yasutaka Takanishi for a lot of work, which con- 
tributed to make the Volume 1 of the Proceedings see the light of day. 

Dear Holger: I am honoured to congratulate you for your 60th jubilee in the name 
of all your friends contributing to this volume and also in the name of the others, 
wishing you many additional fruitful years, health, wealth and success. 
In the name of your friends: 

Norma Susana Mankoc Borstnik Ljubljana, December 2001 
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Unified Internal Space of Spins and Charges 

N. Mankoc Borstnik* 

Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, Ljubljana and 
Primorska Institute for Natural Sciences and Technology, C. Marezganskega upora 2, 
Koper 6000, Slovenia 

Abstract. Can the assumption that the spin and all the charges of either fermions or 
bosons unify, help to find the answers to the open questions of the electroweak Standard 
model? The approach is presented in which polynomials in Grassmann space are used to 
describe all the internal degrees of freedom of spinors, scalars and vectors, that is their 
spins and charges j6po|]. The same can be achieved also by polynomials of differen- 
tial forms. If the space - ordinary and anti-commutative - has 14 dimensions or more, the 
appropriate spontaneous breaking of symmetry leads gravity in d dimensions to manifest 
in four-dimensional subspace as ordinary gravity and all the gauge fields as well as the 
Yukawa couplings. The approach manifests four generations of massless fermions, which 
are left handed SU(2) doublets and right handed SU(2) singlets. 

1 Introduction 

It can not be any doubt that the internal space of spins and charges plays an 
important role in our world, so important as the ordinary space of coordinates 
and momenta does. Without the internal space of spins and charges no spinors 
(fermions), no vectors (gauge fields), no tensors (gravity) would exist, scalars (if 
there are any elementary scalars at low energies) would not interact and accord- 
ingly no matter could exist. We have shown [|6|^U^,Pl|llp|] how a space of anti- 
commuting coordinates can be used to describe spins and charges of not only 
fermions but also of bosons, unifying spins and charges for either fermions or 
bosons and that spins in d-dimensional space manifest (at low energies) as spins 
and charges in four-dimensional space- time, while accordingly gravity in d-dimen- 
sional space - after the appropriate breaking of symmetry, such as 

SO(1, 13) 

i 

SO(1 ,7) ®SO(6) 

/ \ 
SO(1, 7) SU(3)®U(1) 

i i 
SO(1 ,3) ®SU(2) ®U(1) SU(3)®U(1) 

V v ' 

SO(1 , 3) ® SU(2) <g> U(1 ) <8> SU(3) 
* e-mail:norma.s.mankocs@ijs.si 
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- manifests in four-dimensional subspace as the ordinary gravity and all the (known) 
gauge fields. 

The knowledge of the way how our universe (or universes) "has (or have) 
made a choice" of the signature of space-time and of the way of breaking symme- 
tries of (external and internal) space-time from SO(q, d — q) down to S0(1 , 3) x 
SU(3) x SU(2) x U(1) can help to answer the open questions of the Standard 
Model, such as: 

i) Why only left handed spinors carry the weak charge while right handed 
spinors are weak charge-less? 

ii) Why besides the Planck scale (at least) the weak scale occurs? 

iii) Why there are spin (SO(1,3)) and charge (SU(3), SU(2), U(1 )) internal 
degrees of freedom? 

iv) Why there are besides the gravity (gauge) field also the charge gauge 
fields? 

v) Where do the families come from? 

vi) Where the Higgs and the Yukawa couplings come from? 
and many others. 

If the break of let us say SO ( 1 , 1 3) occurred in the way presented in the above 
diagram, then we could say why left handed weak charge doublets and right 
handed weak charge singlets appear in the same multiplet ( in this case not only 
the weak charge scale but also the intermediate scale before the Planck scale) 
exists, and why there are (four rather than three) families of quarks and leptons, 
as well as that there is gravity in d-dimensional space, which manifests at low 
energies as the ordinary gravity and all the (known) gauge fields, the Higgs and 
the Yukawa couplings. 

In this article, we briefly present the approach, which unifies spins and charg- 
es, gravity and gauge fields, leads to multiplets of left handed weak charged 
spinors and of right handed weak chargless spinors, to families of quarks and 
leptons, pointing out to questions like why and how symmetries break, or why 
and how signatures are chosen. 

We present a possible Lagrange function for a free particle and the quan- 
tization of anti-commuting coordinates [B^ ^lE5 ]. Introducing vielbeins and 
spin connections, we demonstrate on the level of a covariant momentum how 
the spontaneous breaking of symmetry might lead to the symmetries of the Stan- 
dard model. This part of the work has been started together with A. Borstnik 
and is continuing with Holger Bech Nielsen. We show how the symmetry of 
the group SO(l, 13) breaks to SO (1,7) (leading to multiplets with left handed 
SU(2) doublets and right handed SU(2) singlets) and to SO (6), which then leads 
to the SO(1, 3) x SU(2) x U(1) x SU(3) x U(1 ). The two U(1 ) symmetries enable 
besides the hyper-charge, needed in the Standard Model, an additional hyper- 
charge, which is nonzero for a right handed SU(2) and SU(3) singlet, like it is a 
right handed neutrino in the Standard model. In the last four years, the author 
of the paper organized together with Holger Bech Nielsen and Colin Froggatt 
annual workshops at Bled, Slovenia, entitled "What comes beyond the Standard 
Model?", the real workshops in which we discuss all the open questions of the 
Standard Model, as well as the proposed approaches which might lead to physics 
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beyond the Standard Model. Very open discussions leaded to many new ideas 
and suggestions, so that we all profit from these workshops a lot. Not only one 
sees connections and correspondences between different approaches, but we also 
open new questions, new problems, trying to find new solutions to the problems. 

Thoughts and ideas, presented in this paper, although started and developed 
by the approach of unification of spins and charges [pi, have been enriched and 
developed under the influence of these discussions. Some of the discussions led 
to common papers with the main opponent Holger Bech Nielsen, some papers 
are in preparation, with the main opponent and with others. 

2 Dirac equation in ordinary space and in space of 
anti-commutative coordinates 

What we call quantum mechanics in Grassmann space is the model for going be- 
yond the Standard Model with extra dimensions of ordinary and anti-commuting 
coordinates, describing spins and charges of either fermions or bosons in a unique 
way i0|§i|l|l]. 

In a d-dimensional space-time the internal degrees of freedom of either spinors 
or vectors and scalars come from the Grassmann odd variables 

d a , q€{0,1,2,3,5,--- ,d}. 

We write wave functions describing either spinors or vectors in the form 

<e a |o>= Y_ Y. a -i>^ a i e ai e a2 ---e ai , 

i=0,l ,..,3,5,..,d {ai <a 2 <...<a i }G{0,l ,..,3,5,..,d} 

(1) 

where the coefficients a ai Q2) . ai depend on commuting coordinates x a , a 6 
{0,1,2, 3, 5, .., d}. The wave function space spanned over Grassmannian coordi- 
nate space has the dimension 2 d . Completely analogously to usual quantum me- 
chanics we have the operator for the conjugate variable Q a to be 

Pa = -i^a. (2) 

The right arrow tells, that the derivation has to be performed from the left hand 
side. These operators then obey the odd Heisenberg algebra, which written by 
means of the generalized commutators 

{A, B} := AB — (—1 ) nAB BA, (3) 

where 

J +1 , if A and B have Grassmann odd character 
nAB \ 0, otherwise, 

takes the form 

{p 9a ,p eb } = o = {e-,e b }, { P fltt ,8 t } = -tn ab . (4) 
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Here r) ab is the flat metric r\ = diag{1 , — 1 , — 1 , ...}. 
We may define the operators 

a a :=i(p 9Q -ie Q ), a Q :=-(p ea + i6 Q ), (5) 

for which we can show that the a a 's among themselves fulfill the Clifford algebra 
as do also the a a 's, while they mutually anti-commute: 

{a a ,a b } = 2ri ab ={a a ,a b }, {a a ,a b } = o. (6) 

We could recognize formally 

either a Q p a |O>=0, or a Q p Q |(D >= (7) 

as the Dirac-like equation, because of the above generalized commutation rela- 
tions. Applying either the operator a a p Q on the left-hand side equation or a a p Q 
on the right-hand side equation we get the Klein-Gordon equation p a p a | c I ) >— 0, 
where we define p a = i-af^"- 

One can check that none of the two equations (0) have solutions which would 
transform as spinors with respect to the generators of the Lorentz transforma- 
tions, when taken in analogy with the generators of the Lorentz transformations 
in ordinary space (L ab = x a p b — x b p a ) 

S ab : = 9 a p eb -e b p 9a . (8) 

But we can write these generators as the sum 

S ab = sab + §ab ) gab . = J [& a > & b] > §ab . = }_ [fla -b] ^ (9) 

with [A, B] := AB — BA and recognize that the solutions of the two equations (^) 
now transform as spinors with respect to either S ab or S ab . 

One also can easily see that the untilded, the single tilded and the dou- 
ble tilded S ab obey the d-dimensional Lorentz generator algebra {M ab , M cd } = 
-i(M ad ri bc + M bc ri ad - M QC u bd - M bd u QC ), when inserted for M ab . 

Kahler formulated spinors [||] in terms of wave functions which are super- 
positions of the p-forms in the (d = 4)- dimensional space. A general linear com- 
bination of p-forms follows from Eq.(Q) if replacing a by dx a A. 

We presented in ref . [ p^ ] the parallelism between our approach and the Kahler 
approach. We also presented in the same reference the generalization of the Kahler 
approach, suggested by our approach. In both approaches two types of operators 
fulfilling the Clifford algebra (the ones of our approach are presented in Eqs.(^^)) 
as well as the two Dirac-like equations (Eq.(^) represent our Dirac-like equation) 
can be obtained. Both approaches offer the generators of the Lorentz transforma- 
tions, describing not only spinors but also vectors (Eq.(||)). In both approaches the 
y a matrices, fulfilling the Clifford algebra and having an Grassmann even char- 
acters (which assures that y a 's transform Grassmann odd object into Grassmann 
odd objects and accordingly do not change the Grassmann character of spinors) 
can be defined 



f a =ia°a a ) with {Y Q ,y b } = 2ri 



(10) 
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The "naive" definition of gamma-matrices (jnaive := fl Q ), which changes the 
Grassmann character of spinors, differs from the Grassmann even definition of 
gamma-matrices, presented in Eq.(|lO|), which keeps the Grassmann character of 
spinors (both fulfilling the Clifford algebra), only when y°-matrix has to sim- 
ulate the parity reflection which is — > —0. In all physical applications (such 
as construction of currents) the two definitions can not be distinguished among 
themselves, since y a 's always appear in pairs. We can check that the y a (Eq.(|lO|)) 
indeed perform the operation of the parity reflection. 



2.1 Scalar product 

In our approach [H0,|fl||[ll|] the scalar product between the two functions < 9 a | O -\ > 
and < 9 Q |3>2 > is defined 



< <Dil® 2 >= 
and cu is a weight function 



d d (to < 6 a |(Di >) < e Q |o 2 > 



a>= j f (B'+d' 1 ) 

i=0,1 ,..,d 



which operates on only the first function < a |Oi > and 



d0 a = 0, 



= i, d a = e^.j 1 © 



According to the above definition and Eq.(|]) it follows 

0,d «i <0C2<..<0Cd 



(11) 



(12) 



2.2 Four copies of Weyl bi-spinors 

We present in this subsection four copies of two- Weyl spinors. 
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Table I: The polynomials of m , representing the four times two Weyl spinors, 
are presented. For each state the eigenvalues of S 12 , S 03 , T' 4 ' :— iQ ^ 1 d 2 a 3 are 
written. The Roman numerals tell the possible family number. We use the relation 
a Q |0>=6 Q . 

We present here for d = 4 the 2 d vectors, which we arrange into four copies 

of two Weyl spinors, one left ( < H 4 > >= -1, T< 4 > = il^lie abcd 5 ab 5 cd ) 
and one right ( < T' 4 ' >= 1) handed in such a way that they are at the same 
time also the eigenvectors of the operators S 12 and the S 03 and have either an 
odd or an even Grassmann character. We have made a choice of (~) operators, 
putting the operators of the type (" ) equal to zero. We present these vectors as 
polynomials of 9 m 's, m S (0,1,2,3). The corresponding Kahler's p-forms fol- 
low if 9 a 's are replaced by dx a A. The two Weyl vectors of one copy of the Weyl 
bi-spinors are connected by the y m (Eq.(p^)) operators, while the two copies of 
different Grassmann character are connected by a a , respectively. The two copies 
of an even Grassmann character are connected by the ( a kind of a time reversal 
operation) 9° — > —9° (or equivalently dx° — > — dx°), if differential forms are 
concerned. 

We present in Table I four copies of the Weyl two spinors as polynomials of 
9 a . Eigenstates are orthonormalized according to the scalar product of Eq.(|l2|) 

Analyzing the irreducible representations of the group SO(1 , 3) with respect 
to the generator of the Lorentz transformations of the vectorial type [^!|8lPUlll3] 
(Eqs.( ||)) one finds for d = 4 two scalars (a scalar and a pseudo scalar), two three 
vectors (in the complex version of the SU(2) x SU(2) representation of SO(1 ,3) 
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denoted by (1,0) and (0,1) representation, respectively, with < P' 4 ' >= ±1) 
and two four vectors. One can find the polynomial representation for this case 
inref.0. 

2.3 Generalization to extra dimensions 

It has been suggested [||] that the Lorentz transformations in the space of 9 a 's in 
(d — 4)-dimensions manifest themselves as generators for charges observable for 
the four dimensional particles. Since both, the extra dimensional spin degrees of 
freedom and the ordinary spin degrees of freedom, originate from the 9 a 's (or the 
forms) we have a unification of these internal degrees of freedom. 
Let us take as an example the model [|J|2|] which has d = 1 4 and at first - at the 
high energy level - SO (1,13) Lorentz group, but which should be broken ( in two 
steps ) to first SO(1,7) x SO (6) and then to SO(1,3) x SU(3) x SU(2). We shall 
comment on this model in section [|. 

2.4 Appearance of spinors 

By exchanging the Lorentz generators S ab by the S ab say ( or the S ab if we choose 
them instead), of Eq.(^), a spinor field appears out of models with only scalar, 
vector and tensor objects. One of the two kinds of operators fulfilling the Clifford 
algebra and anticommuting with the other kind - it has been made a choice of 
a a in our approach and similarly one also can proceed in the Kahler case - are 
put to zero in the operators of the Lorentz transformations; as well as in all the 
operators representing physical quantities. The use of a in the operator y° (and 
equivalently also in the Dirac case) is the exception, only used to simulate the 
Grassmann even parity operation — > —0 (or for p-forms dx a — > — dx a ). 

3 Lagrange function for a free massless particles in ordinary 
and in Grassmann space and canonical quantization 

We present in this section the Lagrange function for a particle which lives in a 
d-dimensional ordinary space of commuting coordinates and in a d-dimensional 
Grassmann space of anti-commuting coordinates X a = {x a , 8 a } and has its geode- 
sies parameterized by an ordinary Grassmann even parameter (t) and a Grass- 
mann odd parameter(£,). We derive the Hamilton function and the correspond- 
ing Poisson brackets and perform the canonical quantization, which leads to the 
Dirac equation with operators presented in section |2|. 

The coordinates X a = X a (x a , 9 a , t, £,) are called the super-coordinates. We 
define the dynamics of a particle by choosing the action (in complete analogy 
with the usual definition of the scalar product in ordinary space [||,|]] 

-J 



dTdtEEX3iX a 4 9 i x Vbil 



where di := (3 T , 9 fJ.T 1 = (t, £,), while E A determines a metric on a two dimen- 
sional super-space T x , E = det(E A ) . We choose r| AA = 0,r|i2 = 1 = T)2i, while 
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r| a b is the Minkowski metric with the diagonal elements (1,-1,— 1,-1, 1). 
The action is invariant under the Lorentz transformations of super-coordinates: 
X' a — A a bX b . Since a super-matrix E 1 a transforms as a vector in a two-dimen- 
sional super-space t 1 under general coordinate transformations of t 1 , E 1 a ti is 
invariant under such transformations and so is d 2 tE. The action is locally super- 
symmetric. The inverse matrix E A i is defined as follows: EVE B ; = 6 B a- 

Taking into account that either x a or a depend on an ordinary time param- 
eter x and that £, 2 = , the geodesies can be described as a polynomial of £, as 
follows: X a = x a + e£,9 a . We choose e 2 to be equal either to +i or to — i so that 
it defines two possible combinations of super-coordinates. Accordingly we also 
choose the metric E i A : E 1 1 = 1 , E 1 2 = — eM, E 2 ! = £,, E 2 2 = N — e£,M, with 
N and M Grassmann even and odd parameters, respectively. We write A = A, 
for any A. 

If we integrate the above action over the Grassmann odd coordinate d£„ the 
action for a super-particle follows: 

1 If 2 M 

dx (-x a x a + E 2 e Q e Q -^px Q e Q ). (13) 

Defining the two momenta 

p -= w = e20a ' Va := I? = ^ {±a Mp9a) - (14) 

the two Euler-Lagrange equations follow: 

*!_0, *!l = e^. (15) 
dx dT z 

Variation of the action (Eq.(|l3|)) with respect to M and N gives the two constraints 

x i :=p a a e =0) x 2 =p a pa=0> a e :=ip e +£ 2 9q) (16) 

while x 3 a := — Pa + e2 ®a = (Eq.Q)) is the third type of constraints of the 
action(|T3|). For e 2 = — i we find that a e a = d a , which agrees with Eq.(^), while 
X 3 a = &a = 0, which makes a choice between a a and a a . 

We find the generators of the Lorentz transformations for the action(|T3|) to be 



M ab = L ab + S ab ^ j_ab = x a p b_ x b p a ^ gab = e a p eb_ e b p 9a = S^ + S ab , (17) 

which agree with definitions in Eq.(^) and show that parameters of the Lorentz 
transformations are the same in both spaces. 
We define the Hamilton function: 

H := x a p Q + e a p e Q — L = jNp Q p a + ^Mp a (a a + ia a ) (18) 

and the corresponding Poisson brackets 

3A 3B 9A 3B dA d& dA 9F5 
{ ' Jp ~9x^9p7~ 9pT9x^ + d9^9p][ + aplse^' ( ' 
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which fulfill the algebra of the generalized commutators [0] of Eq^. 

If we take into account the constraint x 3 a = da = in the Hamilton func- 
tion (which just means that instead of H the Hamilton function H + Y.i + 
Y- a a 3 a x 3a is taken, with parameters oc l ,i = 1,2 and a 3 a = — ^rP a / a = 0, 1,2,3, 
5, .., d chosen in such a way that the Poisson brackets of the three types of con- 
straints with the new Hamilton function are equal to zero) and in all dynamical 
quantities, we find: 

H = Jn P q Pq + ] -M V a a a , x 1 = V a Va = 0, X 2 = V a &a = 0, (20) 



Pa ={Pa,H} P = 0, d a ={a a ,H}p =lMp a , 

which agrees with the Euler-Lagrange equations (|i~5|). 
We further find 

X i ={H,x i }p=0, i = 1,2, x 3 a={H,x 3 a}p=0, a = 0,1,2,3 ) 5,..,d ) 

which guarantees that the three constraints will not change with the time param- 
eter t and that Ivl ab = 0, with M ab = L ab + S ab , saying that M ab is the constant 
of motion. 

The Dirac brackets, which can be obtained from the Poisson brackets of Eq.(p9[) 
by adding to these brackets on the right hand side a term —{A, d c }p- (— j^r\ ce )- 
{a e , B}p, give for the dynamical quantities, which are observables, the same re- 
sults as the Poisson brackets. This is true also for d a , ( {d a ,d b }D = in ab = 
{a a , d b }p), which is the dynamical quantity but not an observable since its odd 
Grassmann character causes super-symmetric transformations. We also find that 
{a a , d b }D =0 = {d a , d b }p . The Dirac brackets give different results only for the 
quantities 9 Q andp 9a andfor d Q among themselves: {e Q ,p eb } P = r| ab ,{e a ,p eb }D = 
jT| ab , {a a ,a b }p = 2ir| ab ,{d a , d b }D = 0. According to the above properties of 
the Poisson brackets, we suggested [^10| that in the quantization procedure the 
Poisson brackets ( |T^ ) rather than the Dirac brackets are used, so that variables d a , 
which are removed from all dynamical quantities, stay as operators. Then fi a and 
d a are expressible with Q a and p 9a (Eq.(^)) and the algebra of linear operators 
introduced in section ^, can be used. We shall show, that suggested quantization 
procedure leads to the Dirac equation, which is the differential equation in ordi- 
nary and Grassmann space and has all desired properties. 

In the proposed quantization procedure — i{A, B} p goes to either a commuta- 
tor or to an anticommutator, according to the Poisson brackets (|9|). The operators 



9 a ,p 9a (in the coordinate representation they become a — > a , p 9 — > ig§rr) 
fulfill the Grassmann odd Heisenberg algebra, while the operators d a and d a 
fulfill the Clifford algebra (Eq.(§)). 



The constraints (Eqs.(16)) lead to the Weyl-like and the Klein-Gordon equa- 
tions 

p Q a a |0 >=0 , p a p a |0 >=0, withp a a a p b a b =p a p a - (21) 

Trying to solve the eigenvalue problem q q |<P >= 0, a = (0, 1,2,3,5, d), we 
find that no solution of this eigenvalue problem exists, which means that the third 
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constraint fi a = can't be fulfilled in the operator form (although we take it 
into account in the operators for all dynamical variables in order that operator 
equations would agree with classical equations). We can only take it into account 
in the expectation value form 

< 0|a a |0 >= 0. (22) 

Since a a are Grassmann odd operators, they change monomials (Eq.(Q)) of an 
Grassmann odd character into monomials of an Grassmann even character and 
opposite, which is the super-symmetry transformation. It means that Eq.(p2|) is 
fulfilled for monomials of either odd or even Grassmann character and that su- 
perpositions of the Grassmann odd and the Grassmann even monomials are not 
solutions for this system. 
We define the projectors 

P± = ^0±^/(-) tf tt), (P±) 2 =P±, (23) 
where T and T are the two operators defined for any dimension d as follows 

t = i a Ila=0,1, 2,3,5,.. V^> t = i a na=0,1,2,3,5,..,d ^V^, with a 

equal either to d/2 or to (d — 1 )/2 for even and odd dimension d of the space, 
respectively. It can be checked that (T) 2 = 1 = (T) 2 . 

We can use the projector P± of Eq.(|23|) to project out of monomials either the 
Grassmann odd or the Grassmann even part. Since this projector commutes with 
the Hamilton function ({P±,H} = 0), it means that eigenf unctions of H, which 
fulfill the Eq.(^), have either an odd or an even Grassmann character. In order 
that in the second quantization procedure fields |0 > would describe fermions, 
it is meaningful to accept in the fermion case Grassmann odd monomials only. 

4 Particles in gauge fields 

The dynamics of a point particle in gauge fields, the gravitational field in d- 
dimensions, which then, as we shall show, manifests in the subspace d = 4 as 
ordinary gravity and all the Yang-Mills fields, can be obtained by transforming 
vectors from a freely falling to an external coordinate system [|l4j]. To do this, su- 
pervielbeins e a H have to be introduced, which in our case depend on ordinary 
and on Grassmann coordinates, as well as on two types of parameters t 1 = (t, £,). 
The index a refers to a freely falling coordinate system ( a Lorentz index), the 
index \i refers to an external coordinate system ( an Einstein index). 

We write the transformation of vectors as follows diX a = e a lx diX^ L , diX^ — 
f%3iX Q , 9 t = (3 T ,3 £ ). From here it follows that e a ^ b = 5 Q b , f%e% = 5%. 

Again we make a Taylor expansion of vielbeins with respect to e 1 ^ = 

e<v + £ 2 £,e b ev , f% = f% - e 2 ee b f%b. 

Both expansion coefficients again depend on ordinary and on Grassmann 
coordinates. Having an even Grassmann character e a ^ will describe the spin 2 
part of a gravitational field. The coefficients e a ,^b define the spin connections 
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It follows that e VS = 6 a b , f %e% = 5% , e % b ^c = e%f ^ cb . 

We find the metric tensor g^ = e a |x e Q ^,, = f^af 1 " 1 . Rewriting the action 
from section || in terms of an external coordinate system, using the Taylor expan- 
sion of super-coordinates and super-fields e a ,^ and integrating the action over 
the Grassmann odd parameter £,, the action follows 

+ e 2 l(e b e Q -e Q e b )e% b x^}, (24) 

which defines the two momenta of the system p H = —VoiL+jS ab e avL b, V% = 
— i9 a e a ^ ( e 2 = — i ). Here poix are the canonical (covariant) momenta of a particle. 
For p^ = pjjf it follows that p^ is proportional to Q . Then d Q = i(Va — i9a), 
while d Q = 0. We may further write 

POn = P|x — 2^ abea ^ b =Vii— 2^ abuJab H ' ^abix = r ( e an.b — e b|ia ), (25) 

which is the usual expression for the covariant momenta in gauge gravitational 
fields [0. One can find the two constraints 

p£p<V = 0=p o ^f%a Q . (26) 

We shall comment on the breaking of symmetries which leads in (d = 4)- 
dimensional subspace as ordinary gravity and all the gauge fields in section [| 

5 Breaking SO(1,13) through SO(1,7) x SO(6) to 
SO(l, 3) x SU(2) x U(1) x SU(3) 

In this section, we shall first discuss a possible breaking of symmetry, which leads 
from the unified theory of only spins and gravity in d dimensions to spins and 
charges and to the symmetries and assumptions of the Standard Model, on the 



algebraic level < p.l\ ). We shall then comment on the breaking of symmetries on the 



level of canonical momentum for the particle in the presence of the gravitational 



field {§2) 



We shall present as well the possible explanation for that postulate of the 
Standard Model, which requires that only left handed weak charged massless 
doublets and right handed weak charged massless singlets exist, and accordingly 
connect spins and charges of fermions. 

5.1 Algebraic considerations of symmetries 

The algebra of the group SO(l , d — 1 ) or SO(d) contains n subalgebras de- 
fined by operators x Al , A = 1 , n; i = 1 , tia, where Ua is the number of elements 
of each subalgebra, with the properties 

[t a \t b; ] =i5 AB f Aijk T A \ (27) 
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if operators x Ax can be expressed as linear superpositions of operators M ab 

T Al = c Ai ab M Qb , c Ai ah = -c Ai ba , A = 1,n, i=1,n A , a,b = 1,d. (28) 

Here f Al ' k are structure constants of the (A) subgroup with tia operators. 
According to the three kinds of operators S ah , two of spinorial and one of vec- 
torial character, there are three kinds of operators t Ai defining subalgebras of 
spinorial and vectorial character, respectively, those of spinorial types being ex- 
pressed with either S ab or S ab and those of vectorial type being expressed by S ab . 
All three kinds of operators are, according to Eq.(^), defined by the same coef- 
ficients c Al a b and the same structure constants f Al i k . From Eq.(^) the following 
relations among constants c Av a b follow 

-4c Ai Qb c Bib c - 6 AB f Aijk c Ak QC = 0. (29) 

When we look for coefficients c Al a b which express operators t Ai , forming a 
subalgebra SU(n) of an algebra SO(2n) in terms of M ab , the procedure is rather 
simple[||l]]. We find: 

T Am _ _i( & Am). k j M (2j-l)(2k-1) + M (2j)(2k) + iM (2j)(2k-l ) _ iM (2j-1 )(2k) j_ 

2 (30) 
Here (& Am )jic are the traceless matrices which form the algebra of SU(n). One 
can easily prove that operators x Am fulfill the algebra of the group SU(n) for any 
of three choices for operators M Qb : S ab ,S Qb ,S ab . 

While the coefficients are the same for all three kinds of operators, the repre- 
sentations depend on the operators M ab . After solving the eigenvalue problem 
for invariants of subgroups, the representations can be presented as polynomials 
of coordinates a , or dx a A, a = 0, 1 , 2, 3, 5, .., 14. The operators of spinorial char- 
acter define the fundamental representations of the group and the subgroups, 
while the operators of vectorial character define the adjoint representations of 
the groups. We shall from now on, for the sake of simplicity, refer to the polyno- 
mials of Grassmann coordinates only. 

We first analyze the space of 2 d vectors for d = 14 with respect to commut- 
ing operators (Casimirs) of subgroups SO ( 1 , 7) and SO (6), so that polynomials of 
6° , 9 1 , 6 2 , 9 3 , 9 5 , 6 6 , 6 7 and e 8 are used to describe states of the group SO(l,7) and 
then polynomials of© 9 ,© 10 ,© 11 ,© 12 ,© 13 and 14 further to describe states of the 
group SO (6). The group SO(1 , 13) has the rank equal to r = 7, since it has 7 com- 
muting operators (namely for example 6> 01 ,S 23 ,S 56 , ...,<S 13 14 ), while the ranks 
of the subgroups SO ( 1 , 7) and SO (6) are accordingly r = 4 and r = 3, respectively. 
We may further decide to arrange the basic states in the space of polynomials of 
0°, s as eigenstates of 4 Casimirs of the subgroups SO(1 , 3), SU(2), and U(1 ) 
(the first has r = 2, the second and the third have r = 1 ) of the group SO ( 1 , 7), and 

the basic states in the space of polynomials of 9 14 as eigenstates of r = 3 

Casimirs of subgroups SU(3) and U( 1 ) ( with r = 2 and r = 1 , respectively) of the 
group SO(6). 
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We presented in Table I the eight Weyl spinors, two by two - one left ( T' 4 ' = 
— 1) and one right ( T' 4 ' = 1) handed - connected by y m , m = 0, 1 , 2, 3 into Weyl 
bi-spinors. Half of vectors have Grassmann odd (odd products of m ) and half 
Grassmann even character. The two four vectors of the same Grassmann char- 
acter are connected by the discrete time reversal operation 0° — » —8° ( ref.[12|), 
while the two four vectors, which differ in Grassmann character, are connected 
by the operation of a a . 

According to Eqs.(^,^8[^9|), one can express the generators of the subgroups 
SU(2) and U(1) of the group SO(l , 7) in terms of the generators S ab . 

We find (since the indices 0, 1 , 2, 3 are reserved for the subgroup SO(1 , 3)) 



T 31 :=1(S 58 -S 67 ), T 32. = ^ iS 57 +<S 68 )( T 33 . = ^56 _ ^ (31) 

One also finds 

t 41 := l(5 56 +5 78 ). (32) 

The algebra of Eq.(|2^) follows (since the operators t Ai have an even Grass- 
mann character, the generalized commutation relations agree with the usual com- 
mutators, denoted by [ , ]). 

{T 3i ,T 3i } = te ijk T 3 \ {t 4 \t 31 } = 0. (33) 

One notices that t 51 := }(5 58 + S 67 ) and t 52 := \{S 57 - S 68 ) together with 
t 41 form the algebra of the group SU(2) and that the generators of this group 
commute with T 3t . 

We present in Table II the eigenvectors of the operators f 33 and (t 3 ) 2 = 
(t 31 ) 2 + (f 32 ) 2 + (f 33 ) 2 , which are at the same time the eigenvectors of f 41 , for 
spinors. We find, with respect to the group SU(2), two doublets and four singlets 
of an even and another two doublets and four singlets of an odd Grassmann char- 
acter. 
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Table II: The eigenstates of the operators for spinors f 33 ,f 41 are presented. We 
find two doublets and four singlets of an even Grassmann character and two 
doublets and four singlets of an odd Grassmann character. One sees that complex 
conjugation transforms one doublet of either odd or even Grassmann character 
into another of the same Grassmann character changing the sign of the value of 
x 33 , while it transforms one singlet into another singlet of the same Grassmann 
character and of the opposite value of f 41 . One can check that a h , h e (5, 6, 7, 8), 
transforms the doublets of an even Grassmann character into singlets of an odd 
Grassmann character. 

One sees that f 5 \ i = 1,2, transform doublets into singlets (which can easily 
be understood if taking into account that f 5t close together with t 41 the algebra 
of SU(2) and that the two SU(2) groups are isomorphic to the group SO (4)). 

One also sees the following very important property of representations of the 
group SO(1 , 7): If applying the operators S ab , a, b = 0,1,2,3, 5, 6, 7, 8 on the direct 
product of polynomials ofTableland Table II, which forms the representations of the group 
SO(1 , 7), one finds that a multiplet o/SO(1 , 7) exists, which contains left handed SU(2) 
doublets and right handed SU(2) singlets. It exists also another multiplet which con- 
tains left handed SU(2) singlets and right handed SU(2) doublets. It turns out 
that the operators S mh , with m = 0, 1,2,3 and h = 5,6,7,8, although having 
an even Grassmann character, change the Grassmann character of that part of the 
polynomials which belong to Table I and Table II, respectively, keeping the Grass- 
mann character of the products of the two types of polynomials unchanged. This 
can be understood if taking into account that S mh = — ja m a K and that the op- 
erator a m changes the polynomials of an odd Grassmann character of Table I, 
into an even polynomial, transforming a left handed Weyl spinor of one family 
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into a right handed Weyl spinor of another family, while d h changes simultane- 
ously the SU(2) doublet of an even Grassmann character into a singlet of an odd 
Grassmann character. 

The symmetry called the mirror symmetry presented in this approach, is not 
broken, as none of the symmetry is broken. We only have arranged basic states to 
demonstrate possible symmetries. 

We can express the generators of subgroups SU(3) and U(1) of the group 
SO (6) in terms of the generators <S ab (according to Eq.(^)). 

We find (since the indices 9, 10, 11, 12, 13, 14 are reserved for the subgroup 
SO(6)) 

t 61 :=1(S 912 -S 1011 ), t 62 :=1(S 911 + S 1012 ), t 63 := 1(S 9 10 - S ] 1 ,2 ), 

T 64. = 1 (<S 9 14_ <s1 13 )) T 65. = ^ ((S 9 1 3 +5 10 14 )) ^6 ._ 1 ^1 1 14 _ ^12 13 }> 



T*7 : =hs u 13 +«S 1214 ), t 68 :=-L(5 9,0 +5 11 12 -25 1314 ). (34) 
2 2V3 

(35) 

One finds in addition 

t 71 :=-l(5 910 +5 11 12 + 5 1314 ). (36) 

The algebra for the subgroups SU(3) and U(l ) follows from the algebra of 
the Lorentz group SO(1 , 13) 

{ T 6i , t 6 '} = if ijk T 6 \ {t 71 , T 6i } = 0, for each i. (37) 

The coefficients ftjk are the structure constants of the group SU(3). 

We can find the eigenvectors of the Casimirs of the groups SU(3) and U(l ) 
for spinors as polynomials of 9 h , h = 9, 14. The eigenvectors, which are poly- 
nomials of an even Grassmann character, can be found in ref.[pT|. We shall present 
here only not yet published [||] polynomials of an odd Grassmann character. 
Table III: The eigenstates of the operators for spinors t 63 , f 68 , f 71 are presented 
for odd Grassmann character polynomials. We find four triplets, four anti- triplets 
and eight singlets. One sees that complex conjugation transforms one triplet into 
anti-triplet, while f 8x transform triplets into anti-triplets or singlets. 
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One finds four triplets and four anti-triplets as well as eight singlets. Besides 
the eigenvalues of the commuting operators f 63 and f 68 of the group SU(3) also 
the eigenvalue of f 71 forming U(1 ), is presented. The operators f 81 := US 9 12 + 
§ioii )> T 82 . = i(j9ii _s 1012 ), t 83 := 1(S 914 + S 1013 ), t 84 := |(5 9 13 - 
51014^ ^ss := ^(S 11 14 +S 1213 ), T 86 :=l(S n 13 -S 1214 ), which transform 
triplets of the group SU(3) into anti-triplets and singlets with respect to the group 
SU(3). 
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The spinorial representations of the group SO ( 1 , 1 3) are the direct product of 
polynomials of Table I, Table II and Table III. 

We can find all the members of a spinorial multiplet of the group S0(1 , 13) 
by applying S ab on any initial Grassmann odd product of polynomials, if one 
polynomial is taken from Table I, another from Table II and the third from Table 
III. In the same multiplet there are triplets, singlets and anti-triplets with respect 
to SU(3), which are doublets or singlets with respect to SU(2), and are left and 
right handed with respect to SO ( 1 , 3) . 

We can arrange in the same sense also eigenstates of operators of vectorial 
character, with bosonic character. In this paper we shall not do that. 



5.2 Dynamical arrangement of representations of SO(1 , 1 3) with respect to 
subgroups SO(1 , 7) and SO(6) 

To see how Yang-Mills fields enter into the theory, we shall rewrite the Weyl-like 



equation in the presence of the gravitational field (26) in terms of components 
of fields which determine gravitation in the four dimensional subspace and of 
those which determine gravitation in higher dimensions, assuming that the co- 
ordinates of ordinary space with indices higher than four stay compacted to un- 
measurable small dimensions (or can not at all be noticed for some other reason). 
Since Grassmann space only manifests itself through average values of observ- 
ables, compactification of a part of Grassmann space has no meaning. However, 
since parameters of the Lorentz transformations in a freely falling coordinate sys- 
tem for both spaces have to be the same, no transformations to the fifth or higher 
coordinates may occur at measurable energies. Therefore, at low energies, the 
four dimensional subspace of Grassmann space with the generators defining the 
Lorentz group SO (1,3) is (almost) decomposed from the rest of the Grassmann 
space with the generators forming the (compact) group SO ( d — 4), because of the 
decomposition of ordinary space. This is valid on the classical level only. 

According to the previous subsection, the breaking of symmetry of SO ( 1 , 1 3) 
should, however, appears in steps, first through S0(1 , 7) x SO (6) and later to the 
final symmetry, which is needed in the Standard Model for massless particles. 

We shall comment on possible ways of spontaneously broken symmetries by 
studying the Weyl equation in the presence of gravitational fields in d dimensions 
for massless particles (Eqs.^l ^)) 

Y q Poq=0, poa=f%Pcm, Vov.=V\i-^S ab w abl i- (38) 



Standard Model case To make discussions more transparent we shall first com- 
ment on the well known case of the Standard model. Before the breaking of the 
symmetry SU(3) x SU(2) x U(l ) into SU(3) x U(l ), the canonical momentum po« 
( a = 0, 1 , 2, 3 and d = 4) includes the gauge fields, connected with the groups 
SU(3), SU(2) and U(1 ). We shall pay attention on only the groups SU(2) and U(1 ), 
which are involved in the breaking of symmetry 



Vooc = Voc~ 9T x A l « - g'YB a , 



(39) 
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where g and g' are the two coupling constants. Introducing — t 1 ± it 2 , the 
superposition follows A ± a = A 1 a =p i A 2 a . If defining A 3 a = == = Z a + 

, 1 A a and B a = Z a H , g/g A«, so that the transfor- 

VT+ts/g') 2 v /l +(9/9') 2 VT+'s/g') 2 



mation is orthonormalized, one can easily rewrite Eq.([39|) as follows 

Poa =P a - |(T+A+ a + T-A-J + - ? 9g/ QA a + - f = Q'Z X . (40) 
2 V9 2 + 9' 2 vT+F 

with ^ 

Q=t 3 +Y, Q'=t 3 -(^) 2 Y. (41) 

9 

In the Standard Model < Q > is the conserved quantity and < Q ' > is not, 
since < Q > is zero for the Higgs fields in the ground state, while < Q ' > is 
nonzero (< Q' >= -1(1 + (^) 2 ))- 

If no symmetry is spontaneously broken, that is if no Higgs breaks symmetry 
by making a choice for his ground state symmetry, the only thing which has been 
done by introducing linear superpositions of fields, is the rearrangement of fields, 
which always can be done without any consequence, except that it may help to 
better see the symmetries. 

Spontaneously breaking of symmetries causes the non-conservation of quan- 
tum numbers, as well as massive clusters of fields. 



Spin connections and gauge fields leading to the Standard Model We shall 
rewrite the canonical momentum of Eq.(p8|) to manifest possible ways of breaking 
symmetries of SO(1 , 13) down to the symmetries of the Standard model. We first 
write 



y a Voa = = f a f ^aP<v - (y m f a m + y h f a H)Po« + (? m f V + y H f a h )vo., (42) 

with a, m e {0, 1 , 2, 3} and cr, h S {5, 14} to separate the d = 4 dimensional sub- 
space out of d = 14 dimensional space. We may further rearrange the canonical 
momentum po H 

P<V=P^S hlH2 o)K,H 2H -^S k ' k2 cu klk2 ^-ls h ' k ' cu hlW , (43) 

with hi e {0,1, .., 8} and lq € {9, 14} so that S h ' Hz define the algebra of the 
subgroup SO (1,7), while S hlh2 define the algebra of the subgroup SO (6). The 
generators S h ' k ' rotate states of a multiplet of the group SO ( 1 , 1 3) into each other. 

Taking into account subsection |5.1| we may rewrite the generators S ab in 
terms of the corresponding generators of subgroups f Al and accordingly, simi- 
larly to the Standard Model case, introduce new fields (see subsection ^^), which 
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are superpositions of the old ones 

gA 3 V = ^(o> 5 8m - u)67n-)) 9 A 32 m = ^(o>57m + 0> 6 8m), 



^ 1 

g A 3 % = -(cu 56 ^-cu78^), (44) 



gA 4 V = ^(tL>56,x + a> 78lx ), 



g A 51 H = ^(tu 58 ^ + u>67n), g A 52 ^ = ^(cu 5 7n - OJes^)- (45) 
It follows then 

ls hlh ^K,h 2H = gT Ai A A V, (46) 

where for A = 3, i = 1 , 2, 3, for A = 4, i = I and for A = 5 i = 1,2. Accordingly, 
the fields A Al u are the gauge fields of the group SU(2), if A = 3 and of U(1) if 
A = 4. Since f 41 and f 5x form the group SU(2) as well, the corresponding fields 
could be the gauge fields of this group. The breaking of symmetry should make 
a choice between the gauge groups U(1 ) and SU(2). 

We leave the notation for spin connection fields in the case that ht e {0, 1 , 2, 3} 
unchanged. We also leave unchanged the spin connection fields for the case, that 
hi = 0, 1 , 2, 3 and h 2 — 5, 6, 7, 8 as well as for the case, that hi e {0, 1 8} and 
ki € {9, .., 14}, while we arrange terms with kt <E {8, 14} to demonstrate the 
symmetry SU(3) and U(1) 

g A 61 H = ~(o> 9 Un - CU 10 lljt), g A 62 R = ~(o> 9 ,m + O> 10 12ji)» 

I 



g A 63 H = -(o> 9 10(i- t"n 12n). g A 64 H = ~(a> 9 14( 4. - CUiO 13m-). 



g A 65 H = ~(o> 9 + a)io i4n), g A 66 H = -(tun u H - a> 12 i3n)» ( 47 ) 



g A 67 R = -(a)! ! 13^ + 0)12 14n), 



1 

2' 
1 



g A 68 H = ^-Wu>9 10^ + 0)11 12^-2ail3 Hp.) 
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= g A 7 V = -V ao("»9 lOji + tuii i2m + o> 13 hm-) = 9\/ ^A^ 71 = g' A ^ 71 . (48) 



We may accordingly define fields g A 8 ^ = 5(0)9 n H + o>i n^.), g A 82 ,^ = 

2-(o)9 — O>io 12]x)) g A 83 ,^ = 5(0)9 14^ + O>io 13m). 9 A 84 R = 2(0)9 13jl — 
O>1014m). g A 85 H = 2(0)11 14m + 0>12 13m). 9 A 86 H = \ (O) 11 1 3^ - 0>1 2 14m) ' 

so that it follows 

\$^ 2 o) klk2H = g'f Ai A' A V, (49) 

with A = 6, 7, 8 and all A^ Ai = A A1 , except for A = 7, i = 1, which is defined 
in Eq.@. While A<%, i G {1 , .., 8}, form the gauge field of the group SU(3) and 
A ^ corresponds to the gauge group U(1), terms gf 71 A 7l H transform SU(3) 
triplets into singlets and anti-triplets. Again, without additional requirements, all 
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the coupling constants g are equal. To be in agreement with what the Standard 
model needs as an input, we further rearrange the gauge fields belonging to the 
two 11(1 ) fields, one coming from the subgroup SO(1 , 7) the other from the sub- 
group SO (6). We therefore define 



Yi =f 



41 



,71 



Y 2 = — T 



41 



,71 



(50) 



and accordingly similarly to the Standard Model case of subsection 5.2 we make 
the corresponding superpositions of the fields A 41 |x and A' 71 

The rearrangement of fields demonstrates all the symmetries of the massless 
particles of the Standard Model and more. For further comments on the coupling 
constants of the fields before and after the break of symmetries see ref . [ p^| 

Taking into account Tables I, II and III one finds for the quantum numbers of 
spinors, which belong to a multiplet of SO(1 , 7) with left handed SU(2) doublets 
and right handed SU(2) singlets and which are triplets or singlets with respect to 
SU(3), the ones, presented on Table IV. We use the names of the Standard model 
to denote triplets and singlets with respect to SU(3) and SU(2). 
Table IV: Expectation values for the generators f 63 and t of the group SU(3) 
and the generator f 71 of the group 11(1), the two groups are subgroups of the 
group SO (6), and of the generators f 33 of the group SU(2) , f 41 of the group 
11(1 ) and T' 4 ' of the group SO(1 , 3), the three groups are subgroups of the group 
SO(1 , 7), for the multiplet (with respect to SO(1 , 7)), which contains left handed 
(< r (4) >=-1)SU(2) doublets and right handed (< r (4) >= 1) SU(2) singlets. In 
addition, values for Y] and Y 2 are also presented. Index i of Ut , d\ , "Vt and et runs 
over four families presented in Table I. 







SU(2) doublets 

f 33 f 41 f 71 Yt Y 2 r' 4 ' 


SU(2) singlets 

f 33 f 41 f 71 Yt Y 2 r< 4 > 


SU(3) triplets 

t 6 3 = 

(h -h ) 
t 6 8 = 


Ui 

di 


1/2 1/6 1/6 1/6 -1 
-1/2 1/6 1/6 1/6 -1 


1/2 1/6 2/3 -1/3 1 
-1/2 1/6 -1/3 2/3 1 


SU(3) singlets 

f 63 = 
f 68 =0 


"Vi 

ei 


1/2 -1/2-1/2-1/2 -1 
-1/2 -1/2-1/2-1/2 -1 


1/2 -1/2 0-11 
-1/2-1/2 -10 1 



We see that, besides Y2, these are just the quantum numbers needed for mass- 
less fermions of the Standard Model. The value for the additional hyper charge 
Y 2 is nonzero for the right handed neutrinos, as well as for other states, except 
right handed electrons. 

Since no symmetry is broken yet, all the gauge fields are of the same strength. 
To come to the symmetries of massless fields of the Standard Model, surplus sym- 
metries should be broken so that all the fields fields cu a b H which do not determine 
the fields A A \, A = 3, 6 (Eqs.(||||)) and A 41 and A 71 should be invisible at low 
energies. 

The mirror symmetry should also be broken so that multiplets of SO(1 , 7) 
with right handed SU(2) doublets and left handed SU(2) singlets become very 
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massive. All the surplus multiplets, either bosonic or fermionic should become of 
large enough masses not to be measurable at low energies. 

The proposed approach predicts four rather than three families of fermions. 

Although in this paper, we do not discuss possible ways of appearance of 
spontaneously broken symmetries, bringing the symmetries of the group SO ( 1 , 1 3) 
down to symmetries of the Standard model (for these discussions the reader 
should look at refs. [ p^p^ ], which will also appear at the proceedings), we still 
would like to know, whether there are terms in the Weyl equation (Eq.[l2|) which 
may behave like the Yukawa couplings. We see that indeed the term y h f a hVoa, 
with h e {5, 6, 7, 8} and cr g {5,6,..} really may, if operating on a right handed 
SU(2) singlet transform it to a left handed SU(2) doublet. We also can find among 
scalars the terms with quantum numbers of Higgs bosons (which are SU(2) dou- 
blets with respect to operators of the vectorial character.) All this is in preparation 
and not yet finished or fully understood. 

6 CONCLUDING REMARKS 

In this paper, we demonstrated that if assuming that the space has d commuting 
and d anti-commuting coordinates, then, for d > 14, all spins in d dimensions, 
described in the vector space spanned over the space of anti-commuting coordi- 
nates, demonstrate in four dimensional subspace as the spins and all the charges, 
unifying spins and charges of fermions and bosons independently, although the 
super-symmetry, which guarantees the same number of fermions and bosons, is 
a manifesting symmetry. The anti-commuting coordinates can be represented by 
either Grassmann coordinates or by the Kahler differential forms. 

We demonstrated that either our approach or the approach of differential 
forms suggest four families of quarks and leptons, rather than three. We have 
shown that starting (in any of the two approaches) with the Lorentz symmetry 
in the tangent space in d > 14, spins degrees of freedom ( described by dynam- 
ics in the space of anti-commuting coordinates) manifests in four dimensional 
subspace as spins and color, weak and hyper charges, with one additional hyper 
charge, in a way that only left handed weak charge doublets together with right 
handed weak charge singlets appear, if the symmetry is spontaneously broken 
from SO(1, 13) first to SO(1,7) and SO (6), so that a multiplet of SO(1,7) with 
only left handed SU(2) doublets and right handed SU(2) singlets survive, while 
the mirror symmetry is broken, and then to SO (1,3), SU(2) , SU(3) and U( 1 ) . 

We have demonstrated that the gravity in d dimensions manifests as ordi- 
nary gravity and all gauge fields in four-dimensional subspace, after the break- 
ing of symmetry and the accordingly changed coupling constant. We also have 
shown that there are terms in the Weyl equations, which in four-dimensional sub- 
space manifest as Yukawa couplings. 

The two approaches, the Kahler one after the generalization, which we have 
been suggested, and our, lead to the same results. 

A lot of work and ideas are still needed to show that the approach, although 
a very promising one, is showing the right way behind the Standard Model. 
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Semitopological Q-Rings* 

M. Axenides** 

Institute of Nuclear Physics, N.C.S.R. Demokritos, 15310, Athens, Greece 

Abstract. Semitopological Vortices (Q-Rings) are identified to be classical soliton config- 
urations whose stability is attributed to both topological and nontopological charges. We 
discuss some recent work on the simplest possible realization of such a configuration in a 
scalar field theory with an unbroken U ( 1 ) global symmetry. We show that Q-Rings corre- 
spond to local minima of the energy, exhibit numerical solutions of their field configura- 
tions and derive virial theorems demonstrating their stability. 

As we celebrate the 60th birthday of Holger Bech Nielsen we can without 
doubt assess his contributions to the development of the theory of strings and 
vortices to bear the strongest possible impact. Indeed his early work on the de- 
velopment of multiparticle dual models was soon after followed by the 
introduction of the string picture in the study of strong interaction physics [0]. At 
the time it improved tremedously our physical understanding of dual models[||]. 
The string concept, of course, was bound to become much more useful in the uni- 
fication of particle interactions with gravity. Aside from Holger 's contribution 
to the development of the "string idea" he much later provided the first covari- 
ant formulation of a vortex in a theory with spontaneously broken abelian gauge 
symmetry [||. The stability of such a gauged vortex is due the presence of a topo- 
logical charge. The cosmic role of such topological defects in the phase transi- 
tions of the early universe has been important. It is in the spirit of this line work 
of Holger 's that we will present a novel class of vortex like configurations that 
share some of the properties of topological solitons as well as those that are non- 
topological in character.Hence their identification as semitopological. The work 
was done in collaboration with E.G.Floratos,S.Komineas and L.Perivolaropoulos[^] 

Non-topological solitons (Q balls) are localized time dependent field con- 
figurations with a rotating internal phase and their stability is due to the con- 
servation of a Noether charge Q[^]. They have been studied extensively in the 
literature in one, two and three dimensions^. In three dimensions, the only lo- 
calized, stable configurations of this type have been assumed to be of spherical 
symmetry hence the name Q balls. The generalization of two dimensional (pla- 
nar) Q balls to three dimensional Q strings leads to loops which are unstable due 
to tension. Closed strings of this type are naturally produced during the collisions 
of spherical Q balls and have been seen to be unstable towards collapse due to 
their tension[|oKll]]. 

* Contributed to the 4th Workshop "What Comes Beyond The Standard Model", 

Bled,Slovenia, July 17-27 2001 in honor of Holger Bech Nielsen's 60thBirthday. 
** e-mail : axenides@mail. demokritos . gr 
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There is a simple mechanism however that can stabilize these closed loops. It 
is based on the introduction of an additional phase on the scalar field that twists 
by 27tN as the length of the loop is scanned. This phase introduces additional 
pressure terms in the energy that can balance the tension and lead to a stabi- 
lized configuration, the Q ring. This type of pressure is analogous to the pressure 
of the superconducting string loops[Q (also called 'springs'^]). In fact it will 
be shown that Q rings carry both Noether charge and Noether current and in 
that sense they are also superconducting. However they also differ in many ways 
from superconducting strings. Q rings do not carry two topological invariants 
like superconducting strings but only one: the winding N of the phase along the 
Q ring. Their metastability is due not only to the topological twist conservation 
but also due to the conservation of the Noether charge as in the case of ordinary 
Q balls. Due to this combination of topological with non-topological invariants 
Q rings may be viewed as semitopological defects. In what follows we demon- 
strate the existence and metastability of Q rings in the context of a simple model. 
We use the term 'metastability' instead of 'stability' because^/imYe size fluctuations 
can lead to violation of cylindrical symmetry and decay of a Q ring to a Q ball as 
demonstrated by our numerical simulations. 

Consider a complex scalar field whose dynamics is determined by the La- 
grangian 

£ = -9^0*3^0-11(101) (1) 

The model has a global U(l) symmetry and the associated conserved Noether 
current is 

= Ira(0*3 H 0) (2) 

with conserved Noether charge Q = J d 3 x Jo- Provided that the potential of (jj]) 
satisfies certain conditions [^,||| the model accepts stable Q ball solutions which 
are described by the ansatz O = f (r)e lwt . The energy density of this Q ball con- 
figuration is localized and spherically symmetric. The stability is due to the con- 
served charge Q. 

In addition to the Q ball there are other similar stable configurations with 
cylindrical or planar symmetry but infinite, not localized energy in three dimen- 
sions. For example an infinite stable Q string that extends along the z axis is de- 
scribed by the ansatz 

O =n P )e iwt (3) 

where p is the azimouthal radius. This configuration has also been called 'planar' 
or 'two dimensional' Q ball[||. 

The energy of this configuration can be made finite and localized in three di- 
mensions by considering closed Q strings. These configurations which have been 
shown to be produced during spherical Q ball collisions JTo| , pl| ] are unstable to- 
wards collapse due to their tension. In order to stabilize them we need a pressure 
term that will balance the effects of tension. This term appears if we substitute 
the string ansatz (^) by the ansatz of the form 



O = f(p)e itu V a(z) 



(4) 
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where <x(z) is a phase that varies uniformly along the z axis. This phase intro- 
duces a new non-zero T z component to the conserved current density (^). The 
corresponding current is of the form 



. dcx 

dz — — Z7T 

dz 



dppf 2 



(5) 



Consider now closing the infinite Q string ansatz (|J) to a finite (but large) loop of 
size L. The energy of this configuration may be approximated by 



E = 



47tL J dp p f 2 

(27tN) 2 7T 



+ 7tL 



dp p f 



a 



L 



dp p f z + 2ttL 



dp pU(f) 



Ii + Iz + I3 + 14 



where we have assumed a(z) = ^^-z and the terms h are all positive. Also Q is 
the charge conserved in 3D defined as 

Q = cu2ttL 



dppf 2 



(6) 



The winding 27tN = J" dz is topologically conserved and therefore the current 
(||) is very similar to the current of superconducting strings. 

After a rescaling p — > a/ATP/ z — > A2Z the rescaled energy may be written 



as 



E = 



1 



-Ii + A 2 I 2 



Ai 



I3 + Ai A2I4 



(7) 



A1A2 A2 

This configuration can be metastable towards collapse since Derrick's theorem[fj"4|| 
is evaded due to the time dependence Jl5|Jl6| ] of the configuration Demanding 
metastability towards collapse in any direction we obtain the virial conditions 



I3 + U = Ii 

I 2 + I4 = Ii + I 3 



(8) 
(9) 



In order to check the validity of these conditions numerically we must first 
solve the ode which f obeys. This is of the form 



f" + -f + (w 2 - (27rN) 2 /L 2 )f - U'(f) = 



(10) 



with boundary conditions f(oo) = and ^(0) = 0. Equation dTfj| ) is identical 
with the corresponding equation for 2D Qballs[[l6| (see ansatz (^) with the re- 
placement of cu 2 by 

2 ^ N)2 .c' 2 (11) 



L 2 



Solutions of © for various cu' and U(f ) = |f 2 - If 3 + f f 4 with B = 4/9 were 
obtained in Ref . p6[ ] . Now it is easy to see that the first virial condition (^) may 
be written as 



a 



dp p f 2 = 2 



dp pU(f ) 



(12) 
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This is exactly the virial theorem for 2D Qballs (infinite Q strings) with N = 
and field ansatz given by @ with a> replaced by cu'. The validity of this virial 
condition has been verified in Ref . Jl6| ] . This therefore is an effective verification 
of our first virial condition (||). 

The second virial condition (0) can be written (using the first virial (^)) as 

2h = h (13) 

which implies that 



2tiN 2 _ J dp pf 



a 



L 2 Jdppf 2 (M) 

This can be viewed as a relation determining the value of L required for balancing 
the tension ie for metastability 

These virial conditions can be used to lead to a determination of the energy 

as 

E = 2(I 1 +I 3 ) (15) 

In the thin wall limit where 2n J" dppf 2 = Af 2 , (A is the surface of a cross section 
of the Q ring) this may be written as 

c ^ Q 2 , (27rN) 2 Af 2 (16) 



2LAf 2 2L 



and can be minimized with respect to f 2 ,. The value of fo that minimizes the en- 
ergy in the thin wall approximation is 



2tiNA 

Substituting this value back on the expression (Q for the energy we obtain 



(17) 



E = ^2 (18) 

This is consistent with the corresponding relation for spherical Q balls which in 
the thin wall approximation lead to a linear increase of the energy with Q . 

The above virial conditions demonstrate the persistance of the Q ring con- 
figuration towards shrinking or expansion in the two periodic directions of the Q 
ring torus for large radius. In order to study the Q rings of any size and its stabil- 
ity properties towards any type of fluctuation we must study the full evolution 
of a Q ring in 3D by performing energy minimization and numerical simulation 
of evolution. This is precisely what we did for a potential energy given by : 

U(cM = ^|0| 2 - 1|(D| 3 + ||(D| 4 (19) 
The ansatz we used that captures the above mentioned properties of the Q ring is 

O =f(p,z) e i[tut+TltH (20) 
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where the center of the coordinate system now is in the center of the torus that 
describes the Q ring and the ansatz is valid for any radius of the Q ring. We have 
also replaced N by n. 

The energy of this configuration is 



1 Q 2 1 

2J?2qV + 2 



2-f2 







dV + 













U(f ) dV 



dV 



The field equation for O is 

O - AO + ® - |(D|® + B|(D| 2 cI) = 

Substituting the ansatz ( p0| ) we find that f (p, z) should satisfy 

3 2 f 1 3f n 2 f 3 2 f , 2 r2 „ 

- T + -- r + - T + (cu 2 -l)f + f 2 -Bf 3 =0 

op z p op p z 3z z 



(21) 



(22) 



(23) 



In order to solve this equation we minimized the energy ( |2l| ) at fixed Q using the 
algorithm 



3f _ _6E _ 
9t ~ 5f ^ 



(24) 



3f 
3r 



n 2 f 3 2 f 
+ ■ 



with boundary conditions f(0,z 
rithm is checked by 



3 2 f 1 3f 

3p 2 p 3p P 2 3z 
+ (cu 2 -l)f + f 2 -Bf 3 

3f(p,z) | 



0, 



(25) 

0. The validity of the algo- 



dG/dT = 6E/6f df/dx: 
In (^i|) w is defined as 



(dE/dx) 2 < 



w 



Jf 2 dV 

In the algorithm, we have used the initial ansatz: 



f(p,z) = const exp const 



(26) 
(27) 

(28) 



where po is a fixed initial radius. The energy minimization resulted to a non- 
trivial configuration f (p, z) for a given set of parameters B, n, Q in the expression 
for the energy. We then used (27) to calculate cu and constructed the full Q ring 
configuration using (p0|). While the details of our numerical analysis can be found 
elsewhere we just report the main results. 

It was verified that the Q ring configurations evolve with practically no dis- 
tortion and are metastable despite their long evolution. Finite size nonsymmetric 
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fluctuations were found to lead to a break up and eventual decay of the Q ring to 
one or more Q balls. Thus a Q ring is a metastable as opposed to a stable config- 
uration. 

The Q ring configuration we have discovered is the simplest metastable ring- 
like defect known so far. Previous attempts to construct metastable ring-like con- 
figurations were based on pure topological arguments (Hopf maps) and required 
gauge fields to evade Derrick's theorem due to their static nature [ p^ , ^8| . This 
resulted in complicated models that were difficult to study analytically or even 
numerically. Q rings require only a single complex scalar field and they appear 
in all theories that admit stable Q balls including the minimal supersymmetric 
standard model (MSSM). The simplicity of the theory despite the non- trivial ge- 
ometry of the field configuration is due to the combination of topological with 
non-topological charges that combine to secure metastability without added field 
complications. 

The derivation of metastability of this configuration opens up several in- 
teresting issues that deserve detailed investigation. They pertain to the various 
mechanisms of formation of Q Rings (Kibble and Affleck-Dine mechanisms, Q 
ball collisions etc.) as well as on the dependence of the winding N on Q. We hope 
to have something interesting to report in the forthcoming 70th birthday celebra- 
tion of Holger. 
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Non-local Axial Anomalies in the Standard Model 

D. Melikhov and B. Stech 

Institut fur Theoretische Physik, Universitat Heidelberg, Philosophenweg 16, D-69120, 
Heidelberg, Germany 

Abstract. We demonstrate that the amplitude (pY|9^(qYvY5q)|0) does not vanish in the 
limit of zero quark masses. This represents a new kind of violation of the classical equa- 
tion of motion for the axial current and should be interpreted as the axial anomaly for 
bound states. The anomaly emerges in spite of the fact that the one loop integrals are 
ultraviolet-finite as guaranteed by the finite-size of bound-state wave functions. As a re- 
sult, the amplitude behaves like - 1 /p 2 in the limit of a large momentum p of the current. 
This is to be compared with the amplitude (YY|9~"(qYvY5q}|0) which remains finite in the 
limit p 2 — > oo. 

The observed effect requires the modification of the classical equation of motion of 
the axial-vector current by non-local operators. The non-local axial anomaly is a general 
phenomenon which is effective for axial-vector currents interacting with spin-1 bound 
states. 

1 Introduction 

It is well-known that the classical equations of motion can be violated in quan- 
tum field theories. One of the best-known examples of such a violation is the axial 
vector current: in the limit of vanishing quark masses the divergence of the neu- 
tral hadronic axial vector current should be zero but in fact is proportional to a 
local operator containing the photon fields. This is the famous Adler-Bell-Jackiw 
anomaly The standard model is free of these local anomalies: the anomalies of 
the hadronic (quark) part of the divergence is compensated by the leptonic part. 
Such cancelations are in fact required for any renormalizable field theory. 

Also the amplitude of the decay n° — > yy, the process which led to the 
discovery of the axial anomaly, is free of anomalies! Only when one relates the 
7T° — > yy decay amplitude via PCAC to matrix elements of the pure hadronic part 
of the axial-vector current, one is able to use the above mentioned local anomaly 
term for an estimate of the decay rate. 

In this article we will show that the standard model has non-local anomalies 
and that these anomalies are the only ones which do not cancel. They arise from 
the coupling of the of the local axial vector current to two vector particles where 
at least one of them is a bound state. Examples are the couplings of the axial vec- 
tor current to a photon and a vector meson, or to two vector bound states. For the 
leptonic axial current an example is the coupling to a photon and orthopositron- 
ium, or two orthopositronium states. 

We will demonstrate the existence of such anomalies for the case where the 
isovector axial current emits a photon and a p-meson in the limit m q — > 0. Since 
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this process necessitates a quark loop, it is evident that there is no cancellation 
by the leptonic part of the axial current. This anomaly turns out to be non local, 
and thus has no consequences for the renormalizability of the standard model. 
An example for a physical process for which such an anomaly is important is 
weak annihilation in radiative B — > py decays: The B meson is annihilated by the 
uysYvb axial current, and the dysYyU part of the axial current for light quarks 
generates the photon and the p meson. 

It is interesting that these anomalies had so far not been discovered even 
though their derivation is not difficult. There are many examples in the literature 
in which - in the limit of massless quarks - the corresponding divergence of the 
axial current has been put equal to zero by relying on the validity of the classical 
equation of motion. 

This article follows our recent analysis presented in Ref . [|| . 

2 The Adler-Bell-Jackiw anomaly 

The analysis of two-photon decays of pseudoscalar mesons in the late 60-s led to 
the discovery of the famous axial anomaly the divergence of the axial vector 
current violates the classical equation of motion and does not vanish in the limit 
of zero fermion masses. For a quark of mass m q and charge eQ q the properly 
modified equation of motion contains a local anomalous term and has the form[] 

9"(qy.y 5 q) =2im q qY 5 q + N c ^]!FF, F^ = e^ a(3 F a P. (1) 

with F^ the electromagnetic field tensor. This modification of the classical equa- 
tion of motion accounts for the fact that the form factor G Y defined by the 2- 
photon matrix element 

(y(qi)T(q2)|3"(qy.y 5 q)|0) = e 2 e q , ^^G^p 2 , qf, q 2 ), (2) 

does not vanish for m q = but turns out to be a constant independent of the 
current momentum p = qi + q 2 : 

G Y (p 2 ,q 2 = q 2 =0) =-2N c (Q q ) 2 /47t 2 . (3) 

In this letter we study the properties of axial currents when one of the photons, 
y(q2)/ is replaced by a vector meson V(q2), e. g. a p-meson. We demonstrate that 
the form factor G v defined according to the relation 

<y(qi )V(q 2 )|3-(q'y.y 5 q)|0) = ee q , etq2e .G v (p 2 , qf , q 2 ) (4) 

has also an anomalous behavior and does not vanish for massless quarks. This 
occurs in spite of the fact that the vector meson is a bound q q state and the cor- 
responding loop graph has no ultraviolet divergence. Moreover, the anomalous 
behavior is observed for both, the neutral and the charged axial-vector currents. 

1 We use the following notations: e = y*4mx e m, y 5 = VY < V'y 2 Y 3 / e 0123 = — 1, 
Sp(y 5 yVyV) = 4ie^ a|3 ,F^ = 3^A. - 9^. e a hcd = e a |3^a a b p c^ for 
any 4-vectors a, b, c, d. 
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The classical equation of motion reads 

^(qVvKsq) =i(m.q + m' q )q'y 5 q + e(Q q > - QqJA^q'y^Tsq, (5) 

where A v is the electromagnetic field. The (yV|...|0) matrix element of the second 
term on the r.h.s. of O vanishes to order e. Therefore, to this order, the classical 
equation of motion (H) predicts G v = for m q ' = m q = 0. We find however to 
order e and for large p 2 

G v ~M v fv/p 2 , (6) 

where My and f v denote the mass and the decay constant of the vector meson. 
Because of the dependence on p 2 , the newly found deviation from the classical 
equation of motion for bound states corresponds to a non-local anomaly. 

We are interested in the region |p 2 | ^> m 2 , Aq CD , in which case the quarks 
in the triangle diagram have high momenta and their propagation can be treated 
perturbatively. We discuss in parallel the yy and y V final states in order to show 
how the anomaly emerges in both cases. As in Ref. [g], we consider the spectral 
representation for the axial-vector current itself before forming the divergence. In 
distinction to where the spectral representation in p 2 was considered, we use 
the spectral representation in the variable q|- This allows us to take bound state 
properties into account. 

3 The absorptive part of the triangle amplitude 

The amplitude of the single-flavor axial current between the vacuum and the 
two-photon and the photon-vector meson states, respectively, can be written in 
the form 

e*P(q 2 )e* a (q 1 )T vaP (q 1) q 2 ). (7) 

The absorptive part t^ a p of Ty^p is calculated by setting the two quarks attached 
to the external particle with the momentum q 2 on the mass shell, see Fig 1. 




Fig.l. Diagrams describing (y(qi )y(q 2 )|qy^ysq|0) and/or (y(qi )V(q2)|qy^y5q|0), with 
p = q i + q 2 . The cut corresponds to the calculation of the absorptive part in the variable 
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is our basis for the spectral representation of T-v^p in terms of the variable q^. The 
coupling at the vertex |3 is y p e Q q if the particle 2 is a photon, and — y p g ( q§ )/ \/N7 
if it is a vector meson.^The coupling g ( q| ) will be further discussed below. By 
taking the trace and performing the integration over the internal momentum in 
the loop it is straighforward to obtain ty a p. The result is automatically gauge- 
invariant 

qftvap(th,c|2) = 0, qfWpCqi^) =0. (8) 

It is therefore possible to write the covariant decomposition of tv a p(qi , q2) in 
terms of three invariant amplitudes 

tv«p(qi , c\z) = -p-^e a pq,q 2 ico + (qf e A , ocP q 2 - qi«e v q 1 p q2 )ici (9) 

+ (q2 e ^(3aq, - q2(3e-vq 2 a qi )lC 2 . (10) 

This Lorentz structure is chosen in such a way that no kinematical singularities 
appear. We take y(qi ) to be a real photon, q 2 =0. Hence, the term containing 
the invariant amplitude Ci does not contribute to the divergence of the current. 
Setting in addition m q = one obtains for Co and C2 with s = q| 



,2 



■ , C(s) s 2 C(s) p 2 

C p,S= — - j-j, C 2 (p ,S = — j-j, (11) 

4n (s — p z J z 47i (s — p z ' z 



where £(s) = 2N c Qq0(s) for the yy process and C(s) = — 2 v / N7Qqg(s)9(s) for 
the y V process. 

Clearly, the absorptive part tvap(qi > q2) of the axial-vector current matrix 
element respects the classical equation of motion, that is 

p 2 Co(p 2 ,s)-sc 2 (p 2 ,s) =0. (12) 



4 The triangle amplitude and its divergence 

The full amplitude T-v a p (qi , q 2 ) has the same Lorentz structure as its absorptive 
part 

T-vap (qi > q2) = — P^eapq, q 2 ^o + (q 2 e-v a p q2 — qiae-vq, pq 2 )lCl (13) 

+ (q^-vpaq, - q2pe-vq 2aqi )iC 2 . (14) 

The absence of any contact terms in Tv a p can be verified by reducing out one 
of the photons and using the conservation of the electromagnetic current. The 
invariant amplitudes Ct can be represented by the following dispersion integrals 



Ci(p 2 ,q 2 =0,q 2 ) = l 

7t 



-2 — ^ds. (15) 

o s — q| - uO 

i 



The full vertex has the form [|] T p = -1/VNJyp - ^ 2m (Id -k 2 )p/( v / s+2m)]g((k 1 + 
k2) 2 ), but the term proportional to (ki — k2)p does not contribute to the trace. The overall 
(— ) sign is the standard choice of the phase of the vector meson wave function which 
leads to a positive leptonic decay constant. 
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All the integrals converge and thus need no subtraction. 
Taking the divergence of Tv a p we find 



Co(P ,S) 2 

— ds - q 2 



o s - q 2 

The form factor G defined in Eqs. (||) and 



C2(P 2 ,S) 

s - q\ 



ds > £q , aq 2 (3 ■ 



now reads 

as 



(s -p 



2U 



ds 



(16) 



(17) 



In the case of the yy process £ (s) is a constant. The integral can be performed and 
gives the well-known value shown in Eq (||). In the case of the yV matrix element 
the integrals converge even better since g(s) which appears in £(s) descibes the 
spatial size of the vector meson. We conclude from Eq. ( p^ ) that the divergence of 
the axial-vector current is nonzero for m q = not only for the yy but also for the 
yV final state! Namely, 



4n 2 



9(s) 



■,212 



ds. 



(18) 



The behavior with respect to p 2 is however different from the yy case and has the 
form G v (p 2 ) - 1 /p 2 for the large values of p 2 where our formula applies. 

For the transition to the y p (isospin-1) arising from the isovector axial current 
we obtain 



G p 



4ti 2 f p M p 



g(s) 



2U 



ds. 



o is -p 
and |p 2 



(19) 



00. fr 



1S 



The parameter k in this equation is non zero for m q 
defined by the relation (p~|dy-vu|0) = f p M p e*. 

Eq. ( |19|) takes into account the soft contribution to the form factor G p . For 
large |p 2 | one should take care of the QCD evolution of the p-meson wave func- 
tion from the soft scale [i 2 - 1 GeV 2 to the scale [i 2 - |p 2 |.[| This can be done most 
directly by expressing k in terms of the p-meson light-cone distribution ampli- 
tudes [0] 



(p(q 2 )|d(x)TAu(0)|0) = -iq 2 A(e*x)f p Mp 



due 1 



e^fpMp 



due 1 



[v) t 



(20) 



(p(q 2 )|a(x)y A y 5 u(0)|0) = --e Mpa e*^q 9 2 x°f P M P 



due 



xuq2X 



(a) 



U 



(21) 



3 We want to point here to the similarity of the form factor G p with the ny transition form 
factor Vrtyy* (p 2 ). For a detailed analysis of the latter we refer to Ref. [||]. Likewise, the 
form factor G pp describing the amplitude (pp|9 v (qyvY5q)|0) has some common feature 
with the pion elastic form factor. 
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The diagrams of Fig 2 explain the procedure we follow: The quark propagator 
connects the axial current with the photon and the two distant space time points 
are bridged by the p meson wave function. This leads to the following expression 
for k 



y y 




Fig. 2. The diagrams for the matrix element (y(c|i)p (q2)|dy^ysu|0), with two space time 
points bridged by the p meson light-cone distribution amplitude. 



k = — 



du 



1 +u 



"9j 



1 



4u 2 ^ u aj - 1 - u 
The leading-twist distribution amplitudes 



(22) 



9 



6u(l — u), 
3 
4 
3 
2 



(1 +(2u-1) 2 ), 
u(l -u)(2u-l) 



(23) 



give the main contribution for large p 2 and lead to the value k = —3/2. Correc- 
tions to this value are calculable in terms of the higher twist distribution ampli- 
tudes [§]. 

Summing up, our results are as follows: 

1. The divergence of the axial vector current and thus the form factor G does not 
vanish in the limit m q — > 0. This holds for the yp final state as well as for the 
yy final state. The observed effect for vector mesons requires a proper mod- 
ification of the equation of motion for the axial current. For large momenta 
p of the axial-vector current the corresponding 'bound state anomaly' can be 
described in terms of a non-local operator appearing at order e (there are no 
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local operators of the appropriate dimension): 

Single flavor current : 

9 (qYVYsq) = 2tm qy 5 q + — — FF 
+ eKQ q D- 1 {3^qy^q) -F^} + 0(D- 2 ) 

Isovector charged current : 

3^(dY-vY5u) = i(ra u + ra d )uy 5 d - ie(Q u - Qd)uyvY 5 d • 

+ eK (Qu + Qd) D _i { 9 ^ ({lY ^ d) . F^} + 0(D- 2 ). (26) 

k can be expanded in a power series of a s . In leading order one obtains the 
value k = —3/2. 

As pointed out in the py anomaly is important in rare decays: for instance, 
in B — > py decays, it substantially corrects the weak annihilation amplitude 
which carries the CP violating phase. 

2. The amplitude for the pp final state also stays finite for ra q = 0. The corre- 
sponding non-local anomalous term for the divergence of the isovector axial 
current 3^, ( q Yv Y5T a q ) appears already at order 0(e). The operator structure 
of the anomalous term is more complicated in this case. One of the possible 
lowest-dimension operators which has a nonvanishing (pp|...|0) matrix ele- 
ment and thus will contribute to the anomalous term (in leading order in 
1 /p ) is the product of the two isovector tensor currents 

e abc e ^«p n -i ( q<VvT * q . qo- apT Cq) . (27) 

Accordingly, for large |p 2 | the (pp|...|0) amplitude of the divergence of the 
axial-vector current is given by the factorized matrix element of the anoma- 
lous term and has a 1 /p 2 suppression 

(P(qi)p(q2)|3"(dy.y 5 u)|0) = e q , e , qi£l G™, 

Gpp ~ f 2 /p 2 + 0(1/p 4 ). (ZO) 

3. We have illustrated the appearance of non-local anomaly due to vector mesons 
in QCD. This anomaly is of a general nature and should be present in any 
theory containing J p = 1 ~~ bound states. For example, the anomaly will also 
contribute to the generation of orthopositronium by the leptonic axial-vector 
current. In contrast to the conventional local anomalies, there is no obvious 
cancellation of the non-local anomalies in the standard model. 
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Could there be a Fourth Generation?* 

CD. Froggatt** and J.E. Dubicki*** 

Department of Physics and Astronomy Glasgow University Glasgow G12 8QQ, Scotland 

Abstract. We investigate the possibility of incorporating a chiral fourth-generation into 
a GUT model. We find that in order to do so, precision fits to electroweak observables 
demand the introduction of light (< M z ) supersymmetric particles. This also enables us 
to provide decay channels to the fourth-generation quarks. Perturbative consistency sets 
an upper bound on the coloured supersymmetric spectrum. The mass of the lightest Higgs 
boson is calculated and found to be above the present experimental lower limit. 

1 Introduction 

Despite the success of the Standard Model (SM), it is far from complete. For 
instance, the origin of the inter-generation mass hierarchy and issues such as 
baryon asymmetry are yet to be resolved. We would expect new physics to con- 
tribute to these areas. At the more fundamental level, nobody knows of any deep 
explanation as to why there should only be three generations of quarks & leptons. 
For these reasons, and others, extensions of the three-generation SM are being 
investigated. We investigate the possibility of incorporating a fourth-generation 
into a GUT model. 

As is well known, all fourth-generation models must adhere to certain experi- 
mental constraints, the first of which stems from precise measurements of the de- 
cay characteristics of the Z-boson performed at LEP. This has set a lower bound 
of Mp > on any non-SM particles that couple to the Z-boson. Ignoring for the 
moment the unnatural hierarchy emerging within the neutrino sector, we assume 
a Dirac mass ~ (^2 ) ^ or * ne heavy neutral lepton. We label the fourth-generation 
doublets explicitly as, 




Physics beyond the Standard Model is severely constrained by precision elec- 
troweak data[[ Assuming SM contributions, fits to LEP data give the radiative 

* We would like to thank David Sutherland for discussions and useful remarks. 

E-mail: c.froggatt@physics.gla.ac.uk 

E-mail: j.dubicki@physics.gla.ac.uk 
1 We assume |V tb | 2 , 1 V TB | 2 - 1 and | V T b | 2 , | V tB | 2 < 1 so that contributions to the Z -> bb 

decay can be ignored. 
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correction parameteij^ S = —0.04 ± 0.1 1 . For a heavy (3> Mz) degenerate fourth- 
generation we obtain AS = and so is ruled out at 99.8% C.L. However, Maltoni 
et al. [Q] have shown that particles with mass M ~ M^ give drastically different 
contributions to S. It is possible for a heavy neutrino N with mass M N - Mi to 
cancel the contributions from the heavy T, B , E with the SM solution 

Me > M N ; M t > M B 

(M E -M N )~3(M T -M B )-60GeV (1) 

and the N must be relatively stable to avoid detection (i.e. mixing matrix elements 
^Ne.ix.T < 10~ 6 ). An extra generation can be accomodated below the 1 cr level (or 
even two generations at 1 .5o~)[| However, it was shown by Gunion et al. ^ that 
the fourth-generation charged lepton must be relatively light (< 60 GeV) in order 
to stay in the perturbative regime below the GUT scale. Although they worked 
within a supersymmetric framework, similar results are expected to hold in the 
SM Decreasing M E below ~ 60 GeV to achieve perturbative unification is not 
consistent with the above SM fits (Eq.(Q)) to the precision data. Therefore, we 
require the cancellations to arise from another sector and so we consider a super- 
symmetric theory. Indeed, it is shown in [Q] that light supersymmetric particles 
can also effect the fit to precision data, allowing for a fourth-generation below the 
2cr level. In particular, neutralinos (x?) and charginos (%f ) with masses < Mz can 
provide the correct sign contributions whilst being consistent with current exper- 
imental limits if nearly degenerate [^] . 

In this paper we investigate the possibility of consistently incorporating a fourth- 
generation into a N = 1 R-parity conserving supergravity model. We assume a 
structure akin to that of the minimal supersymmetric standard model (MSSM3), 
adding a complete chiral fourth-generation and its associated SUSY partners (the 
so called MSSM4). Specifically we require (i) perturbative values for all Yukawa 
couplings at energies up to the GUT scale and (ii) gauge coupling unification. 
These two constraints will be termed collectively as perturbative unification. 

Having satisfied precision data fits, it remains for us to provide solutions that 
(i) unify perturbatively at the GUT scale and (ii) evade the direct experimental 
searches performed at CDF. In section^ we investigate the specific decay chan- 
nels of the fourth-generation quarks and ensure we can provide consistency with 
the experimental direct searches. The solution we present requires the introduc- 
tion of light SUSY particles (M^o , M^± , M B ), so the two-body decays T — > Bx + 
and B — > Bx° are kinematically allowed. These will always dominate over the 

2 S is the well known radiative correction parameter (weak isospin symmetric), normal- 
ized to zero for the SM with M H i fl BS = 100 GeV. 

3 Updated fits to recent LEP data show that an extra generation with Mn = Mu ~ 180 
GeV; Md = Me - 130 GeV gives a x 2 fit for the fourth-generation as good as that of 
the three-generation case [0]. Perturbative unification is not consistent with such large 
masses. 
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one-loop FCNC decays B — > bZ° and two-generation decays B — > cW~ that tra- 
ditional searches have looked for, whilst at the same time suppress the decay rate 
T — > bW + . We note a light, degenerate chargino/neutralino pair is just what is 
needed to provide the necessary cancellations in precision data. 

In section [3]we discuss the influence of specific fourth-generation/ SUSY masses 
on precision data fits in more detail and point out new features that appear when 
trying to satisfy all constraints. In section |] we proceed with a renormalisation 
group study of the four generation MSSM. Under the assumption of a common 
mass scale for the coloured/weak sparticle spectrum (Mcoi/Mwic respectively) 
we will derive upper bounds on M co i. This bound stems from the fact that if we 
are to achieve perturbative values for the top quark Yukawa coupling constant 
h t up to the GUT scale, then the effects of the coloured sparticle spectrum must 
be included in the running of the strong coupling 0.3 at an early stage. Finally, in 
section [| we investigate the lightest Higgs mass in the MSSM4. 



2 Experiment and a Fourth Generation 

In this study we assume the hierarchy within the fourth generation is such that 
Mj > Mb and Me > Mm. We begin our discussion by considering the lep- 
tonic sector. Under the assumption that the mixing between (E, N) and the first 
three generations is negligible, the decay E — > NW* will be dominant. In order 
to evade experimental bounds, the mass difference AMl = Me — Mn must be 
small enough to result in a virtual W* whose decay products are too soft to be 
triggered. Regarding the heavy neutrino, LEPII has set the bound Mm > 70 - 80 
GeV based on the search for N — > IW* (I = e, u. or t) where the mixing matrix el- 
ements Vncix.t > 1 0~ 6 [^]. A heavy neutrino with this mass is inconsistent with 
perturbative unification. However, we have previously assumed that the mixing 
of the fourth generation leptonic sector is in fact negligible (Vn e ^ a < 10~ 6 ) so 
the neutrino is stable enough to leave the detector and in this case only the DEL- 
PHI bound of Mm > 45 GeV from measurements of the Z-width applies. From 
figure (4) in [^] we can see that Me - Mm - 50 GeV is allowed. 

Turning to the quark sector, the requirement of perturbative unification places 
strict upper limits on the masses of the T and B and they certainly must be be- 
low the top quark whose mass is M t = 174± 5.1 GeV. Experimental searches 
for the fourth-generation quarks have mainly concentrated on the B-quark where 
CDF have set a bound of Mb > 199 GeV assuming the branching ratio BR(B — > 
bZ°) ~ 1 . The search is also sensitive to other decay modes; for instance the de- 
cay B — > bh° or B — > cW~ is triggered as long as BR(R — > bZ) is not negligible, 
since the hadronic decay of the h° or W are kinematically similar to those of 
the Z. One possible escape might come about if we note that h° — > NN would 
compete with h° — > bb for Mj^ - M^ (dominating for larger M° ) and would 
provide an invisible signature. For Mb ~ 100 GeV and M^ ~ M^ we would 
also expect BR(B -> bh°) - BR(B -> bZ°) [§. However, CDF still exclude 
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a B-mass in the range 104 GeV -» 152 GeV assuming BR(B -» bZ°) > 50% 
and no sensitivity to the other decay modes Taking M h o ~ 100 GeV (which 
places Mb ~ 105 GeV for B — > bh° to be kinematically allowed) we obtain 
BR(h° — > invisible) - 80% — 90% and so we can expect to rule this possibil- 
ity out. Moreover, we have yet to account for the T decays which turn out to be 
highly constrained if we assume SM-like processes. The channel T — > bW + is 
prohibited for the obvious reason that the T-quark would have been picked up in 
the CDF search for the top quark. Although we might assume that T — > BW* is 
dominant by suppressing Vjb we must notice that, since the TT production cross- 
section is similar to that of BB, we would effectively expect double the event rate 
on the B-quark search. Taking this into account would further strengthen existing 
bounds on the B-quark mass. 

From the ideas presented so far we might conclude that the fourth-generation 
with perturbative unification is not consistent with experimental bounds on the 
(T, B) quarks. We have not, however, considered the possibility of light SUSY par- 
ticles providing decay channels for (T B). In this situation one can constrain the 
light (i.e. < Mz) neutralino/ chargino pair, which is already required by fits to 
precision data, in order to provide the following two-body decays 

T^Bxf; B^Bx? (2) 

Ensuring both decay channels are kinematically accessible, combined with the 
fact that perturbative unification requires^for M t = 170 GeV) Mq < 110 GeV, 
where Mq = Mb — Mj, places severe restrictions on the allowed spectrum. To 
be definite we choose 

M*~50GeV Mq ~ 105— 107 GeV 

M B ~55GeV M t = 170 GeV 

where under the assumption of a nearly degenerate neutralino /chargino pair we 
define the notation AMy = M^± — M^o; M* = M^o ~ M~±. Although these 

x Xi Xi' x Xi Xi ° 

masses seem to be contrived, we note that they can be obtained from reasonable 
assumptions about the supersymmetric sector. In order to obtain AM^ ~ M n + we 
require the hierarchy Mi > M2 where \\±\ is the Higgs mixing parameter 

and is fixed by the requirement of radiative electroweak breaking. Mi and M2 
are the electroweak gaugino masses (for a review of supersymmetry see U). This 
structure can occur naturally when the gaugino masses are dominated by loop 
corrections, originating from superstring models {i.e. the O — II model [|10||). We 
also require substantial Bl/Br mixing to induce the light B mass^. 

4 As we increase M t within the allowed experimental range, the solution becomes harder 
to maintain. Low values of M t are consistent with the fit from the jet + lepton channel 
which is thought to be more precise. 

5 The most recent LEP search for B squarks from B — > bx° decays is insensitive to 
AMg^ = Mg — M.<j < 8 GeV, especially in the case of large mixing [TL2fl. 
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3 Precision Measurements and a Fourth Generation 

It is difficult to provide bounds from precision data without a fully consistent 
study taking into account exact particle masses, any light SUSY spectra present 
and mixings between different flavours. However, it is pointed out in [Q] that a 
highly degenerate neutralino/chargino pair (AM^ - M n +; M^ ~ 60 GeV) can 
provide the necessary contributions needed to cancel that of the fourth- gen- 
eration whilst at the same time being consistent with LEP bounds. Looking at 
their results, we see that the magnitude of the contribution to the fitted param- 
eters from this sector is highly dependent on M^. Indeed, demanding M-^ ~ 50 
GeV in our model might be too restrictive. One possible solution presents itself 
if we assume that the lightest supersymmetric particle (LSP) is in fact the fourth- 
generation sneutrino N . We have no phenomenological bounds on this scenario, 
the X? would decay invisibly via x? — > v r N assuming the 3 — 4 mixing angle is 
non-zero [S]. This does, however, allow us to increase M^ whilst decreasing Mg. 
The following masses are possible 

JVU ~ 55 GeV 

Mo ^ 105 — 107 GeV 
M e ~ 50 GeV g 

M t = 170GeV 

~ 46 GeV 

The sneutrino N is now the LSP and is stable due to R-parity. The bottom squark 
will decay via B — > tjN (where (i, j) = (c,t) or (b,u T )) depending on which de- 
cay is kinematically allowed. Such decays involve the factor VBtViMj, leading to 
a long lifetime^. If the coupling of N to the Z-boson is small then it is possible to 
have < This allows us to decrease Mg even further, hence increase 
and still retain the decay channels in Eq.(^). The maximum value obtained in this 
case is M£ iax < 60 GeV, consistent with precision fits as shown in [El]. 

The fits to precision data performed in [Q] were for a different set of fourth- 
generation masses. Although the non-universal contributions (see Maltoni, Ph.D. 
thesis [Q] for terminology) remain the same in the limit |JVL — Mj | -C Mz (where 
i, j label the two fermions within the same SU(2)l doublet), the universal contri- 
butions will differ. However, we could arrange to have significant SU(2)l break- 
ing in the (T, B) L doublet if My was found to be sufficiently large. This would pro- 
vide universal contributions to the precision parameters [Q], thus compensating 
for the difference between our model with highly degenerate fourth-generation 
fermion doublets and the fits as performed in Q. 



4 Renormalization Group Study of the MSSM4 

Here we investigate the effect of the fourth-generation on the evolution of cou- 
plings to the GUT scale, where we require gauge coupling unification. This places 
upper limits on the masses of the extra particles to ensure their Yukawa couplings 

6 Current bounds looking at stable/long-lived squarks exclude the range 5 GeV < Mg < 
38 GeV if the mixing in the squark mass matrix is large [ |TTj ] . 
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run perturbatively to the unification scale Mu (h 2 (|j.) < 4n, Mz < |x < Mu). 
Starting at the low-energy scale Mz, the electroweak gauge couplings ai (Mz), 

0C2(M Z ) are fixed through the relations - ,L , = l-^r— nonr^ tt^-itF \ f° r i = 

^ v '-' o Ui(Mz) 5 a, m (M z ] >a em (M z ) 

1,2 respectively. The strong coupling <X3(Mz) is taken from the Particle Data 
Group (PDG) [||| to be 0.1181 ± 0.002. We take the best fit values for a ;™, ( 9 ^ z> 
from the PDG. In principle, one should extract the Z-pole couplings assuming 
the full MSSM4, thereby accounting for the fourth-generation fermions and light 
SUSY spectra in a fully consistent way|^. We have also performed our study in the 
MS scheme although it is the DR scheme that is consistent with supersymmetry. 
However, differences between the DR and MS schemes are not significant at the 
low-energy scale as compared to other uncertainties. 

In our analysis we neglect all Yukawa couplings from the first three gener- 
ations except that of the t-quark whose mass we take to be M? ole = 170 GeV. 
As is typical with four-generation models, we require small values of tan (3 (the 
ratio of the Higgs vev's) so as to avoid h.B(Mz) >0(V4n). Once all couplings at 
Mz have been fixed, we integrate up in energy scale using the two-loop renor- 
malization group equations (RGE). The one-loop leading logarithmic threshold 
corrections from the SUSY sector are accounted for in the numerical procedure. 
We select the point where oci ( \i) — ot2 ( p) as the unification scale Mu with cou- 
pling (Xu(Mu). Any deviation in 0.3 (Mu) = <*u(Mu), which we parameter- 
ize as 5 = a 3( M u)-<*u( M u) can ar j se f rom either of two sectors. On the one 
« u 1 m u ) ' 

hand, we have the MS vs. DR mismatch, experimental errors in a em ( Mz ) / sin 9w 
and the variations in the best fit values of a em (Mz) and sinOw as the fourth- 
generation and light SUSY particles (M < Mz) are included. However, more 
importantly, assuming no intermediate scales, high-energy threshold corrections 
from specific GUT/ string models can provide contributions to 6. Following [Q, 
we note that these corrections (for particular models)can be large. Considering all 
the uncertainties together, we conservatively require unification to within ±5% 
(6 = ±0.05). 

The SUSY threshold corrections in the MSSM4 are important as they can influ- 
ence whether or not a particular set of masses (M t , My, Mb , Mm , Me) will retain 
perturbative consistency to the GUT scale. This can be observed analytically if we 
write the one-loop leading logarithmic correction to the strong coupling from the 
SUSY sector 

1 1 _b^M4_ b |M4 /mv ^ 



aJ(Mz) CX3 (M z ) 2tt V M 

where M co i represents an effective threshold scale and rj ^ iSSM4 / SM4 represents 
the one-loop beta function contribution to the strong coupling in the MSSM4/SM4 
respectively. This correction is implemented at the scale Mz and accounts for the 

7 We account for the one-loop leading logarithmic corrections from the SUSY sector when 
running the RGE by employing the step-function approach This procedure is accu- 
rate in the limit of heavy sparticles but fails for masses M< Mz where both logarithmic 
and finite corrections will influence the extraction of the couplings. 
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coloured sparticles with masses M C oi > Mz m the running of the strong cou- 
pling. oc + / aT represents the renormalized gauge coupling just above/below the 
scale M/. Eq.(^|) holds if the masses of the sparticles are < 2 TeV @. It is obvious 
that the higher the scale M co i, the lower ocj (Mz). Looking at the structure of the 
RGE it is evident that the strong coupling plays an important role in keeping the 
Yukawa couplings perturbative. We obtain an upper bound M™^*, above which 
the initial 0.3 is too small to counteract the effect of the fourth-generation cou- 
plings. 

In practice, we perform a full numerical study, accounting for the threshold cor- 
rections from t, T, B, Mi (where i runs over all sparticles with M > Mz) by chang- 
ing the |3-functions and using the step-function approach in the running of the 
gauge couplings We assume separate degeneracies amongst the coloured 
and weak SUSY spectrum 

M wk = = M t = M ft = M w = M H (4) 

MCOI = Mq = Mg 

We take the fixed value M W |. = 500 GeV though in principle it could be anywhere 
up to ~ 1 TeV. Of course, considering separate degeneracies amongst the coloured 
and weak sparticle spectra is an approximation. In fact, our model demands the 
introduction of light sparticles to ensure consistency with experimental searches 
for the fourth-generation and precision data bounds, thus providing significant 
deviations from degeneracy. Nevertheless, the study does serve to show that the 
combined effect of the coloured spectrum is bounded from above in the process 
of retaining perturbative consistency of the model. See figure [[] for a plot of 
vs. tan 6. In particular we conclude that perturbative unification is only possible 
for the small range 1 .5 < tan (3 < 1 .7. Further restrictions on the allowed tan (3 
range arise from the mass of the lightest Higgs boson that must exceed its exper- 
imental lower bound (see section g). 



5 The Higgs Sector of the MSSM4 

The effect of the fourth-generation on the lightest Higgs mass (Mh) is inves- 
tigated. For our particular solution of (M t , My, Mb, Mm , Me, Mi) that is con- 
sistent with perturbative unification, we see the upper bound Mf^ ax increases 
through loop corrections from the extra particles. Normally, the Higgs sector of 
minimal supersymmetry models is fixed by two parameters (i) tan (3 and (u)Ma 
(the pseudoscalar Higgs mass). We fix tan (3 ~ 1 .63 from the requirement of per- 
turbative unification (see figure @)and so we are left with the one free parameter 
Ma - Assuming this is to be included in the heavy Higgs sector (Mh) we vary 500 
GeV < ||x| = M A < 1 TeV to be consistent with Eq.(@). 

In the MSSM4, for M h > 100 GeV, the channel (h —> invisible), where invis- 
ible represents (NN, EE, NN,XX)/ wm open and would dominate over conven- 
tional h — > bb rates. Exclusion limits will now come from the missing energy 
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Fig.l. Plot of Mf a 1 x vs. tan |3 for M T = 107 GeV; M B = 105 GeV; M N = 46 GeV; M E = 50 
GeV; M t = 170 GeV. 

search (e + e~ — > Zh — > Z + missing energy) that currently sets the lower bound 
at 114.4 GeV [jl6p. We therefore need to check, for our constrained set of pa- 
rameters (My, M b , M N , M E , M t , tan |3), that M™ ax falls above this experimen- 
tal lower bound. We shall employ the one-loop effective potential with contri- 
butions from the top/stop and fourth-generation fermions/sfermions. Since we 
are at low tan (3 we can ignore the contributions from the third generation bot- 
tom/ sbottom masses. 

In general, the sfermion mass matrix can be written as 

/ M§ t + M? + A D i Mt (Ai — n.cot |3] \ 
' VLi ^ Mi (At - ncot P) M|. + M^ + A T5i J 

where Ad and A-jy represents the D-term contributions. We set the soft super- 
symmetry breaking parameters for the coloured squarks Msi = M^ t = M co i 
(i = t, T, B) and, for a given tan 6, we vary M CO i within the range as plotted in 
figure |} For the sleptons Msi = M-g t = |u.| (i = E, N). We randomly vary the mix- 
ing parameters (A t , At, Ab , An , At) and retain the maximum value returned for 
the lightest Higgs mass. As shown in previous studies of the Higgs sector, this 
will occur for A ~ \/~6M co i. Large mixings will induce light sparticles (< Mz) 
which are required in the MSSM4 to provide the decay channels to the fourth- 
generation quarks. Since we are at low tan |3 (~ 1 .63) we require relatively heavy 
squarks in order to increase the Higgs mass above its experimental lower bound. 
We find that M^ QX > 114 GeV for M col > 350 GeV. This is consistent with results 
from section ||. The absolute upper bound (that occurs when M CO i = 500 GeV) is 
M? bs ~ 120 GeV. 





Could there be a Fourth Generation? 



47 



6 Conclusions 

We have seen that it is possible to incorporate a fourth-generation into a GUT 
model, requiring the existence of supersymmetric particles in order to provide 
the necessary cancellations in precision fits. We constrain the fourth-generation 
masses to be M T = 107 Gev; M B = 105 GeV; M N = 46 GeV; M E = 50 GeV. The 
MSSM4 SUSY spectrum needs to be relatively light in order to retain perturbative 
consistency to the unification scale (M co i < 500 GeV). On the other hand (M co i > 
350 GeV) is set from the requirement that the lightest Higgs mass must be above 
its present experimental lower bound. In order to provide decay channels to the 
fourth-generation quarks, it might be that the LSP is the sneutrino N. This has 
implications for dark matter constraints. The SUSY spectrum needed to satisfy 
all constraints can not be obtained from MSUGRA scenarios with universal pa- 
rameters at the unification scale. 
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Another Complex Bateman Equation 
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Department of Mathematical Sciences University of Durham, Durham DH1 3LE 

Abstract. A further class of complex covariant field equations is investigated. These equa- 
tions possess several common features: they may be solved, or partially solved in terms of 
implicit functional relations, they possess an infinite number of inequivalent Lagrangians 
which vanish on the space of solutions of the equations of motion, they are invariant under 
linear transformations of the independent variables, and thus are signature-blind and are 
consequences of first order equations of hydrodynamic type. 

1 Introduction 

This paper is another in a series devoted to an investigation of simple equations 
exhibiting covariance of solutions. These equations have arisen in the study of 
generalisations of the Bateman equation [El], in the equations arising from con- 
tinuations of the String and Brane Lagrangians to the situation where the target 
space has fewer dimensions than the base space [^| and a complex form of these 
equations [|| . The simplest example of these is a complexification of the Bateman 
equation. What we have called the Complex Bateman equation is the following 
equation for a real function cf) defined over the space of variables (xi , *2! Xi , £2); 

4>x, 4>x, 4>x 2 x 2 + <k 2 <l>* 2 4>%i*i - 4>x, fozfolXi - 4>x 2 4>x, 4>x,x 2 =0. (1) 

(Here subscripts denote differentiation). This equation was shown to be com- 
pletely integrable |§][^], with solution given by solving for the following con- 
straint upon two arbitrary functions of three variables, F, G; 

F(cf); xi ,xi) = G(4>; xi ,x 2 ). 

From the form of the solution, or from the equation itself, it is manifest that if $ 
is a solution, any function of 4> will also be a solution and thus that the equation 
exhibits covariance. It is also invariant under separate diffeomorphisms of the 
pairs of variables (xi , X2) and (xi , X2). In fact a subclass of solutions is given by 
the sum of 'holomorphic' and 'antiholomorphic' functions 

4> = f(xi ,x 2 ) + g(xi ,x 2 ). 

A general characteristic of such equations is that they possess an infinite number 
of inequivalent Lagrangians. The equations of motion are partially, or sometimes 
fully solveable in implicit form, as in the examples cited. The fully integrable 
eqautions arise from kinematical first order equations of hydrodynamic type. 

* E-mail:david. fairlie@durham.ac.uk 
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2 Another Complex Bateman Equation 

Now there is another possibility for complexification; we could take instead 

4>x4> x 4>tt - $x4>t4>tx - 4>t4 > x<t>tx + 4>t4H4> X x = 0, (2) 

together with its complex conjugate. These equations also exhibit covariance; cf) 
may be replaced by any function of itself, and the same for ij) and the equations 
remain invariant. Where do these equations come from? Take the hydrodynamic 
equations 

3u 3u 

91 =U ^' (4) 

and set u = v = — i, and the equation (Q), together with its complex con- 
^x 4^ 

jugate are reproduced. Indeed, all that is necessary is to set in an alternative re- 
duction, u = $ an d v = 4> and the same equations arise in consequence. These 
equations admit an infinite number of conserved quantities If S n denotes the 
symmetric polynomial of degree n in u, v, then 

^rS n = ^-(uvS n _i) (5) 
ot ox 

is a conservation law. This is easily proved by induction and from the itera- 
tive definition: S n = u n + vS n -i- These equations can be integrated by the 
usual hodographic method of interchanging dependent and independent vari- 
ables, where they become 

9x 9t 

which can be solved in terms of two arbitrary functions f , g to give 

t = f'(u) + g'(v); x = f(u) -uf'(u) + g(v) -vg'(v) 

with primes denoting differentiation with respect to the argument. If u = cj) and 
v = <j) this parametrisation is a solution to the alternative complexification. Note 
that the requirements that (t, x) be real imposes a further constraint upon the 
functions (f, g). Of course, if (4>, <f>) are treated as independent real functions, no 
such restriction exists. Now the second order equations (||) are Poincare invariant; 
indeed are covariant under general inhomogeneous linear transformations of the 
independent variables. This must be true also for the first order equations @, (0). 
They are clearly translation invariant; if (t, x) transform as 

t' = at + bx, x' — ct + dx, 

then invariance will be maintained if 

, du — c . dv — c 
u = — , v - 



a — bu a — bv 
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3 2-dimensional Born-Infeld equation 

We may remark parenthetically that the same equations(|]),(|J) also yield the gen- 
eral solution to one form of the so-called Born-Infeld equation in two dimensions 
in light-cone co-ordinates [^][0|; 

c)44> tt + fy tt - (A + 2cp x 4> t )4> xt = 0. 

This is achieved by setting 

9 cj> \/A 9cf) VAuv 

dx y/u—y/v' 9t yJ\L—^/v' 

The integrability constraints upon these equations is just the Born-Infeld equa- 
tion itself. Thus the primacy of the first order hydrodynamic equations is again 
manifest. This is a phenomenon which has been noticed before; that the same first 
order equations yield different second order ones depending upon the assump- 
tions made about the dependency of the unknown functions in the first order 
equations upon the functions which enter into the second order equations 



4 Lagrangian 



The construction of a Lagrangian for follows along the lines of 
an auxiliary field \\> and consider the singular Lagrangian 



C 



9$ 9tJj 3t|) 9c£\ If 
"9t 9x ~~ 9t dx J M " 

7 OX 



Introduce 



(8) 



The equation of motion corresponding to variations in the field t]> is simply equa- 
tion (|). Similarly for the variations with respect to 4> we obtain 



i^x^x^tt - •4'x4 ) t4>tx - 4>t4>x4Hx + ihch4>x 



(9) 



i.e. a similar equation with cf) replaced with tj>. But the third equation, correspond- 
ing to variations with respect to <J) is just 



9 

9t 



(4>t^x - 4>xih) 





9 


(*x). 


~ dx 



o. 



(10) 



This is satisfied if ij) is a function of 4~>; then equation ^ is the same as equation 
(§). Incidentally, we see here a situation which has been remarked upon before 
in the context of free field equations [^|, and equations arising from Born-Infeld 
Lagrangians, namely that the Lagrangian itself is a constant, or else a divergence 
on the space of solutions of the equations of motion. It is also evident that the fac- 

4£ /94> a<t>\ 

tor may be replaced by any homogeneous function of I -g— , J of weight 

3x ^ ' 

zero, without affecting the equations of motion. 
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5 Multi-field Lagrangian 

The Lagrangian can be constructed along similar lines to that for the single field; 
one choice is 



c = 



3(X 1 ,X 2> X3) 



9txj , * 2 ) 
3(4)' , <j> 2 
3(xi , x 3 ) 



- CC. 



Variation with respect to gives a combination of the equations of motion for 
cj) 1 and 4> 2 and their complex conjugates; variations with respect to the fields 
4) 1 and 4> 2 yields other linear combinations which together imply the following 
equations,where j takes the values (1, 2); 



det 









*1, 


44 


44 








4> Xl 


4> X2 


*X 3 


*i, 


4>i 


4>x lXl 


4 ) x,x 2 


4 ) X,X 3 


4> X2 


44 




4'x 2 x 2 


4 ) X 2 X 3 


4>i 3 


4> X3 


4>X 3 X, 


4>x 3 x 2 


4>X 3 X 3 



= 0. 



(12) 



and that 9 is a function of 4* 1 and cf> 2 , in much the same manner as xj) is a function 
of 4> in the single field case. As in the case of a single pair of complex fields, these 
equations follow from a set of hydrodynamic equations. 



du 1 
3xi 
dv i 
3xi 



v 1 — + v 2 — =0t=l,2 
3x 3 3x2 



u 



1 Bv 1 

3x 3 



u 



3v x 

3X2 



i = 1,2. 



(13) 
(14) 



Once again, these equations remain the same up to a constant factor under a 
general linear transformation of co-ordinates; this may be seen most easily if they 
are written in a homogeneous notation by introducing vectors ly-, r^; \i — 0,1,2 

such that u 1 = , v 1 = — so that the equations may be written as 



3v 



3u x 
3xm 



making the invariance up to a factor of the hydrodynamic equations under linear 
transformations of the co-ordinates and the vectors n manifest. 

Set u 1 = 4h and choose the following set of these equations and their deriva- 
tives: 



44 +v 1 44+v 2 4, 1 



x 3 







44 +v 1 4> 2 2 +v 2 4) 2 = o 



v2 44 X3 +vi, 



44 X2 +v 1 4>' 2X2 +v 2 4>I 2X3 +vlM 2 + V " 2 4) X3 



4 ) XlX3 



+ v 1 * x , X3 +v 2 ct. x 



x 3X3 



44 = o 

o 
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Eliminate the first derivatives (<$>\ 2 , <p^ 3 ) from the final three equations , solve the 
first pair of equations for (v 1 , v 2 ) and subsitute in the undifferentiated terms of 
the result, setting v 1 = cj) 1 , v 2 — <f> 2 . The ensuing equation is just one member of 

©• 

6 The fundamental hydrodynamic equations 

All second order integrable equations of the type discussed here and in earlier 
work [Ul,^] are consequences of the general first order equations, for which an 
implicit solution may be constructed following Leznov [pUfl. Consider a 2n di- 
mensional Euclidean space with indepenent co-ordinates (x-1, Xi, i = 1 . . .u) and 
construct the differential operators 

3 n ' 9 3 n ^ 3 

° 3x n + ^ U 3xj ' ^ 3x n + ^7 V 3xj ^ 15 ^ 

j=1 j=i ' 

Since Dxt = 0, Dxt = 0, D may be considered a holomorphic differential op- 
erator and D an antiholomorphic operator.. Now imposing the zero curvature 
condition, [D , 0] =0 requires that 

Dv 1 = v\ n +^u i vj c . = 0, Du l = u\ n +^v j u\. = 0. (16) 

These are the general first order equations mentioned above. Since D, com- 
mute, these equations imply that Dv x is a solution to the same equation as v 1 
satisfies. The integration of the equation DDv 1 = requires that Dv l is a general 
anti-holomorphic function, hence, 

Dv l = v\ n +^vhi. = VMv k ;x0, Du 1 = u^+^u'u*. = U^u^xO- (17) 

Indeed f (D)v 1 , for arbitrary differentiable f is also a solution to the equation for 
v 1 . Suppose now we take (n — 1 ) functions 4> v constrained by the (n — 1 ) relations 

QHct^Xk) = P l (* J ;Xic), i=l ...n-1. (18) 

The arbitrary functions Q\ P v depend upon (2n — 1) co-ordinates. They imply 
straightforwardly 

44, = (P^ -Q^r'Q^. <t>k = — (P^j -Ql^r^k- ( 19 ) 

Suppressing the vector indices, suppose u is a function u(4>,x] and v is a function 
u(4>, x). Then the equations ( |l6| ) imply that 

n-1 tl-1 

D4>i = ^ + Y_ v k 4v = 0, Op = 4^ + Y_ uH* k = 0, (20) 
i l 

In other words this requires that cf)' be both holomorphic and antiholomorphic 
in this definition of holomorphicity. Subsituting the derivatives from ( |l^ ) and 
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multiplying on the left by the matrix (P(j, — Q $ ) the equations become 

n-1 n-1 

qL + Y_ v k Qk = o, pi n + y_ u k pi k = o, (2i) 
i i 

which leads to the identifications 

v = -(Qxr'Qx n , u^-fP^P^. (22) 

In consequence any function of x x) will be annihilated by D (D). In partic- 
ular 

D4> = DcJ> = DQ = DP = DQ = DP = 0. (23) 

The last two results follow a forteriori from the equality P = Q. The functions 4>' 
satisfy the multi-field complex Bateman equation &. 

7 Partially integrable covariant equations. 

It appears likely that in the case where the difference between the number of 
dimensions of the base space exceeds that of the target space by more than one, 
the equations of motion are only partially integrable, though this is by no means 
a definitive conclusion. In the case of one field dependent on three co-ordinates, 
the equation which results from the Euclidean Lagrangian 



c = ^ + 4>2 + 4,2 

is 

2 2 2 2 2 2 

4 ) tt(4 ) x + c M+4 ) xx(4 ) u +4 ) t)+4 ) ijij(4 ) t+4 ) x ) = 24) tx 4> t 4 ) x+24v4 ) u4 ) t+24'xu4>x4v 

This equation possesses a large class of solutions given implicitly by solving 
tF(4>) +xG(4>) +yK(4>) = constant 
for <p. It comes from the following first order system; 

uu x + vv x = u t + Vy + v 2 u t — uv(u y + Vt) + u 2 v y 

UV X - VU X = V t - Uy , 

where u = v — This construction suggests further analysis to try to de- 

•Px (Px 

termine whether the system is fully integrable or not. It is also not known whether 
there exist other Lagrangian formulations of these equations. 
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Abstract. Finite groups are of the greatest importance in science. Loops are a simple gen- 
eralization of finite groups: they share all the group axioms except for the requirement 
that the binary operation be associative. The least loops that are not themselves groups 
are those of order five. We offer a brief discussion of these loops and challenge the reader 
(especially Holger) to find useful applications for them in physics. 

1 Introduction 

Many physical systems have symmetries, and groups are the natural mathemat- 
ical objects to describe those symmetries (finite groups for discrete symmetries 
and infinite continuous groups for continuous symmetries). If the elements of a 
group act independently, then the group is abelian; if not, it is non-abelian and 
commutativity amongst the group elements is lost. For discrete groups, this cor- 
responds to an asymmetry of the group multiplication table about its principal 
diagonal, i. e., ab ^ bu for all a and b € G. However, group multiplication is 
associative by definition, 

(ab)c = a(bc) , (1) 

and the concept of nonassociative operations [[[]] has played a limited role in sci- 
ence. Nevertheless, it has not been totally absent. Its main point of entry into 
physics has been through octonions. Also called octaves or Cay ley numbers, they 
define the only division algebra aside from the real, complex and quarternionic 
numbers. An early, but seemingly fruitless, application of non-associativity in 
physics is an octonionic version of quantum mechanics formulated by Jordan, 
von Neumann, and Wigner[^||Q]. Attempts have been made to use octonions in 
particle physics to describe quark structure and other aspects of internal struc- 
ture. For reviews see [^,||^]. There are also an eight-dimensional octonionic in- 
stantons rf|l5|] and applications to superstrings [ p~0|JTT] | . Here we observe that the 
minimal non-associative structures are not octonions, but objects called loops. Let 
us first define them. 



Non- Associative Loops for Holger Bech Nielsen 



57 



2 Loops 

A loop of order n is a set L of n elements with a binary operation [[TJ| such that 
for a and b elements of L, the equations 

ax = b and y a — b (2) 

each has a unique solution in L. Furthermore, a loop possesses an identity element 
e which satisfies: 

ex = xe = x V x e L (3) 

The conditions Eqs.(§) and ^ imply that the multiplication table is a Latin square 
Jl3]jni|l5| ]. The multiplication table of a finite group is such a Latin square, which 
was defined by Euler as a square matrix with n 2 entries of n different elements, 
none occuring twice in the same row or column. 

Any Latin square whose first row and column are identical defines a loop 
whose upper-left entry is the identity element. It follows that any Latin square 
uniquely defines a loop, although different Latin squares may define isomorphic 
loops. This is because a Latin square remains a Latin square under any permu- 
tation of its columns. Thus, one can rearrange any Latin square so that one row 
is identical to one column. Once this is done, that row and column label the ele- 
ments of the loop and their common element is the identity element. 

A system whose multiplication table has non-identical first row and column 
is a quasi-group which is like a loop but which lacks the identity element of Eq. 
(^). We do not consider these structures here. 

In contrast to a group multiplication table, the binary opearation defined by 
a Latin square need not be associative. However, all loops corresponding to Latin 
squares with n < 4 satisfy equation (1). They yield the groups I, Z2 and Z3 at 
orders 1, 2, and 3, and either Z2 x Z2 or Z4 at order 4. 

The situation becomes more interesting at n = 5, for which there are five dis- 
tinct loops. One of these is the group Z5 . The remaining four are non-associative 
loops. For n = 6, there are two groups, Z2 x Z3 and D3, and 107 non-associative 
loops. 

The number of non-associative loops rises very rapidly with n and is known 
only for small values. The number of reduced Latin squares (those in the form 
with identical first row and first column as in all the examples below) is known 
to be 9,408; 16,942,080; 535,281,401,856; 377,597,570,964,258,816 and 

7, 580,721,483, 160, 132,811, 489, 280 

at orders n = 6; 7; 8; 9 and 10 respectively. For n=ll the number of reduced Latin 
squares, and hence the (smaller) number of non-associative loops which corre- 
sponds to the number of isomorphism classes of Latin squares which contain at 
least one reduced Latin square per class, is not yet known (see e.g [|13[Eq]). 

Loops are known to arise in the geometry of projective planes [jlq], in com- 
binatorics, in knot theory and in non-associative algebras, but have yet to 
play a role in physics. Thus we present all the n = 5 cases and some (not all!) of 
the n=6 non-associative loops as a challenge to Holger and others, who may find 



58 P.H. Frampton 



them to be interesting and useful for reasons too subtle to have been revealed to 
us. 

We begin by presenting all of the five n = 5 multiplication tables (see p. 129 
ofQlql) in a form familiar from group theory. Case (la) is the group Z5 (the fifth 
roots of unity) whilst the other four are inequivalent non-associative n=5 loops. 
Case (lb) is special in that the square of any element is the identity element. As 
we discuss is §3, all 5-loops define commutation algebras that satisfy the Jacobi 
identity. 



X 


1 2345 


Xl 


1 23 45 


X2 


1 2345 


1 


1 2 345 


1 


1 2 3 45 


1 


1 2 3 4 5 


2 


2345 1 


2 


2 14 5 3 


2 


2 1534 


3 


3 4 5 1 2 


3 


3 5 124 


3 


3 4 2 5 1 


4 


4 5 12 3 


4 


435 1 2 


4 


4 5 12 3 


5 


5 1234 


5 


5 4 2 3 1 


5 


5 3 4 1 2 


(la) 




(lb) 




(lc) 



x 3 


1 2345 


x 4 


1 2345 


1 


1 2 3 4 5 


1 


1 2 34 5 


2 


2 14 5 3 


2 


2 3 45 1 


3 


3 45 1 2 


3 


3 5 2 1 4 


4 


4 5 2 3 1 


4 


4 15 3 2 


5 


5 3 124 


5 


54 123 


(Id) 






(le) 
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Here we present three illustrative examples of the 107 distinct non-associative 6- 
loops. Each of these defines a commutation algebra that satisfies the Jacobi iden- 
tity: 



xf 


1 2345 6 


A 


123456 


1 


1 2345 6 


1 


1 23456 


2 


2 143 65 


2 


2 16534 


3 


35 1642 


3 


3 6 1245 


4 


465 123 


4 


452 163 


5 


5362 14 


5 


5346 12 


6 


6425 3 1 


6 


64532 1 



xf 


1 2345 6 


1 


12345 6 


2 


2 15643 


3 


34 15 62 


4 


436 125 


5 


5 623 14 


6 


65423 1 
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The following two examples of non-associative n=6 loops define commutator 
algebras that fail to satisfy the Jacobi identity: 



X 6 

x 4 


123456 


xf 


1 2345 6 


1 


123456 


1 


123456 


2 


2 14563 


2 


265 134 


3 


3 6 1245 


3 


3 142 65 


4 


456 132 


4 


43 65 12 


5 


532614 


5 


54 1623 


6 


645 3 2 1 


6 


65234 1 



3 Physics Challenge 

In this section we suggest a few possible applications of loops to physics. We chal- 
lenge the reader to develop a useful application to physics from these notions or 
any others. First, it may be useful to point out that the condition of associativity 
which is required of groups is a natural condition for symmetry transformations, 
since it is an automatic consequence of the composition of mappings. Such map- 
pings between particle states, or between states in a Hilbert space, give rise to the 
familiar symmetry groups. Groups themselves act as transformation groups on 
themselves, and this action is consistent with the group action because of associa- 
tivity. For a finite group, for instance, the multiplication table of the group gives 
a representation of the group as a set of n permutations 

gd9j) =7ti(gj) = 9i x gj (4) 

and clearly 

(9i x 9j)(g) = 9i(9j(g)) (5) 

is a consequence of associativity. For a loop multiplication table, we again get a 
set of permutations, but the multiplication by composition of the permutations is 
not consistent with the loop multiplication, for the same reason. Thus our intu- 
ition about groups as transformations may be a hindrance in interpreting loops 
in physical applications. 

I. One could imagine defining a group product in a way similar to the defi- 
nition of a q-deformed bosonic commutator algebra where a fermionic anticom- 
mutator piece is added, i.e., here we would consider an associative group algebra 
product a • b deformed by a non-associative loop algebra piece a * b to generate 
an algebra with product 

a<g)b = (1 -e)a-b + ea*b. (6) 
This may be a way of introducing dissipation or decoherence into a system. 
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II. We could try to start with a space S and factor out a loop L similar to an 
orbifold construction, where a finite group is factored out. Such an S/L loopifold 
could have application in string theory although its implementation is made non- 
trivial by the absence of matrix representations of the loop. 

III. It is also a consequence of nonassociativity that a representation of a loop 
in terms of linear transformations is never faithful; since matrix multiplication 
associates, the nonassociativity must be annihilated in any map from the loop to 
operators on a vector space. In order to bypass this obstacle, it is useful to con- 
struct an object familiar to finite group representation theory, a loop (or group) 
algebra. We take formal linear combinations of the elements of the loop (with co- 
efficients in R or C), with multiplication carried out termwise according to the 
loop multiplication table. This procedure defines a vector space whose basis el- 
ements are the loop elements and a natural (but non-associative) multiplication 
operation between vectors. We denote the non-associative algebra corresponding 
to a group L as A(L). 

In particular, the loop elements themselves act as linear transformations on 
A via either left- or right-multiplication. If L is a group, this action admits the de- 
composition of A into subspaces corresponding to the irreducible representations 
of the group. For non-associative loops, the situation is less clear because matrix 
multiplication does not follow the loop multiplication. 

However, the algebra associated with a loop has another interesting prop- 
erty. To any A(L) (associative or not), we may define the bracket of two elementa 
a , e A as 

[a, b] = a x b - b x a . (7) 

It is evident that this operation yields an element of A, and furthermore that it is 
antisymmetric: [a, b] = — [b, a]. However, for non-associative loops it is far from 
evident that the bracket operation satisfies the Jacobi identity, 

[a,[b l c]] + [b,[c,a]] + [c,[a,b]]=0. (8) 

Eq. (^) is always satisfied if L is a group. Every finite group, through the commu- 
tator algebra thus defined, corresponds uniquely to a Lie algebra. 

What we find fascinating is that some (but not all) non-associative loops do 
yield bracket operations that satisfy the Jacobi identity, thereby defining commu- 
tator algebras that are Lie algebras. Curiously (and as indicated above), all of the 
non-associative loops with n = 5 are of this class, but only some of those with 
n=6. 

One could imagine using loops as objects to replace flavor or horizontal sym- 
metries in particle physics, or using them as "pregroups." For example, let us 
rewrite the x j loop of Table (lb) in the form 

x i 1 a b c d 

1 1 ab c d 

a a 1 c d b 

b b d 1 a c 

c c b d 1 a 

d d c a b 1 
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For this case, the bracket operation of the loop algebra satisfies the Jacobi identity. 
The structure of the algebra is revealed in terms of the linear combinations 

K= (a + b + c + d)/2 
u, = (a + b-c-d)/2 
u 2 = (a-b + c-d)/2 

and 

u 3 = (a-b-c + d)/2 

The bracket operation reveals that K (and the identity element) commute with 
the other operations and the U| satisfy the su 2 algebra [ut,Uj] = — 2e ijk u k . The 
nonassociativity lurks still in the products of these elements, resembling a twisted 
version of the Pauli matrices; in this basis they are given byKxUi=UjXK = 
— ui/2, ui XU2 — 3u 3 /2, U2 x ui = — U3/2, and cyclic permutations of these. We 
also have the relations K 2 = 1 + ^ and u 2 = 1 — K/2. It is interesting to note that 
the combination 1 — K/2 commutes and associates with the other elements, and 
the relation Iiu? = 3(1 — K/2) suggests an interpretation as a Casimir operator 
for the SU2; we leave this and other details for the interested reader to interpret 
and, hopefully, apply to physics. 
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Particles, Fluids and Vortices 

J.W. van Holten* 
NIKHEF, Amsterdam NL 

Abstract. Classical particle mechanics on curved spaces is related to the flow of ideal flu- 
ids, by a dual interpretation of the Hamilton-Jacobi equation. As in second quantization, 
the procedure relates the description of a system with a finite number of degrees of free- 
dom to one with infinitely many degrees of freedom. In some two-dimensional fluid me- 
chanics models a duality transformation between the velocity potential and the stream 
function can be performed relating sources and sinks in one model to vortices in the other. 
The particle mechanics counterpart of the dual theory is reconstructed. In the quantum 
theory the strength of sources and sinks, as well as vorticity are quantized; for the duality 
between theories to be preserved these quantization conditions must be related. 



1 Particles 

The free motion of a classical particle with unit mass, moving in a smooth space 
with metric gij (x) is described by the Lagrangean 

1- = \ 9ij (x)xSc' , (1) 

where as usual the overdot represents a time-derivative. The Euler-Lagrange equa- 
tions imply that the particle moves on a geodesic: 

^=«i + r Jlc W = o. (2) 

The canonical formulation of this theory is constructed in terms of the momenta 

Pi = ^ = 91^, (3) 

and the hamiltonian 

H = lgVPj- (4) 

The time-development of any scalar function F(x, p) of the phase-space co-ordin- 
ates is then determined by the Poisson brackets 

dt~ 1, *~dxid Vi 3 Pi 3xi' () 
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In particular the Hamilton equations themselves read 

opt ox 1 

A third formulation of the classical theory is provided by Hamilton's principal 
function^ S(x, t), which is the solution of the partial differential equation 

|| = -H(x,p = VS). (7) 
For the case at hand this Hamilton-Jacobi equation takes the form 

— =--g ij V t SVjS. (8) 

at 2 1 

Particular solutions S are provided by the action for classical paths x x (t) obeying 
the Euler-Lagrange equation (§), starting at time T = at an initial point x l {0), 
and reaching the point x x [t) — x v at time X — t: 



S(x,t) 



L 



dxL(x,x) 

o 



(9) 



An example of the class of theories of this type is that of a particle moving on the 
surface of the unit sphere, S 2 . A convenient co-ordinate system is provided by the 
polar angles (0, (p), in terms of which 

L(9,cp] = ]- fe 2 + sin 2 9cp 2 ) , (10) 



2 

for a particle of unit mass. The corresponding hamiltonian is 

with the momenta and velocities related by 

Pe =0, P<p =sin 2 6cp. (12) 

The second equality ( flit ) relates the hamiltonian to the Casimir invariant of angu- 
lar momentum, the components of which are constants of motion given by 

T x = -sincppe - cos cp ctgBp^, ] v = coscppe -sintpctgBp^, J z =y v . 

(13) 

The geodesies on the sphere are the great circles; they can be parametrized by 

cos 6 (t) = sin <x sin cu (t — t* ) , tg (cp(x) — cp*) = cos a tg cu(t — t*), (14) 

where a is a constant, and t* and cp* are the time and longitude at which the orbit 
crosses the equator: 9* = 7t/2. On these orbits the angular frequency is related to 
the total angular momentum by 

tu 2 =2H=J 2 , (15) 



1 The terminology follows ref.[ 
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Observe that, for an orbit reaching the point with co-ordinates (9, cp) at time T*+t, 
the following relations hold: 



cos w = sin 9 cos(cp — cp*J, 

COS0 

sin a = 



sin cot = \h — sin 9 cos 2 ( cp — cp* ) , 



(16) 



1 — sin 2 9cos 2 (cp — cp, 



The last equation implicitly describes the orbit 9(cp), defining a great circle which 
cuts the equator at 9 = 9* = tc/2 and cp = cp*, at an angle a defined by the 
direction of the angular momentum: 



Iz 



cos a, 



Jj 



sin a, 



J2 + J2 

■>X T )y ■ 



The Hamilton-Jacobi equation for this system reads 



as 
at 



+ 



i 



dS 



sin 2 9 V 9( P 



The solution corresponding to the orbit (|14|) is 

1 



S(9, cp,t) = — arccos^ LsmH coslcp — cp„jj , 



(17) 



(18) 



(19) 



which satisfies the equations 
3S 



90 



as 



Pe 



3cp 



cucos9 cos(cp — cp*) 
1 — sin 2 9cos 2 (cp — cp, 

tu sin 9 sin(tp — cp*) 
1 — sin 2 9 cos 2 (cp — cp*) 



(20) 



as cu 2 
aT -~ H --— ■ 

In this approach, the expressions on the right-hand side are obtained by defining 
w via the last expression, in agreement with (fl6|). The same principle of energy 
conservation/ time-translation invariance implies that S does not depend on t*. 

2 Fluids 



The Hamilton-Jacobi equation (|S|) can itself be obtained in a straightforward way 
from a variational principle: introduce a Lagrange multiplier field p(x) and con- 
struct the action functional 



A(p,S) 



dt 



(21) 
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The square root of the (time-independent) background metric has been included 
to make the integration measure invariant under reparametrizations. Of course, 
we could absorb it in the definition of Lagrange multiplier field, but then p would 
transform as a density rather than as scalar. 

The Hamilton-Jacobi equation follows by requiring the action to be station- 
ary wr.t. variations of p: 

On the other hand, the stationarity of A(p, S) w.r.t. S gives 

This equation can be interpreted as the covariant equation of continuity for a 
fluidQ with density p and local velocity 

Vi = ViS ^+V 1 (pv i )=0. (24) 

In this interpretation the gradient of the Hamilton-Jacobi equation gives the co- 
variant Euler equation 

^+ v JV J v t = 0, Vjvt = ^ - V n W- (25) 



Eq.(24) states that the fluid flow is of the potential type. Indeed, in the absence 
of torsion the Riemann-Christoffel connection Py k is symmetric and the local 
vorticity vanishes: 

ViVj - VjVt = 0. (26) 
For the fluid flow to be incompressible, the velocity field must be divergence free: 

V • v = AS = 0, (27) 

where A = g 1 ' Vi Vj is the covariant laplacean on scalar functions over the space. 
It follows that the number of incompressible modes of flow on the manifold 
equals the number of zero-modes of the scalar laplacean. For example, for flow 
on the sphere S 2 (or any other compact Riemann surface) there is only one incom- 
pressible mode, the trivial one with v 1 = everywhere. 

For a given geometry gtj (x), the solution of the Hamilton-Jacobi equation 
(53) provides a special solution of the Euler equation (p5|); for a conservative sys- 
tem: 3S/3t = — H = constant, it implies 3vt/3t = and v> VjVi = 0. Accordingly, 
this solution describes geodesic flow starting from the point (9*, cp*). 

To turn this into a complete solution of the fluid-dynamical equations (^i|), 
( j25| ) it remains to solve for the density p. The equation of continuity takes the form 

^ + V i (pV i S)=0. (28) 
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It follows that a stationary flow, with p not explicitly depending on time t, is 
possible if 

V-(pVS)=0. (29) 

In addition to the trivial solution p = po = constant, v = VS/m = 0, it is possible 
to find non-trivial solutions of equation ( |2^ ) for spatially varying density p. As 
an example, we consider flow in a 2-dimensional space; in this case one can in- 
troduce a generalized stream function T(x, t), dual to the fluid momentum, and 
write 

pV^-S = — e^VjT. (30) 
Then for theories of the type (g): 

e l *VtSVjT e ij ViSVjT 



^g(VS) 2 2ryg 



(31) 



With H constant, the factor 2H in the denominator can be absorbed into the defi- 
nition of T = T/2H, and hence the density is given by 

P = -J= e^'ViSVjT = -L e^ViVjT, (32) 

for the pseudo-scalar function T the gradient of which is dual to pVS. Note also, 
that eq.@ implies VS • VT = v • VT = 0. 

As an illustration, we again consider the unit sphere S 2 . The velocity field is 
given by the momenta (^) per unit mass: 

w cos cos(cp — cp*) tusin9 sin(cp — cp*) 

ve = ; =, v v = —== (33) 

1 — sin 2 9 cos 2 (cp — cp*) \/1 — sin 2 9 cos 2 



cos^icp — cp*j y I — sin t) cos ( cp — cp* J 

Taking into account that on the sphere the non-vanishing components of the con- 
nection are 

V = r^-sinGcosG, (34) 

a straightforward calculation shows that indeed 

Vj v j =cu 2 , WjVi^O, ^=°- (35) 

The first two equations actually imply v' (VtVj — Vj v\) — 0, in agreement with the 
absence of local circulation. From these results it follows, that the flowlines are 
geodesies (great circles) given by eq.(|l6|), and stationary. 

For the gradient of the stream function T to be orthogonal to the velocity field 
(33 1, it must satisfy the linear differential equation 

v-VT = <=> tg(cp - cp„) VcpT = sin9 cos9V e T. (36) 

The general solution can be obtained by separation of variables, and is a function 
of the single variable: T( 9, cp) =f (y), withy = tg0 sin(cp — cp») = ctga.Forsuch 
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a scalar field 



(37) 



The corresponding density p is then 



P(0, <P) = ^ = — s nv)\ v=a , K ■ (38) 

cos 9 cusinacosB y 8 



!Q p*srna p* 

P 0, <P = h~ = / 39 ) 

cos0 !■* a mil 



The simplest, most regular solution is obtained for p(<x) = p* sin oc: 

p* 

^1 — sin 2 9 cos 2 ( cp — cp 

This solution corresponds to 

T(9, cp) = cup, <x(9, cp) f(y) = cup, arcctgu. (40) 

Observe, that in this case T, like oc, is an angular variable; indeed, a increases by 
27m on any loop winding around the point (9 = 7r/2;cp = cp*)u times. 

The solution ( p9| ) possesses singular points at 9 = n/2, cp = cp* + rut, corre- 
sponding to a source for n = 0, and a sink for n = 1 . This can be established from 
the expression for V • v: 

a)sin9cos(cp - cp*) 
V ■ v = —= (41) 

J 1 — sin 2 9 cos 2 (cp — cp* ) 

which becomes (+oo, — oo) at the singular points. However, a more elegant way 
to establish the result, is to make use of the stream function ( fiO| ) and consider the 
flux integral 

<D(D =<f dlpv n> (42) 

. r 

representing the total flow of material across the closed curve V per unit of time. 
Consider a contour Y winding once around the singularity at (9 = 7t/2; cp = cp*); 
on such a curve a increases from to 2n. Then 

®(D = o VgEijP vi dx' =o ViTdx 1 =2ttcup*. (43) 
. r Jr 

This represents the total flow of matter from the hemisphere centered on the 
source at (9 = 7t/2; cp = cp„) to the hemisphere centered on its antipodal point, 
the sink at (9 = 7t/2; cp = cp* + n). 

3 Vortices 

The dual relationship between the velocity potential S and the stream function T 
suggests to study the dynamics of a fluid for which T is the velocity potential: 

Vi = — ViT. (44) 

p* 
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The constant p* has been included for dimensional reasons. Like before, this ve- 
locity field is stationary: 3vi/3t = 0, but it is not geodesic. Indeed, the velocity 
field describes motion under the influence of an external potential; specifically: 

v-Vvi-lvtv 2 =-Lvi(VT) 2 = _L Vi (pVS) 2 . (45) 
2 2pj 2pi 

Here p(x) and S(x) denote the previously defined functions mapping the man- 
ifold to the real numbers — e.g. ( p^ ) and ( ^ ) for fluid motion on a sphere — ir- 
respective of their physical interpretation. Now again, as (VS) 2 = 2H = w 2 = 
constant, it follows that 



w 



: 



V ■ Vv i = TZ? V i P 2 = - V i h - ( 4 6) 

Combining eqs.(|||) and @): 



2P 2 



lv 2 =-(h-h ) = ^, (47) 

where h represents the external potential. Because of the potential nature of the 
flow, eq. (p3|), the local vorticity again vanishes: ViVj — VjVt = 0, but as eq.(^3[) 
shows, this is not necessarily true globally. Indeed, in singular points of the orig- 
inal geodesic fluid flow (with sources / sinks), the dual flow generally has vor- 
tices/ anti-vortices. 

Continuing our example from the previous sections, we can illustrate these 
results in terms of flow on the unit sphere, for which T/p* = wot and = cuVtix: 

w sinfcp — cp*) cusin9cos8cos(cp — cp*) 

ve = 2 > v <p = 5 ' ' ' 

1 — sin 9 cos 2 (cp — cp*) 1 — sin 9 cos 2 (cp — cp*) 

It follows, as expected, that 

v 2 = cu 2 (Va) 2 = ^ = . (49) 

P* 1 — sin 9 cos 2 (cp — cp*) 

A further remarkable property, is that the dual flow is divergence free: 

V-v = & Aa = 0, (50) 

where-ever v is well-defined; obviously, the result can only be true because of the 
two singular points (9 = n/2; cp = cp*) and (9 = n/2; cp = cp* +n), where Vi and 
its divergence are not well-defined, i.e. topologically the velocity field is defined 
on a cylinder, rather than a sphere. These two points are centers of vorticity, as 
follows directly from eq.(^), which in the present context can be rewritten as 

o Vidx 1 = 2nw, (51) 
. r 

for any closed curve V winding once around the singular point (n/2, cp*); as this 
curve also winds once around the other singular point in the opposite direction, 
they clearly define a pair of vortices of equal but opposite magnitude. 
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As the flow is divergence free, it follows that in this case there can be non- 
trivial incompressible and stationary flow modes: for constant density pi one has 

^ = 0, V Pl =0 V-(piv)=0, (52) 

and the equation of continuity is satisfied. 

The nature of the flow lines defined by eq.(^8|) is clear: they are parallel circles 
of equidistant points around the centers of vorticity. On these circles the velocity 
is constant in magnitude, implying by ( fl^ ) that sin 9 cos( cp — cp*) = cos |3 = con- 
stant. For example, for cp* = we get x — cos |3 = constant; the flow line then is 
the circle where this plane of constant x cuts the unit sphere. On these flow lines 

ve = — tui sin(cp — cp*), v v — — tv-\ cos (3 cos 6, (53) 

with 

v 2 u> w 

CUl = = J ; = 2 — • v**' 

w 1— sin 9 cos 2 (cp — cp*) sin |3 

4 The dual particle model 

Having clarified the nature of the (incompressible) flow described by the dual 
velocity potential T/p*, we now reconstruct the corresponding particle-mechanics 
model for which T/p* is Hamilton's principal function. From eqs.(^), ( fl^ ) we 
observe that the hamiltonian is of the form Hi = K + h, with for the specific case 
at hand a kinetic-energy term: 

K = 1 gVpj -> \ fpe + - 05) 

2 2 \ sin 9 / 

and the potential (normalized for later convenience such that 2H = cucui ): 

h(9,cp) =h - -r-2 >— ^— 1-- rx; 77 r • (56) 

2p 2 2 V 1 -sin 2 9cos 2 (cp - cp*)/ 

The corresponding lagrangean Li = K—h produces the Euler-Lagrange equations 

^ . . . . 2 cu 2 sin 9 cos 9 cos 2 (cp — cp*) 
Pe = 9 = sin 9 cos 9 cp 1 



{1 — sin 2 9 cos 2 (cp — cp*)) 2 



(57) 

d /. i \ cu 2 sin 2 9 sin(cp — cp*) cos(cp — cp*) 

P cp = — sin 9 cp ) — = . 

v dtV V J (1 -sin 2 9cos 2 (cp-cp*)) 2 

These equations have solutions 

cos 9 — sin (3 sin cu 1 1 , tg (cp — cp*) = tg |3 cos cuit, (58) 
with |3 a constant, implying the relation 

sin 9 cos(cp — cp*) = cos (3. (59) 
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Solving for the velocity (and taking into account the unit mass) 

p e = v e = — cui sin(tp — cp»), p v = sin 2 v v — — a>i cos |3 cos 6, (60) 



in agreement with ([53|). From these results we can compute Hamilton's principal 
function 



Si(9,«p,t) = 



drLi [9(t), <p(t)] = ^arcctg 2 (tgesin(cp-(pj). (61) 



This function indeed satisfies the Hamilton-Jacobi equations 

3$i _ 9Si 
^ a — Pe> "r — Pep i 

with (pe , Pep ) as given by eq.(|60|), and 



3Si 

at 



CUO>1 



3Si 



3Si 



sin 2 6 V 9( P 



-h(6,(p). 



(62) 



(63) 



Using the relation ctg a = tg sin( cp — cp„) = ctg a> 1 1, the equations (|62|) can be 
recast in the form 

Vi = cuVitx = — VtT. (64) 
p* 

Hence T/p* can indeed be identified with Hamilton's principal function of this 
system. 

Repeating the arguments of sect. ^, the action ( pl| ) for the Hamilton-Jacobi 
theory is now generalized to: 



A(p,S i; H) = 



dt 



d"x N ?(>( ^ + Ig^ViSiVjSi +h 



(65) 



Reinterpretation of Si as a velocity potential for fluid flow: v = VSi, leads back 
directly to the inhomogeneous Euler equation 



9vi 

— + v • Vv t = -Vth, 
ox 



(66) 



which for stationary flow becomes eq.(|46|). Variation of this action w.r.t. Si gives 
the equation of continuity for p, as before; note that in this action h plays the role 
of an external source for the density p. 



5 Quantum theory 

The quantum theory of a particle on a curved manifold is well-established. For 
the wave function to be well-defined and single-valued, the momenta must sat- 
isfy the Bohr-Sommerfeld quantization conditions 



o pidx 1 = 2nnh, 
> r 



(67) 
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for any closed classical orbit T. For the free particle of unit mass on the unit sphere 
the left-hand side is 

■T 

v 2 dT = cu 2 T = 27tu), (68) 



where T = 27t/cu is the period of the orbit. Hence the quantization rule amounts 
to quantization of the rotation frequency (the angular momentum): cu = nh. 
For the dual model, the same quantity takes the value 



cp vtdx 1 — 
r 



v 2 dT = — = — =cua>iTi —2mx>, (69) 
o sin |3 



and again cu = nh. As the quantization conditions in the two dual models are the 
same, the duality can be preserved in the quantum theory. 

If this is to be true also in the fluid interpretation, the quantization conditions 
must be respected at that level as well. Now the first quantization condition for 
the integral ( |68| ) is interpreted in the fluid dynamical context as a quantization 
of the fluid momentum, cf. eq.(^). The second quantization condition ( ^9| ) has a 
twofold interpretation: first, according to eqs.(|42|), ( p3| ) it quantizes the strength of 
the fluid sources and sinks in the model of free geodesic flow; the agreement be- 
tween the two quantization conditions is then obvious: in order for the strength 
of the source / sink to satisfy a quantization condition, the amount of fluid trans- 
fered from one to the other must be quantized as well. 

In the context of the dual model however, the condition imposes the quanti- 
zation of vorticity in the quantum fluid [||. In the more general context of quan- 
tum models of fluids in geodesic flow on a compact two-dimensional surface and 
their duals described by the stream functions, this observation shows that du- 
ality at the quantum level requires the quantization of sources in one model to 
be directly related to the quantization of vorticity in the dual one. This situation 
closely parallels the relation between the quantization of monopole charge [Q| 
and the quantization of the magnetic flux of fluxlines [0] in three dimensions. 
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Results on 2D Current Algebras 

J.L. Petersen* 

The Niels Bohr Institute, Copenhagen, Denmark 

1 Introduction 

A brief account i presented on a series of results on free field realizations of 2-d 
affine current algebras and techniques for working out correlators, in particular 
in non-integrable, admissible representations. A more detailed account of these 
results, obtained in collaboration with Yu Ming and Jorgen Rasmussen, may be 
found i refs. §§§||§|[§].The use of these may be used in inderect ways in 
formulations of perturbative string vacua and 2-D quantum gravity. 



2 SL(2) current algebra and admissible representations 

We first describe how to obtain conformal blocks for SL(2) WZW theories in the 
case of non-integrable representations, in particular for admissible representa- 
tions [g|H]]. For previous approaches to the same problem, see 
Our goal here is to obtain a formulation based on the Wakimoto [ 24 1 free field 
realization. The problem is complicated by the need for introducing a second 
screening charge in the case of admissible representations. This screening opera- 
tor involves a fractional power of a free (antighost) field [p5|]. Also the formalism 
requires the introduction of several other fractional powers of free fields. We have 
found that everything may be treated rather neatly by means of fractional calcu- 
lus. 

As an explicit verification that our formalism works we have managed to 
prove that the conformal blocks we obtain satisfy the Knizhnik-Zamolodchikov 
equations p6[p7| ]. As an additional bonus we have been able to provide a proof 
of an interesting suggestion by Furlan, Ganchev, Paunov and Petkova [ |l3| ] for 
how conformal blocks of the SL(2) WZW theory reduce to the conformal blocks 
of minimal models. 

When the level of the affine SL(2) k algebra is k, we define 

k + 2 = t = p/q (1) 

where p, q are positive coprime integers for admissible representations. For these 
there are degenerate representations whenever the spin is given by 

2j r ,s + 1 =T-st, r = 1,...,p-1, s=0,...,q-1 (2) 
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The Wakimoto free field realization is in terms of a scalar field, cp and a pair of 
bosonic dimension (1,0) ghosts, ((3,y): 

(p(z)cp(w) ~ log(z- w), |3(z)y(w) ' 



z — w 
J + (z) =(3(z) 

J 3 (z) =-yP(z)- Wlacp(z) 



J-(z) =- T 2 [3(z)+k9 Y (z)-\/2tY3cp(z) (3) 



They satisfy 



J+(z)J-(w) J 3 (w) + 

z — w 

J 3 (z)J±( W ) ~±— — T^w) 

z — w 

T3f.u3 f ..^ k /2 



r(z)r( W ) 



(z- w) ; 



(4) 



The central charge of the corresponding virasoro algebra is c = 3k/t. Fateev and 
Zamolodchikov [ ]28| ] introduced a very useful formalism for primary fields, which 
we shall adopt. In general there is a multiplet of primary fields (^(z). One com- 
bines these introducing an extra variable, x as follows 

4>,(z,x) = 22<)>nz)x>'- m (5) 

m 

For integrable representations 2j is an integer and we simply get a polynomial in 
x. However, for fractional spins we have highest weight, lowest weight or contin- 
uous representations, and the x-dependence can be arbitrarily complicated. The 
new primary field satisfies the following OPE 

J Q (z)cMw,x) - — — D>j(w,x) 
z — w 

D+ = -x 2 9 x +2xj, D x = -x9 x + j, D X =9 X (6) 
One easily verifies that 

4>j(z,x) = [1 +xy(z)] 2i e-'\ / * cp(z) (7) 
Finally there are the two screening charge currents [^] 

Si(z) = p(z)eVT<PW, S_ t ( Z ) = P(z)- t e- t ^W = (S 1 ( Z ))- t (8) 

We see in these last equations the need for being able to treat fractional powers of 
free fields. 
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Our treatment of Wick contractions in such situations is based on the follow- 
ing identity, trivially valid for — t a positive integer, but non-trivial for general 
t: 

1 _ 



P(z)- t F(Y(w))=:[P(z) 

= L 

neZ 



- t F(y(w)) 



-t 
n 



z — w 

:fi n (z)(z-w) t+n d-t-?Hy(w)): (9) 



Examples of the use of fractional calculus [gfj are me Riemann-Liouville operator 

1 



3- Q f(z) = 



r(a) J 



z-t 



a-1 



f(t)dt, a > 



and 



r(b + i; 



,b-a 



T(b-a+l) 

In addition we shall need unconventional (asymptotic) expansions like 



e x = 3X = V — 

z— r(n — a ■ 
neZ 



r 



(10) 

(ii) 

(12) 



for x an operator. 

According to the above, we may treat the free field representation of an N- 
point conformal block in terms of the following integral of a free field correlator: 



N-1 



n=2 



oll^^We-^Wxon^PlwOe^^'lj,) (13) 

J k=l J 1=1 



Here r and s are the number of screening charges of the first and second kind 
respectively. There has to be a precise choice of the bra and ket as described in 



By carefully analyzing the three-point function it is possible to verify fusion 
rules previously obtained using different techniques. For a three point function, 
jl + jz ~ H — r ~ s "t with r, s being the number of screening operators. Using 
2ji + 1 = i*i — Sit, we find the following fusion rule, already written down in 



1 + |ti - r 2 | < r 3 < p — 1 — |ri + r 2 - p| 
si — S2 1 < S3 < q — 1 — |si + S2 — q + 1 

referred to as their rule I. There is in addition a rule II: 

1 + |p - ti - t 2 | < t 3 < p - 1 - h - r 2 | 
1 + |q — si - s 2 — 1| < s 3 < q — 2 - |si — s 2 | 



(14) 



(15) 
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It is less trivial to obtain, but it follows both by a proper analysis of the three point 
function and from our analysis of the 4-point functions (see below). 

It is now possible to write down the general N -point function with M = r + s 
screening charges as follows 



W N 



M A 

i=1 



M N-1 

m<n i=l m=1 i<j<M 



P N-1 M N-1 , 

II ^r(2j m + l)u 2 '- , e 1 /^n ( ^ J^j-ic (16) 



where zi = xi = 0. Here the powers ki are —1 and t respectively for screening 
charges of the first and second kind. 

The N-point conformal block, W(zn , xn , zi , X] ), is a function of N pairs 
of variables, (zi,Xt). An interesting proposal of Furlan, Ganchev, Paunov and 
Petkova | |l3| ] is that when we put xt = Z\ this block agrees up to normalisation 
with a corresponding block in the minimal conformal theory with the same p , q 
as for the SL(2) k theory with t = k + 2 = p/q. We have proved this result and 
clarified how it is related to hamiltonian reduction based on T + (z) ~ 1 . 

We have used the above formalism to carry the computation of 4-point con- 
formal blocks generalizing the techniques developed by Dotsenko and Fatteev 
|2^]. The calculation is very technical and we refer to ref.[||] for details. Here we 
merely give the highly non trvial result for the operator algebra. Expressing the 
OPE of two primary fields as 

4) i2 (z,z;x,x)4) j J0,0;0,0)^X (z ^H( i ,) + H( j2 )-Hm C A l A 2 4> i (0,0;0,0) (17) 

We use the following notation (Fusin rule I) 

A = 2jt + 1 

m =j(j + i)/t 

j = h + h -r + st 
G(x) = r(x)/r(1-x) = 1/G(1-x) (18) 

Then we have found 

T t^ G(1 — s + (1 — A] + i)/t)G(1 - s + (1 - A 2 + i)/t) 
X ll G(l+s-(l+A + i)/t) (W) 
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3 Generalization to algebras other than SL(2) 

Much of the above formalism may be generalized to arbitrary algebras. Whereas 
a large amount of work has already been carried out on free field realizations, 
"^|||||||||5@|6p|||||l, we have provided several proofs and new 



results, not least concerning extremely compact general expressions and other 
missing elements, see ref.[|^] for more details. 

Consider some simple Lie algebra of dimension d and rank r. Let{j a } denote 
a set of Lie algebra generators with {e a , f a } being raising and lowering operators 
and {hi} being Cartan algebra generators. The set of (positive) roots is denoted 
( A + ) A. Simple roots are {at}, i = 1 , r, and 9 denotes the highest root. The root 
dual to oc is oc v = 2cc/cc 2 . Commutation relations are 

[ja,jb] =fab C jc (20) 

or in the Cartan- Weyl basis 

[fu, cj = (ocY, cc)e a , [hi, fj = a )f a (21) 

for which 

[e«,f«J = h Bt = G t '(a}',a v )h J , Dn,e,] = At 5 ej (22) 

with ej = e ttj for <Xj a simple root, and where Aij = aY ■ <Xj = Gy0Cj/2 is the 
Cartan matrix. The Car tan-Killing form is tr(j a jb) = K a b- Dynkin labels A^ of 
weight A are defined by 

A = A k A ,k) , A k = (a^,A) 

with {A ,k) } the fundamental weights. We introduce a coordinate x a (generalizing 
the x for SL(2)) for every positive root. Further, define Lie algebra elements 

e(x) =x a e a , f(x) =x a f« (23) 

and group elements 

g+(x) =e e,x » . g_(x) = e f < x » 
and thus on a representation space corresponding to highest weight A 

|A,x) = g_(x)|A) 

It proves extremely convenient to introduce the adjoint matrix representation of 

e(x) 

C h a (x) = C(x) a b = (xPC p ) a b = -xPf pa b (24) 

In a Cartan- Weyl basis, a = ( a, i, — (3 ) , a, (3 S A + and C b = unless either a = — a 
or b = +a with a G A + . We then obtain differential operator realizations from 

j a |A,x) = L_ a (x,3,A)|A,x) (25) 

Let 

oo _ 

= - H ^ (26) 

n=0 
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be the generating function of Bernulli numbers. Using the power series expan- 
sion, we may replace u by any matrix. Define the d x (d — r)/2 matrix 

V(x) a p = (e- c ' x VB(-CM)/ (27) 
with a, (3 , y all denoting positive roots. It is then possible to express the result as 

J a (x,d,A)=V(x) a ii d ti + (e- c M) i a A i (28) 

Notice that (e -c )* = (a e A+) and (e~ c )! = 6j\ For any given algrebra, such 
expressions truncate at finite order. 

The free field realization for the current algebra is in terms of bosonic ghost 
fields, ((3 a (z), y |3 (z)) and Liouville type scalars 4>i(z),i = 1,...,r satisfying the 
OPEs 

= — , <(H(z)(|> J (w) = Gtjln(2-w) (29) 

z — w 

with Gij the metric of eq.(p2|). For central charge, k and k v = 2k/0 2 and dual 
Coxeter number h v , we introduce 

t=^(k v + h v ) 



We seek a set of currents satisfying 

T Q (z)T b (w) = — ^ ! '"I (30) 



Kabk , fab C Jc(w) 



The free field realization may then be obtained from the differential operator re- 
alization simply be the replacements 

3 a ^|3 a (z) , x a ^ T a (z) , At -» Vtdfciz) 

followed by normal ordering. In addition, for negative roots we must add to 
T_ a (z) the "anomalous" term 

F™-(yU),3 Y (z)) = 3 Y a (z)F ap ( Y (z)) 
2k 

F«p(Y) = ^(V- 1 (7))^+(V- 1 (Y))^9<rV^Y)a6V- a (Y) (31) 

where V~ 1 is to be understood as obtained from inverting the sqare matrix with 
positive roots only in rows and columns. All roots indicated, a, P, 6, \i, are posi- 
tive. 

For every simple root, there are two kinds of screening operators, the ones 
of the "second kind" with fractional powers of free fields. We find that they are 
given by (see also [ |l^2|j3^|36yi2| , ^8| | for many results in special cases and differ- 
ent mostly more cumbersome forms) 

sj 1) (z)=:E a .(-Y( Z ),-(3(z))) :: *> (z)/2 ^ : 

s{ 2) (z) = (sj 1) (z))- 2t /«J (32) 
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We have provided proofs previously missing, in particular in the last case, that 
these satisfy the required conditions. 

Finally we have provided explicit expressions for free field realizations of 
primary fields. For a general representation with highest weight A = Ai<A' k ' the 
primary field may be simply expressed in terms of those for the fundamental 
repesentations by 

r 

*A(yH,x) = [][*A"'W w )^)] Al (33) 

k=l 

For the fundamentals we here restrict to the case of SL(N), where roots are con- 
veniently labelled by double indices (ij ) , i, ) = 1 N. Introduce the matrix rep- 
resentation F (N) (x) of f (x) eq.(|23|) by (F^'),. 1 = x 1 ' for i < j and otherwise, and 
let 

e F ' N, <*>(I(k)) 

denote the k x k matrix obtained from the N x N matrix e F ' N) ,x ' by using the 
first k rows, and the columns given by the set I(k) = {ii , iy}. Then we have 
obtained 

cp A(k . (y(z),x) = Y_ det (e F ' N) (x) (Kk))) det (e F ' N) (I(k))) (34) 

Kk} 

Using the screening charges, we may then obtain correlators for non-integer Dynkin 
labels of the form 

A lc = A ki T i -G lci sH = tk-t k ^-t (35) 

with t = -iji. — k v + h v . These are degenerate representations for r 1 , s x ,tk,3k 
integer, [Bp, and admissible representations for $ rational, [p9||. 
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1 Introduction 

The modern physics of electroweak and strong interactions is described by the 
Standard Model (SM), unifying the Glashow-Salam-Weinberg electroweak the- 
ory and QCD - theory of strong interactions. 
The gauge group of the SM is: 

SMG=SU(3) c xSU(2)lxU(1)y, (1) 

which describes the present elementary particle physics up to the scale ~ 1 00 GeV. 

Considering the physical processes at very small (Planck scale) distances, 
physicists can make attempts to explain the well-known laws of physics as a 
consequence of the more fundamental laws of Nature. Random Dynamics (RD) 
was suggested and developed in Refs.[[(|]-[^]] as a theory of physical processes 
proceeding at small distances of order of the Planck length Ap : 

A P = Mp, 1 , where M Pl = 1 .22 • 1 1 9 GeV. (2) 

Having an interest in fundamental laws of physics leading to the description of 
the low-energy SM phenomena, observed by the contemporary experiment, we 
can consider two possibilities: 

1. At very small (Planck scale) distances our space-time is continuous, and there 
exists a theory with a very high symmetry. 

2. At very small distances our space-time is discrete, and this discreteness in- 
fluences on the Planck scale physics. 

The item 2 is a base of the RD theory. 

The theory of Scale Relativity (SR) [^] is also related with item 2 and has a 
lot in common with RD. In the SR the resolution of experimental measurements 
plays in quantum physics a completely new role with respect to the classical one 
and there exists a minimal scale of the space- time resolution: e m i n = Ap, which 
can be considered as a fundamental scale of our Nature. In this case, our (3+1)- 
dimensional space is discrete on the fundamental level. This is an initial point of 
view of the present theory, but not an approximation. 

The lattice model of gauge theories is the most convenient formalism for the 
realization of the RD ideas. In the simplest case we can imagine our space-time as 
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a regular hypercubic (3+l)-lattice with the parameter a equal to the fundamental 
scale: 

a = A P = 1 /M PV ~~ 1 Cr 33 cm. (3) 

But, in general, we do not know (at least on the level of our today knowledge) 
what lattice-like structure plays role in the description of the physical processes 
at very small distances. 

2 G-theory, or Anti-Grand Unification Theory 
(AGUT) 

Most efforts to explain the Standard Model (SM) describing well all experimental 
results known today are devoted to Grand Unification Theories (GUTs). The su- 
persymmetric extension of the SM consists of taking the SM and adding the cor- 
responding supersymmetric partners |j^]. The Minimal Supersymmetric Standard 
Model (MSSM) shows [^] the possibility of the existence of the grand unification 
point at (iGUT ~ 10 16 GeV. Unfortunately, at present time experiment does not 
indicate any manifestation of the supersymmetry. In this connection, the Anti- 
Grand Unification Theory (AGUT) was developed in Refs.[|@]-[§]] and [|§-@] 
as a realistic alternative to SUSY GUTs. According to this theory, supersymmetry 
does not come into the existence up to the Planck energy scale (g). The Standard 
Model (SM) is based on the group SMG described by Eq. ([!]). AGUT suggests that 
at the scale |j*g ~ M-Pl = Mpi there exists the more fundamental group G contain- 
ing N gen copies of the Standard Model Group SMG: 

G = SMG] x SMG 2 x ... x SMG Ngen = (SMG) Ngen , (4) 

where N gen designates the number of quark and lepton generations. 

If N g e n —3 (as AGUT predicts), then the fundamental gauge group G is: 

G = (SMG) 3 = SMG lstgen . X SMG 2 ndgen. x SMG 3 rdgen. i (5) 

or the generalized one: 

G f = (SMG) 3 x U(1) f , (6) 

which was suggested by the fitting of fermion masses of the SM (see Refs.JT^]). 
Recently a new generalization of AGUT was suggested in Refs.[p7||: 

G ext = (SMG xU(I)b-l) 3 , (7) 

which takes into account the see-saw mechanism with right-handed neutrinos, 
also gives the reasonable fitting of the SM fermion masses and describes all neu- 
trino experiments known today. 

By reasons considered in this paper, we prefer not to use the terminology 
"Anti-grand unification theory, i.e. AGUT", but call the theory with the group 
of symmetry G, or Gf, or G ext , given by Eqs.Q-(^), as "G-theory", because, as 
it will be shown below, we have a possibility of the Grand Unification near the 
Planck scale using just this theory. 
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The group Gf contains the following gauge fields: 3 x 8 = 24 gluons, 3x3 = 9 
W-bosons and 3x1+1 =4 Abelian gauge bosons. The group G ext contains: 
3 x 8 = 24 gluons, 3x3 = 9 W-bosons and 3x1+3x1 =6 Abelian gauge bosons. 

There are five Higgs fields in AGUT, extended by Froggatt and Nielsen [[l5|] 
with the group of symmetry Gf given by Eq.(^). These fields break AGUT to the 
SM what means that their vacuum expectation values (VEV) are active. The au- 
thors of Refs.[[l5| used three parameters - three independent VEVs with aim to 
find the best fit to conventional experimental data for all fermion masses and 
mixing angles in the SM. The result was encouraging. 

The extended AGUT by Nielsen and Takanishi fll^], having the group of sym- 
metry G e) ct (see Eq.(^])), was suggested with aim to explain the neutrino oscilla- 
tions. Introducing the right-handed neutrino in the model, the authors replaced 
the assumption 1 and considered U(48) group instead of U(45), so that G ext is a 
subgroup of U(48): G ext - U(48). This group ends up having 7 Higgs fields. 

In contrast to the "old" extended AGUT by Froggatt-Nielsen (called here 
as Gf-theory), the new results of G e xt _ theory are more encouraging, and it is 
possible to conclude that the G-theory, in general, is successful in describing of 
the SM experiment. 



3 Multiple Point Principle 

AGUT approach is used in conjunction with the Multiple Point Principle pro- 
posed in Ref.[^. According to this principle Nature seeks a special point — the 
Multiple Critical Point (MCP) — which is a point on the phase diagram of the 
fundamental regulirized gauge theory G (or Gf, or G ex t)/ where the vacua of all 
fields existing in Nature are degenerate having the same vacuum energy den- 
sity. Such a phase diagram has axes given by all coupling constants considered in 
theory. Then all (or just many) numbers of phases meet at the MCP. 

MPM assumes the existence of MCP at the Planck scale, insofar as gravity 
may be "critical" at the Planck scale. 

The usual definition of the SM coupling constants: 

ai = 3 — T5 — - a 2 = — — n — » «3 = a s = — , (8) 
3 cos 2 sin 2 e^rs 47t 

where a and a s are the electromagnetic and SU(3) fine structure constants, re- 
spectively, is given in the Modified minimal subtraction scheme (MS). Here Qj^s 
is the Weinberg weak angle in MS scheme. Using RGE with experimentally es- 
tablished parameters, it is possible to extrapolate the experimental values of three 
inverse rurining constants oc^ 1 (here [i is an energy scale and i=l,2,3 corre- 
spond to U(l), SU(2) and SU(3) groups of the SM) from the Electroweak scale to 
the Planck scale. The precision of the LEP data allows to make this extrapolation 
with small errors (see [||]). Assuming that these RGEs for ct^ 1 contain only the 
contributions of the SM particles up to |J. « [ipi and doing the extrapolation with 
one Higgs doublet under the assumption of a "desert", the following results for 
the inverses \ 3 (here oty = §<Xi ) were obtained in Ref.||| (compare with [^]): 

o^ 1 (ji Pl )w55.5; ocjVpO ~49.5; oc^ 1 (u Pl ) « 54.0. (9) 
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The extrapolation of <Xy \ U P to the point |x = upi is shown in Fig.l. 
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Fig.l: The evolution of three inverse running constants (/^). The 
extrapolation of their experimental values from the Elec- 
troweak scale to the Planck scale was obtained by using the 
renormlization group equations. AGUT works in the region 

< A* < KPi- 



According to AGUT, at some point jx = |j_g < M-pi the fundamental group G 
(or Gf, or G ex t) undergoes spontaneous breakdown to its diagonal subgroup: 

G > Gdiag.subgr. = {g, g, 9 1 1 9 S SMG}, (10) 

which is identified with the usual (low-energy) group SMG. The point \Iq ~ 1 1 8 
GeV also is shown in Fig.l, together with a region of G-theory, where AGUT 
works. 

The AGUT prediction of the values of <Xi(|i) at \i = \ip\. is based on the MPM 
assumption about the existence of phase transition boundary point MCP at the 
Planck scale, and gives these values in terms of the corresponding critical cou- 
plings oti.crtt til, 8, §]: 

fti(M-pi) — °^ |Crit _ tXt ^ rtt f or i — 2 y 3, (also for i > 3), (11) 

I i a 1,CTit crit c iimi /n\ 

2 N gen(N g en + 1) ° 

4 Lattice Theories 



The philosophy of MPM leads to the necessity to investigate the phase transi- 
tion in different gauge theories. A lattice model of gauge theories is the most 
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convenient formalism for the realization of the MPM ideas. As it was mentioned 
above, in the simplest case we can imagine our space-time as a regular hyper- 
cubic (3+l)-lattice with the parameter a equal to the fundamental (Planck) scale: 
a = Ap. 

The lattice SU(N) gauge theories was first introduced by K.Wilson [ |l9| ] for 
studying the problem of confinement. He suggested the following simplest ac- 
tion: 

p 

where the sum runs over all plaquettes of a hypercubic lattice and U v belongs to 
the fundamental representation of SU(N). The simplest Wilson lattice action for 
U(1 ) gauge theory has the form: 

S w = (3 Y_ cos p » where U v = e i@v . (14) 
v 

For the compact lattice QED: (3 = 1 /e^, where eo is the bare electric charge. 
The Villain lattice action for the U( 1 ) gauge theory is: 

S v = (|3/2)^(e p -27tk) 2 , keZ. (15) 
v 

The critical value of the effective electric fine structure constant a was ob- 
tained in Ref.[|20|] in the compact QED described by the Wilson and Villain actions 
dl4| ) and (fT5|), respectively: 

„lat n -)A _l_ A A1 K si I at 



a£?t m 0.20 ± 0.015 and 1 .25 ± 0.10 at p T = Pent « 1-011. (16) 

e g 2 

u a. — — and a = — (17) 

Here 4n An y ' 

are the electric and magnetic fine structure constants, containing the electric charge 
e and magnetic charge g . 

The result of Ref.[^| for the behavior of cx((3) in the vicinity of the phase 
transition point |3y is shown in Fig.2(a) for the Wilson and Villain lattice actions. 
Fig.2(b) demonstrates the comparison of the function <x(|3) obtained by Monte 
Carlo method for the Wilson lattice action and by theoretical calculation of the 
same quantity. The theoretical (dashed) curve was calculated by so-called "Parisi 
improvement formula" [pT[: 



a(P) = [47tpWp]-'. (18) 

Here W p =< cos P > is a mean value of the plaquette energy. The correspond- 
ing values of W p are taken from Ref.p^]. 

The theoretical value of a cr it is less than the "experimental" (Monte Carlo) 
value ©: 

a CT it(in lattice theory) w 0.12. (19) 



According to Fig.2(c) 

crit. ,theor 



<*cr\t..theor. « 8. (20) 
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Fig.2(a,b): (a) The renormalized electric fine structure constant 
plotted versus f3 / fir for the Villain action (circles) 
and the "Wilson action (crosses). The points are ob- 
tained by the Monte— Carlo simulations method for 
the compact QED; 

(b) The behavior of the effective electric fine struc- 
ture constant in the vicinity of the phase transition 
point obtained with the lattice "Wilson action. The 
dashed curve corresponds to the theoretical calcu- 
lations by the "Parisi improvement method". 



This result does not coincide with the lattice result (p^|) / which gives the following 
value: 

acrtt-.lat." 5 - (21) 

The deviation of theoretical calculations of <x(6) from the lattice ones, which is 
shown in Fig.2(b,c), has the following explanation: Parisi improvement formula 
( jl8| ) is valid in Coulomb-like phase where the mass of artifact monopoles is in- 
finitely large and photon is massless. But in the vicinity of the phase transition 
(critical) point the monopole mass m. — > and photon acquires the non-zero 
mass rao 7^ (on side of the confinement). This phenomenon leads to the "freez- 
ing" of a: the effective electric fine structure constant is almost unchanged in the 
confinement phase and approaches to its maximal value a = a max . The authors 




Fig. 2c: The behavior of the inverse effective electric fine structure 
constant in the vicinity of the phase transition point plot- 
ted versus fB / /Sjp for the simple Wilson lattice action. The 
dashed curve corresponds to the theoretical calculations 
by the "Parisi improvement method". 



of Ref.[p3|] predicted that in the confinement phase, where we have the forma- 
tion of strings, the fine structure constant a cannot be infinitely large, but has the 
maximal value: oc mQX w ~ 0.26 due to the Casimir effect for strings. 



5 Lattice Artifact Monopoles and Higgs Monopole Model 

Lattice monopoles are responsible for the confinement in lattice gauge theories 
what was confirmed by many numerical and theoretical investigations. 

In the previous papers [[§]-[0]] the calculations of the U(l) phase transition 
(critical) coupling constant were connected with the existence of artifact monopo- 
les in the lattice gauge theory and also in the Wilson loop action model [|l|] . 

In Ref.[llJ we (L.V.L. and H.B.Nielsen) have put forward the speculations of 
Refs.[|||,[|lC ]] suggesting that the modifications of the form of the lattice action 
might not change too much the phase transition value of the effective continuum 
coupling constant. The purpose was to investigate this approximate stability of 
the critical coupling with respect to a somewhat new regularization being used 
instead of the lattice, rather than just modifying the lattice in various ways. In [ pl[ ] 
the Wilson loop action was considered in the approximation of circular loops of 
radii R > a. It was shown that the phase transition coupling constant is indeed 
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approximately independent of the regularization method: a cr it ~ 0.204, in corre- 
spondence with the Monte Carlo simulation result on lattice: <x cr u « 0.20 ± 0.015 
(see Eq.@). 

But in Refs.[|@I-[0] / instead of using the lattice or Wilson loop cut-off, 
we have considered the Higgs Monopole Model (HMM) approximating the lat- 
tice artifact monopoles as fundamental pointlike particles described by the Higgs 
scalar field. The simplest effective dynamics describing the confinement mech- 
anism in the pure gauge lattice U(l) theory is the dual Abelian Higgs model of 



scalar monopoles [|24J, [[[12j-[|lJ]] (shortly HMM). This model, first suggested in 



Refs.[}24j|, considers the following Lagrangian: 

L = -^l I F^(B) + l|(3 tl -iB ti )0| 2 -U(0) ) where U(O) = ^ 2 |0| 2 + ^|0| 4 

(22) 

is the Higgs potential of scalar monopoles with magnetic charge g, and B H is the 
dual gauge (photon) field interacting with the scalar monopole field O. In this 
model A is the self-interaction constant of scalar fields, and the mass parameter 
[i 2 is negative. In Eq.(^) the complex scalar field O contains the Higgs (cf)) and 
Goldstone (x) boson fields: 

= 4>+i x . (23) 
The effective potential in the Higgs Scalar ElectroDynamics (HSED) was first cal- 



culated by Coleman and Weinberg [ 25 1 in the one-loop approximation. The gen- 
eral method of its calculation is given in the review [g6[|. Using this method, we 
can construct the effective potential for HMM. In this case the total field system of 
the gauge (B H ) and magnetically charged (O) fields is described by the partition 
function which has the following form in Euclidean space: 

Z= [DB][DO][D(P + ]e- s , (24) 

where the action S = J d 4 xL(x) + S g f contains the Lagrangian ( p2| ) written in 
Euclidean space and gauge fixing action S g f . Let us consider now a shift: 

(D(x) = <5 b + 6(x) (25) 

with Ob as a background field and calculate the following expression for the 
partition function in the one-loop approximation: 



Z = 



[DB] [DO] [DO+] exp{-S(B, O b ) - 



d 4 x[^lo^ b a(x)+h.c.]} 



exp{-F(O b) g 2 ,^,A)}. (26) 



Using the representation (23), we obtain the effective potential: 

Veff = F(cb b ,g 2 ,H 2 ,A) (27) 



given by the function F of Eq.(g6|) for the constant background field Ob = <$>b — 
const. In this case the one-loop effective potential for monopoles coincides with 
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the expression of the effective potential calculated by the authors of Re f.Pq ] for 
scalar electrodynamics and extended to the massive theory (see review p6|]). 

Considering the renormalization group improvement of the effective Coleman- 
Weinberg potential [^],|^] ; written in Ref.[fl4y for the dual sector of scalar elec- 
trodynamics in the two-loop approximation, we have calculated the U(l) critical 
values of the magnetic fine structure constant: 

*crit = 9crit/47t « 1 .20, (28) 

and electric fine structure constant: 

<*crit = ^Tit « °- 208 

by the Dirac relation: 

eg — 2n, or aa = -. (30) 
The values (28),(p9[) coincide with the lattice result (|l6|). 



6 Monopoles strength group dependence 

As it was shown in a number of investigations, the confinement in the SU(N) 
lattice gauge theories effectively comes to the same U(l) formalism. The reason 
is the Abelian dominance in their monopole vacuum: monopoles of the Yang- 
Mills theory are the solutions of the U(l)-subgroups, arbitrary embedded into 
the SU(N) group. After a partial gauge fixing (Abelian projection by 't Hooft |Q) 
SU(N) gauge theory is reduced to an Abelian U ( 1 ) N ~~ 1 theory with N — 1 different 
types of Abelian monopoles. Choosing the Abelian gauge for dual gluons, it is 
possible to describe the confinement in the lattice SU(N) gauge theories by the 
analogous dual Abelian Higgs model of scalar monopoles. 

Considering the Abelian gauge and taking into account that the direction 
in the Lie algebra of monopole fields are gauge independent, we have found 
in Ref . fllil ] an average over these directions and obtained the group dependence 



relation between the phase transition fine structure constants for the groups U( 1 ) 
andSU(N)/Z N : 

n PnTT _! 

a N,crit - yy N _ i a U(1),crif V L ) 

7 AGUT-MPM prediction of the Planck scale values of the U(l), 
SU(2) and SU(3) fine structure constants 

As it was assumed in Ref.[^], the MCP values <Xi. jCT i t in Eqs.© and @ coin- 
cide with the critical values of the effective fine structure constants given by the 
generalized lattice SU(3)-, SU(2)- and U(l)-gauge theories. 

Now let us consider ct^ 1 (~ oc -1 ) at the point [i — [Iq ~ 10 18 GeV shown in 
Fig.l. If the point |^ = |j.g is very close to the Planck scale |J. = [ipi, then according 
to Eqs.(||) and (|l^), we have: 

a -i wa -i raa -i ~ g^Vcj , . 

"-1st gen. ~ a 2nd gen. ~ a 3rd gen. ~ r ~ 7 ' \ OA > 
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what is almost equal to the value 



obtained by the Parisi improvement method (see Fig.2(c)). This means that in the 
U(l) sector of AGUT we have a near the critical point. Therefore, we can expect 
the existence of MCP at the Planck scale. 

It is necessary to mention that the lattice investigators were not able to obtain 
the lattice triple point values at iCr tt (i=l,2,3 correspond to U(1),SU(2) and SU(3) 
groups) by Monte Carlo simulation methods. These values were calculated theo- 
retically by Bennett and Nielsen in Ref.^ using the Parisi improvement method 

«Y,cHt«9-2±1> < CTit «16.5±1, ai; 1 CTtt wl8.9±l. (33) 
Assuming the existence of MCP at \i = [ipi and substituting the last results in 
Eqs.(||) and (Q, we have the following prediction of AGUT [||: 

cXyVpi) «55±6; oq 1 (u. Pl ) » 49.5 ± 3; oq 1 (u Pl ) » 57.0 ± 3. (34) 

These results coincide with the results (^|) obtained by the extrapolation of exper- 
imental data to the Planck scale in the framework of pure SM (without any new 
particles) [§, §. 

Using the relation ([31]), we obtained the following relations: 



^2,crit ■ ^s.crit = 1 : : 3/V2 = 1 : 1 .73 : 2.12. (35) 

Let us compare now these relations with the MPM prediction. 
For oq \ rit ~ 9.2 given by the first equation of (33), we have: 

"Y.lrit : «2,cHt = o^Lrtt = 9.2 : 15.9 : 1 9.5. (36) 

In the framework of errors the last result coincides with the AGUT-MPM predic- 
tion (|33"|). Of course, it is necessary to take into account an approximate descrip- 
tion of confinement dynamics in the SU(N) gauge theories, which was used in 
our investigations. 



8 The possibility of the Grand Unification Near the Planck 
Scale 

We can see new consequences of the extension of G -theory, if G -group is broken 
down to its diagonal subgroup Gdiag/ i.e. SM, not at |j.g ~ 10 18 GeV, but at [Iq ~ 
10 15 GeV. In this connection, it is very attractive to consider the gravitational 
interaction. 



8.1 "Gravitational finestructure constant" evolution 

The gravitational interaction between two particles of equal masses M is given 
by the usual classical Newtonian potential: 

,2 / w x 2 



M 1 _ t M y 1 a g (M) 
r \M Pl y t 



Vg=-G— =- — r = -^T^. (37) 
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2 4 6 8 10 12 14 16 18 20 



log, w n(GeV) 

Fig. 3: The intersection of the inverse "gravitational finestructure 
constant" ctg 1 ^) with a^ 1 (fi) occurs at the point (x , a^ 1 ): 

x « 18.3, cXq 1 « 34.4, 
where x = log 10 fi(GeV) . 



which always can be imagined as a tree-level approximation of quantum gravity. 
Then the quantity: 




plays a role of the "gravitational finestructure constant" and the evolution of its 
inverse quantity is presented in Fig.3 together with the evolutions of oC{\ 3 (u) 
(here we have returned to the consideration of <x\ instead of Oiy). 

Then we see the intersection of oc.g 1 ( [i) with a^ 1 ( \i) in the region of G-theory 
at the point: 

where 

x «18.3, and cc^ 1 w 34.4, (39) 

and x — log 10 \i. 
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8.2 The consequences of the breakdown of G-theory at Ug ~ 1 1 5 or 1 1 6 
GeV 

Let us assume now that the group of symmetry G undergoes the breakdown to 
its diagonal subgroup not at ug ~ 10 18 GeV, but at u-g ~ 10 15 GeV, i.e. before the 
intersection of 1 and a J 1 (u) at \i - 10 16 GeV. 

As a consequence of behavior of the function otT 1 ( |3 ) near the phase tran- 
sition point, shown in Fig.2c, we have to expect the change of the evolution of 
a^ 1 (u.) in the region u. > |j.g shown in Fig.l by dashed lines. Instead of these 
dashed lines, we must see the decreasing of ct^ 1 ( when they approach MCP, if 
this MCP really exists at the Planck scale. 

According to Fig.2c, in the very vicinity of the phase transition point (i.e. 
also near the MCP at u. = \i-vx), we cannot describe the behavior of a^ 1 (u.) by the 
one — loop approximation RGE. 

It is well known, that the one-loop approximation RGEs for oc^ 1 (u.) can be 
described in our case by the following expression 

arV^arVp^ + ^log^), (40) 

where bt are given by the following values: 

bt = (bi,b 2 ,b 3 ) = 



I 4N 9^ 1 K, 22 M 4N 9^ 1 X . 11K . 4Ng^ .... 

( 1 ttNs, t N v 1 -N s , 11N V j—). (41) 

5 I v 3 5 b 5 

The integers N gen , Ns, Ny are respectively the numbers of generations, Higgs 

bosons and different vector gauge fields of given "colors". 

For the SM we have: 

N gen = 3, N s =Nv = 1, (42) 

and the corresponding slopes describe the evolutions of oc^ 1 (u) up to u = |j.g 
presented in Fig.3. 

But in the region < < HPi/ when G-theory works, we have Ny = 3 
(here we didn't take into account the additional Higgs fields which can change 
the number N s ), and the one-loop approximation slopes are almost 3 times larger 
than the same ones for the SM. In this case, it is difficult to understand that such 
evolutions give the MCP values of a^ 1 (|Xpi), which are shown in Fig.4. These 
values were obtained by the following way: 

3 

°hVpi)~6--a u | 1)crit «13, 
ajVpi) «3- V3a u J 1) crU » 19, 



3_ 



(m )~ 3 -— a u 1)iCrit «24, (43) 
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The case of the breakdown of G-theory at [la- The evolution 






of 


a 1 2 3(m) decreases near the Planck scale and approaches 






to 


the Multiple Critical Point at fj, = It 


is shown the 



possibility of unification of all gauge interactions, includ- 
ing gravity, at O^GUT ^ 34.4 and Xgut S3 18.3, where 
x = log 10 fj,{GeV): 

a) no S3 10 16 (GeV), v S3 1; 
6) /t G S3 10 15 (GeV), S3 2.4. 



where we have used the relation (^lj) with 



aum.crit = — -tk— ~ 0.77a CTit , (44) 
cos 2 



taking into account our HMM result (|29(): a cr i t ~ 0.208, which coincides with the 
lattice result (|^) and gives: 

inherit « 3.7. (45) 

In the case when G-group undergoes the breakdown to the SM not at ~~ 1 
GeV, but at (J-g ~ 10 GeV, the artifact monopoles of non-Abelian SU(2) and 
SU(3) sectors of G-theory begin to act more essentially. 

According to the group dependence relation ( |3l| ) (although now it is neces- 
sary to expect that it is very approximate) we have, for example, the following 
estimation at u_g ~ 10 15 GeV: 

<*u(d(M-g)~7 - for SU(3)i st gen., etc. (46) 

which is closer to MCP than the previous value ocy 1 ~ 9, obtained for the AGUT 
breakdown at [ic ~ 10 18 GeV. 
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It is possible to assume that |3-functions of SU(2) and SU(3) sectors of G- 
theory change their one-loop approximation behavior in the region jo. > 10 16 
GeV and 3 ( (x) begin to decrease, approaching the phase transition (multiple 
critical) point at \i = |J.pi.. This means that the asymptotic freedom of non-Abelian 
theories becomes weaker near the Planck scale, what can be explained by the in- 
fluence of artifact monopoles. It looks as if these (3-functions have singularity at 
the phase transition point and, for example, can be approximated by the follow- 
ing expression: 

dor 1 (3(a) , ,, a , . . . 

= ss A(l ) near the phase transition pomt. (47) 

dt a a CT it 

This possibility is shown in Fig.4 for "v « 1 and ~v « 2.4. 

Here it is worth-while to comment that such a tendency was revealed in the 
vicinity of the confinement phase by the forth-loop approximation of |3-function 
in QCD (see Ref.[|9|). 



8.3 Does the [SU(5)] 3 SUSY unification exist near the Planck scale? 

Approaching the MCP in the region of G-theory (|j.g < H < M-pi)/ cx^~ 3 ( M-) show 
the necessity of intersection of 1 (|i) with 1 at some point of this region 
if |j.g ~ 10 15 or 10 16 GeV (see Fig.4). If this intersection takes place at the point 



Xq, 0C ) given by Eq.(39), then we have the unification of all gauge interactions 



(including the gravity) at the point: 

(xgut; a GUT ) « (18.3; 34.4), (48) 

where x = log 10 |j.(GeV). Here we assume the existence of [SU(5)] 3 SUSY unifica- 
tion having superparticles of masses 

Mwl0 183 GeV. (49) 



The scale Hgut = M, given by Eq.(^9J), can be considered as a SUSY breaking 
scale. 

Figures 5(a,b) demonstrates such a possibility of unification. We have inves- 
tigated the solutions of joint intersections of oCg 1 (u.) and all a^ 1 (u.) at different 
xgut with different v in Eq.(^). These solutions exist from v w 0.5 to "V « 2.5. 

The unification theory with [SU(5)] 3 -symmetry was suggested first by S. Raj- 
poot@. 

It is essential that the critical point in this theory, obtained by means of 
Eqs.(|Tl|), ( pl[ ) and (p5|), is given by the following value: 

*5,crit ~ J £ V 2 



a 5 ' «3-- A /-« 34.0. (50) 



The point (50) is shown in Fig.5(a,b) presented for the cases: 

1. -v « 1, Uquj « 34.4, xgut ~ 18.3 and [i G ~ 10 16 GeV shown in Fig.5(a); 

2. -v « 2.4, ag UT w 34.4, x G ut ~ 18.3 and |i G « 10 15 GeV shown in Fig.5(b). 



Phase Transition in Gauge Theories and the Planck Scale Physics 97 



a* 



a* 



50 


1 1 1 ] r ■ 

I 


1 1 - 
1 


40 

30 


"2_______ 

a s\^- — 




a GUT 


1 
1 

1 ' 
1 


20 






GCv 1 


1 
1 
1 
1 


10 


1 






1 

I , 




10 12 


14 16 


18 


/'• f*i 20 


50 


1 1 

Q f 1 


1 1 
1 


TT 


i 1 - 
i 

1 (b) 

1 
1 


40 


a^J____ ~ ™ : 


1 ■ — . 
1 




30 




1 
1 


1 
1 


20 




1 
1 
1 
1 


CKcr 1 


1 
1 
1 
1 


10 




1 
1 




1 

1 1 I 




10 12 


14 CO 16 

iogio(i(Gey) 


IS 


20 


.5: The possibility of the 


[SU(5)] 3 SUSY 


unification of all gauge 



interactions, including gravity. A dashed line corresponds to 
the evolution of ck 5 approaching to the critical point at 

«F 1 (A«P!) ~ 34.0: 



a) /i G ks 10 16 (GeV), v ss 1; 
6) /i G « 10 15 (GeV), 1/ ss 2.4. 



We see that the point ( pCj ) is very close to the unification point ctgjjj ~ 
34.4 given by Eq.(^). This means that the unified theory, suggested here as the 
[SU(5)] 3 SUSY unification, approaches the confinement phase at the Planck scale. 
But the confinement of all SM particles is impossible in our world, because then 
they have to be confined also at low energies what is not observed in the Nature. 

It is worth-while to mention that using the Zwanziger formalism for the 
Abelian gauge theory with electric and magnetic charges (see Refs.[pl|]-[p3[|] and 
the possibility of unification of all gauge interactions at the Planck scale was 
considered in Ref . [fl8|] in the case when unconfined monopoles come to the exis- 
tence near the Planck scale. They can appear only in G-theory, because RGEs for 
monopoles strongly forbid their deconfinement in the SM up to the Planck scale. 
But it is not obvious that they really exist in the G-theory. This problem needs 
more careful investigations, because our today knowledge about monopoles is 
still very poor. 

The unified theory, suggested in this paper, essentially differs in its origin 
from the case considered in Ref .[^], because this theory does not assume the exis- 
tence of deconfining monopoles up to the Planck scale, but assumes the influence 
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of our space-time lattice artifact monopoles near the phase transition (critical) 
point. 

Considering the predictions of this unified theory for the low-energy physics 
and cosmology, maybe in future we shall be able to answer the question: "Does 
the [SU(5)] 3 SUSY unification theory really exist near the Planck scale?" 

9 Reminiscences on a collaboration with Holger 

I got acquainted with Holger Bech Nielsen in the year 1970, during the Roches- 
ter Conference in Kiev (Ukraineq, USSR). His brilliant achievements and childish 
charm conquered me. 

When I returned to Tbilisi after this conference, I received from Holger a very 
kind letter, a pretty present and 9 papers of his activity. Then I understood that 
the destiny gave me a chance to meet a genius. 

That time my homeland of Georgia entered the USSR, and KGB immediately 
called me to their office and began to cross-examine me, and shout: "Why 'Dear 
Lara' in his letter, but not 'Dr. Laperashvili'?! You will be thrown out of the Insti- 
tute immediately! Be glad that you will be not imprisoned!" Then I did not inform 
Holger about this story, because I was afraid that he would stop writing to me. 
The director of my Institute of Physics in Tbilisi saved me: he told to KGB, that I 
am spy and I must find out a very important information (it was lie, of course). 
And Holger continued to write his letters and to send his papers. Six times he 
invited me to NBI, but without success: I was forbidden to go abroad and to col- 
laborate with him. And only after "perestroyka", in 1987, 1 began to collaborate 
with him and came to Copenhagen first in 1990. After 20 years! But our friendship 
is alive up to now. 

Dual strings, analogue model, Koba-Nielsen representation, KNO-scaling, 
ANO-strings, spaggheti vacuum, Cheshire cat bag model, Nielsen-Ninomiya 
theorem, Random dynamics, Multiple Point Principle, Anti-grand unification, 
Higgs-Top masses prediction, etc.,etc.,etc, - Holger 's name is well-known in the 
World of Physics. He is in the first row of the World physicists. There was time 
when all leading (high energy) physical journals were full of his name in the titles 
of different papers. 

I am sure that Holger made for Denmark so much that his name and activity 
will add a glorious golden page to history of Denmark. And of course, Russia, 
Georgia and the World Science will not forget his outstanding services to physics. 
I am sure that he will produce excellent physics in the nearest future, because he 
is fresh, active and full of non-trivial ideas. 

The collaboration with Holger Bech Nielsen was the best time of my life. 
Always busy with physics, he has unimaginable flux of energy, fantastic intuition 
and a great volume of knowledges. He operates by complicated mathematical 
structures. 

Holger is a poet in physics. 

Simultaneously, he has the severity in physical investigations and childish 
perception of human life. He is a bright person, full of kindness and progressive 
ideas. Many people in different countries love him. 
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Dear Holger, I congratulate you for your 60th jubilee and wish you many, 
many years activity in future, health, wealth and a Great Success. Thank you 
very much for collaboration with me. I send you my Great Love. 
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Some Remarks on the "Classical" Large N Limit 
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Abstract. It has been proposed some time ago that the large N —limit can be understood as 
a "classical limit", where commutators in some sense approach the corresponding Poisson 
brackets. We discuss this in the light of some recent numerical results for an SU(N) gauge 
model, which do not agree with this "classicality" of the large N —limit. The world sheet 
becomes very crumpled. We speculate that this effect would disappear in supersymmetric 
models. 



In the 1970ies Holger and I discussed the confinement problem quite a lot. 
In those days the picture of confinement by means of a string connecting the 
quarks became more and more accepted. It was known that the Nambu-Goto 
action could be written as a square root of the Poisson bracket squared. One of 
the subjects we discussed was whether somehow the Poison bracket was related 
to the commutator like in quantum mechanics, in some limit. However, we did 
not gain any insight in this. After the work by Eguchi and Kawai it became of 
course clear that the relevant limit was the large N limit, where the action could 
be expressed in terms of commutators. Later on Hoppe [[!) constructed a very 
suggestive relation between Poisson brackets and commutators. This was then 
used in ingenious ways by a number of authors in connection with the 

Eguchi-Kawai formulation of the large N theory. 

To give a brief review of this approach, let us consider SU(N) with the gauge 
field A H and the generators Ik, k = (ki , k2), 

(A^=^< (l k )>, (1) 

k 

where are expansion coefficients to be integrated in functional integrals. This 
expansion can be compared to the one for a string variable X,j.(cr, x), 

X h (ct,t) =2>£e ikff , (2) 

k 

with cr = (cr, x). The variables X H are the Weyl transform of the matrices A^. The 
generators can be constructed explicitly in terms of the Weyl matrices, and have 
the commutation relation 

N /2n \ 
[Ik, ,lk 2 ] sm — k! x k 2 lk,+k 2 -»i(ki x k 2 ) lk,+k 2 (3) 



2tt V N 
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for N — > oo, where ki x k 2 — (ki )i (^2)2 — (k2)i (ki )z- This leads to a comparison 
of the commutator with the Poisson bracket. For example, one has 

[A H , A v ] = t— Y_ a H a ^ sin ( x P ) l k, +k 2 , (4) 

k,p 



and 



Using 



{X^(cr),X„(cx)} = iY_ a$ a? (k x p) e 1 '^'. (5) 
k.P 



tr Iklp cx 6k+ P ,o. 

one can derive that 



d 2 cje i(k+ P )<J oc6 k+p , 0) (6) 



tr [A,, A v f -> const. 



d 2 a{X H ,X,} 2 , (7) 



where the constant (depending on N) can be removed by a suitable normaliza- 
tion, and where it was used that 

^ sm ( Tr^i x k 2 J — > ki x k 2 (8) 



2tt V n 

for N — > 00, even inside the relevant sum over modes. In this sense the commuta- 
tor approaches the Poisson bracket [||,[§||]. 

Unfortunately this argument, leading to a nice connection between strings 
and fields in the large N— limit, depends on performing the limit N — > 00 inside 
sums over modes like in eq. This is evidently valid if the low modes dominate, 
since the difference between sin x and the linear function x is considerable when 
x is not close to 0. 

Recently Anagnostopoulos, Nishimura and the author [^] started to investi- 
gate numerically the mode-distribution for large N in a four-dimensional SU(N) 
model with the partition function 

dA ^ exp ^J^trlA^A^ 2 ^ . (9) 



This model has been shown to exist The expectation value of the commu- 

tator squared 

M=<tr[A^A v ] 2 >, (10) 

can be found from a scaling argument ^ to be N 2 — 1, and we checked our 
numerical method by seeing that it leads very precisely to this result. 

Then we computed the corresponding expectation value of the Poisson bracket 
by introducing the Weyl-transform[]X H (cF) (defined on a torus) of A^ and we 
therefore found the behavior of 



< P >=< 



d 2 CT {Xn(ff), X. (a)} 2 >. (11) 



1 The Weyl transform X^cr) of the matrix A H is given by X^(ff) oc y^ k exp(ikcr) tr(l k A,J 
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Here the normalization is such that if only the very lowest modes are kept, one 
has the result < M >=< P >. For the details of this construction, we refer to the 
paper [Q]. One important point discussed in this paper is that the Poisson bracket 
is gauge dependent, and hence < P > is in general gauge dependent, in contrast 
to the gauge invariant < M >. We have selected a gauge fixing analogous to the 
Landau gauge. Of course, one can anyhow argue that if the commutator really 
approaches the Poisson bracket, gauge invariance of the former implies gauge 
invariance of the latter, at least approximatively 

In our numerical calculations we took N = 15,25, and 35. The results were 
quite embarrasing. Keeping only a few modes, we got < M >«< P >, but includ- 
ing all modes these two quantities differ considerably: Whereas very precisely we 
find <M>=N 2 -1,it turns out that < P > rather grows like 0(N 4 ). Therefore 
the two expectation values do not at all agree, and the discrepancy actually in- 
creases with N. We can therefore say that in the model with the action (^), the 
Poisson bracket is not approximated by the commutator. 

We also computed the average of the area of the world sheet, 



< A >=< 



a^{X^(o-),X,(o-)}2 > . (12) 



It turns out that < A >«< M >. Thus, the expectation value of the Nambu- 
Goto action does indeed behave approximately as the commutator. This does, 
however, not show that the action (|9|) is approximately equal to the Nambu-Goto 
action, since we have computed only one moment of A. 

The negative result given above shows that the question of whether the large 
N limit is a classical limit, depends on the dynamics, and a simple model like (g) 
is not able to suppress the higher modes, and hence one cannot use the limit 
(^) inside sums over modes. The string picture resulting from the behavior of 
Xfxla) turns out to be an extremely crumpled string, with an almost infinite (« 
33) Hausdorff dimension. This is another way of seeing that higher modes are 
important, since dominance of lower modes would lead to a smooth string. 

As pointed out in [^] the model @ does exist in a world-sheet version, since 
the commutator can be replaced by star products of the X's, so the action can be 
rewritten as 





' 1 


dX H exp 


4g 2 N . 



d 2 a (X H ( a) * X v ( a) - X^ ( cx) * X ^ (a) Y 



(13) 



which exists because the corresponding matrix model exists. The world sheet 
defined by this theory (with an infinite number of derivatives in the star product) 
is the same as before, and is thus extremely crumpled. 

Thus the dynamical question of the existence of a "classical" large N— limit 
boils down to finding a rather smooth string, since any "violent" string would re- 
quire higher modes. Presumably the QCD string is smooth. In this case the action 
(^) is replaced by a quenched action, which hopefully can manage to suppress the 
high modes. Another possibility is to take a supersymmetric version of (||), since 
supersymmetry in general makes strings more smooth: For example, in the vac- 
uum energy the tachyon is removed by supersymmetry. For the simple bosonic 
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string it was shown long time ago by Alvarez [||] that the bosonic string is very 
crumpled, and by the author [ [l0| ] that the tachyon produces this crumpling. Al- 
varez' formula for the free energy [pfl, 



2not' 



WR2-R2, (14) 



can actually be interpreted [ fjo| ] as a formula for the lowest tachyonic energy, with 
the tachyon mass squared being proportional to — R§. The phenomenon of crum- 
pling should thus not occur in the supersymmetric case. We hope that these ques- 
tions will be investigated numerically for supersymmetric versions of (^) in the 
future. 

Finally, I would like to take this opportunity to thank Holger for pleasant 
years of discussions of very many aspects of physics. 
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Abstract. We begin by outlining the ancient puzzle of off shell currents and the infinite 
size particles in a string theory of hadrons. We then consider the problem from the mod- 
ern AdS/ CFT perspective. We argue that although hadrons should be thought of as ideal 
thin strings from the 5-dimensional bulk point of view, the 4-dimensional strings are a 
superposition of "fat" strings of different thickness. 

We also find that the warped nature of the target geometry provides a mechanism for 
taming the infinite zero point fluctuations which apparently produce a divergent result for 
hadronic radii. 

1 Meeting Holger 

When I was a school kid during the early 1950's we used to have to read a mag- 
azine called "The Reader's Digest". It was full of corny articles about patriotic 
platitudes which were very boring but it always had an interesting section called 
"My Most Unforgettable Character ". It was usually about a somewhat eccentric 
but admirable character that the writer had especially fond memories of. Well, for 
me (and I suspect anyone else who knows him), Holger will always be one of the 
most unforgettable characters I've ever met. 

Holger and I first met through the mail in 1970. He had seen a paper that 
I wrote claiming that the Veneziano amplitude described the scattering of some 
kind of elastic strings. Unfortunately I no longer have the hand written letter but 
I can still see his distinctive curly handwriting and the signature - Holger Bech 
Nielsen. Most of all I remember his almost child-like enthusiam and simplicity. 
He too had been working on a similar idea [J Unlike so many messages that I've 
received over the years, this one had nothing to do with staking a claim or as we 
say, pissing on territory. The letter straightforwardly expressed his excitement 
and joyously shared his own ideas. It was completely clear to me that I had met 
a larger than life, most unforgettable character. 

That year I invited Holger to spend a month visiting me in New York and 
what a month it was. We ate too much, drank too much and yelled too loud but 
the physics excitement was palpable. I have never had more fun doing physics 
than during that time. At the time, string theory was of course a theory of hadrons. 
Mesons were strings with quarks at their ends. Both of us were disturbed by 

* SU-ITP 00-25 

1 Nambu had also been working on the same ideas but I don't believe that Holger or I 
were aware of it. 
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something that we thought was a very serious shortcoming of the theory. At 
about that time the electromagnetic properties of hadrons were under intensive 
investigation at SLAC and other places. Electromagnetic form factors of nucleons 
were already well measured. SLAC had measured deep inelastic electroproduc- 
tion and Feynman and Bjorken had explained the data with their parton ideas. 
The problem that puzzled Holger and me was that we could see no way to define 
the local electromagnetic current of a hadron using string theory. Every time we 
tried we got nonsense. Hadrons came out infinitely big and in a sense, infinitely 
soft. Holger and I had a wonderful time thinking about the problem. I think it 
is fair to say that many of the themes of my future work trace back to that brief 
month and to Holger 's profound influence on me . 

What was our solution to the problem? Both of us were inclined to think of 
the string as an idealized limit of a discrete system. In my case I viewed it as a sys- 
tem of partons in the light cone frame. Holger had a more covariant view which 
he had been discussing with Aage Kraemmer in Copenhagen. According to this 
view, the Koba Nielsen disc [[[]] was really the continuum limit of an infinitely 
dense planar Feynman diagram or more precisely, a sum over such diagrams. At 
that time we had no idea why planar diagrams should dominate. That had to 
wait for 't Hooft [|]. 

One of the ideas in the paper was that the geometry of a large planar dia- 
gram defined a kind of metric which could be gauge fixed to what would now be 
called the conformal gauge. We realized that a proper treatment should include a 
sum over diagrams which could be represented as a path integral over a diagram 
density. Holger wanted to treat this degree of freedom as an additional dimen- 
sion which a decade later, following the work of Polyakov, became the Liouville 
field. I thought it was a bad idea since I could not see how an additional infinite 
direction could fit into hadron physics. For this reason we decided that the inte- 
gration should be dominated by some specific density that didn't fluctuate too 
much. 

In fact the form factor problem forced us to conclude that the continuum 
limit was just too extreme. Hadrons might be described by fairly dense systems 
of partons but not a continuum. There had to be a cutoff which limited the zero- 
point fluctuations that blew the string up to infinite size and also removed 
the hard effects of discrete partons. Together with Kraemmer we wrote a paper 
[^] formulating what we called the Dual Parton Model which tried to keep the 
good features of strings without passing to the extreme limit. It is a great plea- 
sure to come back to this problem which so occupied our thoughts during that 
month thirty one years ago and to contribute some new thoughts for Holger 's 
Festschrift. 

2 The Puzzle of Infinite Size 

The obvious difficulties with hadronic string theory involved the spectrum which 
invariably included massless vectors, scalars and tensor particles. There were 
also the subtle problems of local currents that Holger and I had wrestled with. 
Technically speaking, there was no possibility of continuing string theory off 
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the mass shell to construct the matrix elements needed to describe the interac- 
tion of hadrons with electromagnetism and gravitation [^|. The natural candi- 
dates, vertex operators like expikX can not be sensibly continued away from 
specific discrete "on shell" values of k 2 . Closely connected with this was the 
divergence encountered in attempting to compute the hadronic electromagnetic 
or gravitational radius [||] Thus string theory was abandoned as a theory of 
hadrons and replaced by QCD. The success of string theory in understanding 
Regge Trajectories and quark confinement was understood in terms of an approx- 
imate string-like behavior of chromo-electric flux tubes. According to this view, 
hadronic strings are not the infinitely thin idealized objects of mathematical string 
theory but are thick tubes similar to the quantized flux lines in superconductors Q . 
The ideal string theory was relegated to the world of quantum gravity. 

However more recent developments have strongly suggested that an ideal- 
ized form of string theory may exactly describe certain gauge theories which are 
quite similar to QCD [^] We have returned full circle to the suspicion that 
hadrons may be precisely described by an idealized string theory, especially in 
the 't Hooft limit The new string theories are certainly more complicated than 
the original versions and it seems very plausible that the problems with the mass- 
less spectrum of particles will be overcome. Less however has been studied about 
the problems connected with local currents. In this contribution I will show that 
the new insights from the AdS/CFT correspondence provide a solution to the 
form factor problem. 

I begin by reviewing the problem. For definiteness we work in the light cone 
frame in which string theory has the form of a conventional Galilean-invariant 
Hamiltonian quantum mechanics. The degrees of freedom of the first-quantized 
string include D — 2 transverse coordinates X m (cr) and the Lagrangian for these 
variables is 



L 4*71 



1 r 27tP_ 



do-(XX-(aV 2 X'X') (1) 

o 

where X and X' mean derivative with respect to light-cone time T and string pa- 
rameter cr. The light-cone momentum P is conjugate to the light like coordinate 
x~ . All irrelevant constants have been set to unity. 

An important feature of the light-cone theory involves the local distribution 
of P_ on the string. The rule is that the distribution of P_ is uniform with respect 
to cr. In other words the longitudinal momentum dP_ carried on a segment of 
string do" is exactly &a/2n. 

Let us now consider the transverse density of P_ . In a space- time field theory 
this would be given by 



P(X) 



dx~T__(X,x-) (2) 



where T is the energy momentum tensor of the field theory. Matrix elements of p 
between strings of equal P_ define form factors for gravitational interactions of 
the string and are entirely analogous to electromagnetic form factors. 
The natural object in string theory to identify with p(X) is 

-Lfdff[6(X-X(cr))] (3) 
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In other words p(X) receives contributions from every element of string localized 
at X. The Fourier transform of p(X) 



p(k) = 



dcrexpikX(cr) (4) 



defines a system of form factors by its matrix elements between string states. 
The mean square radius of the distribution function is given by 



R 2 



X 2 p(X)) (5) 



and can be rewritten in terms of p. 

R 2 = -9 k 9 k (p)| k=0 . (6) 

Eq.(||) is the standard definition of the mean-square radius in terms of the mo- 
mentum space form factor. 

The squared radius is also given by 

(X(a) 2 ). (7) 

where the value of cr is arbitrary. 

For a field theory with a mass gap, such as pure QCD it is possible to prove 
that R 2 is finite. This follows from the standard analytic properties of form fac- 
tors. The problem arises when we attempt to apply the world sheet field theory 
to compute (X(cr) 2 ). An elementary calculation based on the oscillator represen- 
tation of X gives a sum over modes 

oo - 

(X 2 ) - <x'Y_ ~ = a'logoo. (8) 
o n 

A related disaster occurs when we compute the form factor which is easily seen 
to have the form 

(p(k)) =exp-k 2 (X 2 ). (9) 

Evidently it is only non-zero at k 2 = 0. 

In a covariant description of string theory the problem has its roots in the 
fact that the graviton vertex operator is only well defined on the mass shell of 
the graviton, k 2 = 0. Vertex operators to be well defined must correspond to per- 
turbations with vanishing world sheet (3 function. This implies that they should 
correspond to on shell solutions of the appropriate space-time gravitational the- 
ory. For the kinematical situation in which the graviton carries vanishing k± the 
transverse momentum must vanish. Thus no well defined off shell continuation 
of the form factor exists. 

One might wonder if the divergence of X 2 is special to the case of a free world 
sheet field theory. The answer is that the divergence can only be made worse by 
interactions. The 2-point function of a unitary quantum field theory is at least as 
divergent as the corresponding free field theory. This follows from the spectral 
representation for the two point function and the positivity of the spectral func- 
tion. Thus it is hard to see how an ideal string theory can ever describe hadrons. 
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3 Light Cone Strings in AdS 

There are good reasons to believe that certain confining deformations of max- 
imally supersymmetric Yang Mills theory are string theories albeit in higher di- 
mensions. The strings move in a 5 dimensional space ^ that is asymptotically AdS. 
In the 't Hooft limit these theories are believed to be free string theories. Evidently 
if this is so there must exist a well defined string prescription for form factors in 
the 4-D theory. 

What we will see is that although the theory in bulk of AdS is an ideal thin- 
string theory the 4-D boundary field theory is not described by thin strings. That 
may seem surprising. Suppose that in the light-cone frame the thin 5-D string has 
the form 

X(ff),Y(oO (1) 

where X are the transverse coordinates of 4-D Minkowski space and Y is the ad- 
ditional coordinate perpendicular to the boundary of AdS. Then it would seem 
natural to consider the projection of the string onto the X plane to define a thin 
string. According to this view the mean-squared radius would again be (X 2 ) and 
we would be no better off than before. Before discussing the resolution of this 
problem let us work out the bosonic part of the light-cone string Lagrangian in 
AdS. I will make no attempt to derive the full supersymmetric form of the theory 
in this paper. I believe the resolution of the form factor problem does not require 
this. On the geometric side I will also ignore the 5-sphere component of the ge- 
ometry implied by the usual R-symmetry of the N = 4 supersymmetry. 
The metric of AdS is given by 

dsW^**" -f 2 - dY2 (2) 

I have defined the overall scale of the AdS (radius of curvature) to be R. 

In order to pass to the light cone frame we must also introduce the world 
sheet metric htj . In the usual flat space theory it is possible to fix the world sheet 
metric to be in both the light cone gauge o"o = T = x + and also the conformal 
gauge hoo = — Hii, hoi =0. However this is not generally possible since it 
entails 3 gauge conditions which is one too many. The special feature of flat space 
which permit the over-fixing of the gauge is not shared by AdS. Thus we must 
give up the conformal gauge if we wish to work in light-cone gauge. 

Let us fix the gauge by choosing 2 conditions 

ff = x + 

Hoi = 0. (3) 

Let us further define 



^ = E. (4) 
hoo 

2 Strictly speaking the target space is 1 dimensional with the form AdS 5 times a compact 

space such as S5 . In this paper the compact factor plays no role. 
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Setting a' = 1, an elementary calculation gives the Hamiltonian 
H 



R 4 



dcr (^PxPx + PyPy + ^[d a Xd a X + d a Yd a Y)j . (5) 

The precise version of the supersymmetric Hamiltonian was given in [^] . This is a 
more or less conventional string action with the unusual feature that the effective 
string tension scales like 1 /Y 4 . Thus the tension blows up at the AdS boundary 
Y = and tends to zero at the horizon Y = oo. This of course is a manifestation of 
the usual UV/IR connection . 

The Hamiltonian could be obtained from an action 

dadT (xX + YY -^(a«rX3 ff X + a ff Y3 ff Y)) . (6) 

This action thought of as a 1 + 1 dimensional field theory is not Lorentz invariant 
in the world sheet sense. However it is classically scale invariant if we assume 
X and Y are dimension zero. The Hamiltonian has dimension 1 and therefore 
scales as the inverse length of the cr circle which is ~ P . We recognize this scale 
symmetry as space-time longitudinal boost invariance under which H and P_ 
scale oppositely and X,Y are invariant. No doubt the actual Lagrangian when 
properly super-symmetrized retains this symmetry when quantized. 

Let us consider the equal time correlation function (X(O)X(cr)) in the field 
theory defined by @ or more precisely in its 

supersymmetrized version. By inserting a complete set of eigenstates of the 
energy and (world sheet) momentum we obtain 



(X(O)X(cr)) = Y_ 



• l ^|(X(0)lH,p)| 2 -^dHd P 



FfE.p^-^dEdp (7) 
with F > 

The measure of integration dEdp/p 2 follows from the fact that X has "engi- 
neering" dimension zero under the longitudinal boost rescaling. Furthermore the 
assumption that the scale invariance is preserved in the quantum theory requires 
F(E,p) = F(E/p) for large p, E. It follows that as long as F does not go to zero 
in this limit that the correlation function diverges as cr — > 0. This would imply 
X 2 = oo 



4 Dressing the Vertex with Y Dependence 

Let us consider the problem from the point of view of the vertex operator exp ikX. 
One problem that I have emphasized is that it is not a solution of the on-shell con- 
dition. We can try to fix this by replacing it with a solution of the wave equation 
for a graviton in AdS space. The relevant equation is 



(d^ + Y 3 9 Y Y- 3 3 Y ) O=0 



(1) 
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or 

52 _ , v 2 



where |x runs over the four dimensions of flat Minkowski space. 

The particular solutions we are looking for are independent of the x ± and 
have the form 

O = expikXF(k, Y) (2) 

where F satisfies 

Y 3 9 Y Y- 3 9 Y F(k,Y) = k 2 F (3) 
Thus the on shell vertex operator has the form 

dffexpikX(o-)F(k,Y(o-)) (4) 

The factor F is a dressing of the vertex, necessary to make its matrix elements well 
defined for k ^ 0. 

Let us consider the mean square radius of the hadron defined by eq.Q. 

R 2 =-9 k 3 k (F(k,Y)expikX)| k= o (5) 

R 2 = (X 2 F(0,Y) -2iX-F'(0,Y) -F"(0,Y)) (6) 
where F" = d^d^Y. 

For a state of zero angular momentum in the X plane the term linear in X 
vanishes and we have 

R 2 = (X 2 F(0,Y) — F"(0,Y)). (7) 

One possibility for resolving the infinite radius problem is a cancellation of 
the two terms in eq.(0). To compute F and F" we Taylor expand F(k, Y) in powers 
of k and substitute into eq.([|). There are two linearly independent solutions. 

F(k,Y) =Y 4 + ^k 2 Y 6 + --- (8) 

and 

F(k,Y) = 1 -lk 2 Y 2 + -.. (9) 

Only the second of these is relevant to the problem of vertex operators. To see 
this we need only note that the vertex at k = is just the operator that measures 
P_ . For states with P = 1 this operator in just the identity. This implies that 
F(0,Y) = 1. 

Thus we find 

R 2 = ([X 2 +Y 2 ]) (10) 

and R 2 is the sum of a divergent term and a positive term. The mean radius con- 
tinues to be divergent. Evidently cancellation is not the answer. 

The dressing of the vertex by the factor F(k, Y) obviously modifies the ex- 
pression (^) for the transverse density p. If we define the Fourier transform of F 
with respect to k to be F(X, Y) eq.@ is replaced by 

doF(X-X((r),Y). (11) 
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This means that an ideal thin string in the AdS bulk space is smeared out by the 
holographic projection onto the boundary. This is of course the familiar UV/IR 
correspondence at work. Bulk strings near the boundary are projected as very 
thin strings in the 4-D theory but those far from the boundary are fat. The extra 
term (Y 2 ) in eq.(^0|) represents this fattening. Evidently I have only made things 
worse by including the dressing. 

Before discussing the solution to the problem let us make some remarks 
about confining deformations in the context of AdS/CFT. Bulk descriptions of 
confining deformations of super Yang Mills theory have an effective infrared 
"wall" at a value of Y which represents the confinement scale. In these cases the 
metric @ is modified in the infrared region. 



where, as in the conformal case, h ~ 1 /Y 2 for Y — > 0. Assume that h has a mini- 
mum at the confinement scale, Y = Y* . 

The light-cone hamiltonian is easily worked out, 



Consider a string stretched along the X direction and choose cr so that d^X — 
1 . The potential energy of the string is then given by 



which has a minimum at Y = Y* . Thus a classical long straight string will be in 
equilibrium at this value of Y. This classical bulk string corresponds to a field 
theory configuration which, according to the UV/IR connection, is thickened to 
a size ~ Y*, that is, the QCD scale. 

Quantum fluctuations will cause the wave function of the string to fluctuate 
away from Y* . The implication is that the QCD string is a superposition of differ- 
ent thickness values extending from infinitely thin to QCD scale. Indeed different 
parts of the string can fluctuate in thickness over this range. The portions of the 
string near Y = will be very thin and will determine the large momentum be- 
havior of the form factor. 

5 Finiteness of (X 2 ) 

I believe that despite the argument given at the end of Section 3 the value of 
(X 2 ) is finite. This can only be if the function F = Y. l(X(0)|E,p)| 2 vanishes in the 
scaling limit of large E, p. I will first give an intuitive argument and follow it with 
a more technical renormalization group analysis that is due to Joe Polchinski. 

First suppose the string is "stuck" at some value of Y. In that case the action 
for X in eq.@ is a conventional string action except that the string tension is 
replaced by 1 /Y 4 . The divergence in X 2 would then be given by 




(12) 



H = dcr (P X P X + P y Py + h(Y) 2 (3 ff Xd (T X + 3 a Y9 a Y)) 



(13) 



V(Y) = h(Y) 2 



(14) 



(X 2 ) =Y 2 |loge|. 



(1) 
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If we ignore quantum fluctuations of Y we could replace Y by Y* . But Y fluctuates 
as well as X and can be expected to fluctuate toward the boundary as e tends to 
zero. This is just the usual UV/IR connection in AdS. Therefore as we remove the 
cutoff the fluctuations of X are diminished because the string moves into a region 
of increasing effective stiffness. If for example the average value of Y 2 tends to 
zero as 1 1 / log e | or faster then the fluctuations of X would remain bounded. To 
see that this happens we consider the renormalization running of the operator 
X 2 . 

Begin with the bare theory defined with a cutoff length e on the world sheet. 
We can then ask how a given operator in this bare theory is described in a renor- 
malized version of the theory with a cutoff at some longer distance I. A general 
operator 4>(X, Y) runs to lower momentum scales according to the renormaliza- 
tion group equation 

(ia/3l)<t>(X,Y,l) = (a'/2)V 2 4>(X,Y,l). (2) 

For example, consider flat space and the operator X 2 . We look for a solution of 
eq.(||) with 

4>(X,e)=X 2 . (3) 

The solution is 

4>(X,l] =X 2 + a'logl/e. (4) 

Thus if we regulate the theory at some fixed scale, for example I ~ 1 , the matrix 
elements of X 2 blow up as send e — > 0. 

By contrast, consider the the case of AdS space where 

V 2 = R~ 2 (Y 2 9£ + Y 5 3 Y Y- 3 3 Y ) . (5) 

For a solution of the form X 2 + f (l)Y 2 this becomes 

(l3/3l)f = (2a'/R 2 )(1 -f) . (6) 

With f (e) =0 the solution is 

f(l) = 1 - (e/l) l0L ' /Rl . (7) 

So if we fix the scale I and take the cutoff length e to zero the matrix elements 
tend to finite limits and the problem of infinite radii is resolved . If, however, we 
expand in powers of a' there are logarithmic divergences. 

Note that the operator X 2 runs to a fixed point X 2 + Y 2 which is just the 



operator in eq.fllO|) which represents the mean squared radius R 2 . 

The reader may wonder how the finiteness of X 2 can be explained in covari- 
ant gauges such as the conformal gauge in which the world sheet theory has the 
form of a relativistic field theory. A standard argument insures that the singular- 
ity in a two point function can not be less singular than a free field; in this case 
logarithmic. The argument is based on the positivity of spectral functions which 
in turn assumes the metric in the space of states is positive. In general this is not 
the case in covariant gauges. 
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6 Discussion 

The original attempt to describe hadrons as idealized strings was frustrated by 
the infinite zero point oscillations in the size of strings. Early ideas for modifying 
string theory such as replacing the idealized strings by fat flux tubes or as col- 
lections of partons which approximate strings fit well with QCD but seemed to 
preclude an idealized mathematical string description. 

More recent evidence from AdS/CFT type dualities suggest that idealized 
string theory in higher dimensions may provide an exact description of the 't 
Hooft limit of QCD-like theories. I have argued that an ideal bulk string theory in 
five dimensions is fully compatible with a fat non-ideal string in four dimensions. 

The fifth dimension can be divided into two regions. The "wall" region near 
Y = Y* corresponds to the confinement scale A. If we ignore high frequency 
fluctuations, the string spends most of its time in this region. The usual UV/IR 
spreading gives the string a thickness of order A. High frequency fluctuations of 
small sections of string can occur which cause it to fluctuate toward Y = 0, the 
region corresponding to short distance behavior in space- time. These fluctuations 
will control the large momentum behavior of form factors as well as deep inelas- 
tic matrix elements. Such fluctuations give the string a parton-like makeup. We 
have also seen that these fluctuations stiffen the effective string tension so much 
that the infinite zero point size that Holger and I worried about so long ago is 
now eliminated. 
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Relativity and c/ y/3 



S.I. Blinnikov, L.B. Okun and M.I. Vysotsky 
ITEP, 117218 Moscow, Russia 

Abstract. We define the critical coordinate velocity v c . A particle moving radially in the 
Schwarzschild background with this velocity, v c = c/ y/3, is neither accelerated, nor de- 
celerated if gravitational field is weak, r e <C r, where r g is the gravitational radius, while 
r is the current one. We find that the numerical coincidence of v c with velocity of sound in 
ultrarelativistic plasma, u s , is accidental, since two velocities are different if the number of 
spatial dimensions is not equal to 3. 

To Holger Nielsen 

We dedicate this note to our friend Holger Nielsen on the occasion of his 60th 
birthday. It has been a great pleasure to discuss with him exiting physical ideas 
at ITEP, CERN, the Niels Bohr Institute and in other places around the world. We 
wish to Holger many new discoveries and a long happy life. 

1 Motivation 

According to General Relativity (GR) clocks run slowly in the presence of gravi- 
tational field, as a result, the coordinate velocity of photons decreases. This is the 
reason for the delay of radar echo from inner planets predicted and measured by 
I. Shapiro Propagation of ultrarelativistic particles is described similarly to 
that of photons. That is why the retardation must take place not only for photons 
but also for ultrarelativistic particles. In this respect the latter drastically differ 
|^] from nonrelativistic bodies, velocity of which evidently increases when they 
are falling radially onto a gravitating body (e.g., onto the Sun). Obviously, there 
should be some intermediate velocity v c which remains constant for a particle 
falling in gravitational field of the Sun (or another star). The numerical value 
v c = c/\/3 will be found in Sect. 2. When a particle moves radially with this ve- 
locity in weak field it "ignores" gravity: it is neither accelerated, nor decelerated. 
For nonradial trajectories gravity is never ignored: the trajectories are bent for 
any velocity. 

It is well known that u s = c/a/3 is the speed of sound in ultrarelativistic 
plasma and the question arises whether the equality u s = v c has some physical 
reason, or it is a numerical coincidence. To answer this question we find in Sect. 3 
expressions for v c and u s in spaces with number of dimensions n different from 
3. Since for n^3 we get v c ^ u s we come to the conclusion that their coincidence 
at n = 3 does not have deep physical reason. 
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2 Derivation of v c = c/ "</3 

To simplify formulas, we put light velocity c = 1 , restoring it when it is necessary. 
In what follows G is gravitational constant; gravitational radius r g of an object 
with mass M equals 

r g =2GM. (1) 

Let us start from definitions used in GR. The expression for interval in the 
case of radial motion (d0 = dcp = 0) has the well known Schwarzschild form: 

ds 2 = goodt 2 - g TT dr 2 = dx 2 - dl 2 , (2) 

where goo = (grr) -1 = 1 — ^f-. The local velocity v of a particle measured by a 
local observer at rest is: 

v = il = f^V /2 ^ = -L^ , (3) 
dT \9oo/ dt goo dt 

while observer at infinity, where goo(oo) = g rr (oo) = 1, measures the so-called 
coordinate velocity at r: 

dr / goo 



v = -=v^— ) =g 00 v. (4) 

In order to determine the time of radial motion from a to b, the infinitely distant 
observer should calculate the integral 



b dr 



(5) 



that is why the coordinate velocity is relevant for radar echo. 

For a particle moving in static gravitational field one can introduce con- 
served energy (see ref. g, eq. 88.9): 



E=m^|. (6) 
Vl — v z 



The expression for E through 



v: 



E = (7) 
V 7 ! - (v/goo) 2 



allows us to determine v(r) from the energy conservation: 

E(r = oo)=E(r) , (8) 

v 2 = goo - 9oo + Qoo v i = 9oo H - gooO - v£j] . (9) 
For the local velocity v measured by a local observer we obtain: 

v 2 = l-good-v 2 ) , (10) 
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so, while v always increases for a falling massive particle, reaching c at r = r g , 
the behaviour of v is more complicated. Substituting goo = 1 — -f - into @, we get 
for weak gravitational field (r > T g ): 

v 2 =v^ + ^(l-3v^) . (11) 

For the motion of a nonrelativistic particle (v,^ <C 1 ) the well-known expression 
is reproduced: 

2 2 2MG 
v 2 = + — . (12) 

For Vqo = v c = 1 1 \f% the coordinate velocity of particle does not change, while it 
grows for v^ < v c and diminishes for Voo > v c . At r = 3r g according to eq.([ll|) 
the coordinate velocity becomes equal to v c . However, for r = 3r g our weak field 
approximation fails. 

Let us dispose of the assumption of the weak field. Coming back to expres- 
sion @ and substituting there goo = 1 — we observe that for > v c the 
coordinate velocity always diminishes and becomes zero at r = r g , while in the 
case Vqo < v c it grows up to the value 

v 2 max =4/(27(1 -v^) 2 ) , (13) 

which is reached at 

and after that diminishes to zero at r = r g . It is interesting to note that the velocity 
v measured by local observer equals v c at the point where v = v max . 

Thus, if the coordinate velocity is only mildly relativistic, Voo > c/V3, then v 
already decreases at the free fall. 

As an example of a non-radial motion let us consider the deflection of light 
from a star by the Sun and compare it with the deflection of a massive particle. It 
is well known that the angle of deflection 9 of photons grazing the Sun is given 
by 

0y = P- (15) 

where is the radius of the Sim. In the case of massive particles the deflection 
angle is larger: 

e = e Y (i + r 2 ) , (16) 

where (3 = /c < 1 .(See ref.[[§, eq. 25.49, and ref.[[|, problem 15.9, eq.13.) 



3 Speed of sound u s and critical speed v c in n dimensions 



For ultrarelativistic plasma with equation of state P = e/3, where P is pressure 
and e is energy density (including mass density), we have for the speed of sound 
u s : 



9P 



(17) 



ad 
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We use eq. 134.14 of ref . [[51], and correct misprint there, or eq. 126.9 from [0]; "ad" 
means adiabatic, i.e. for constant specific entropy. In order to obtain the expres- 
sion for u s in the case when n ^ 3, where n is the number of spatial dimensions, 
let us start with equation of state. 

One can use a virial theorem to connect pressure P and thermal energy £ of 
an ideal gas using classical equations of particle motion (cf. []|]). We have for a 
particle with momentum p and a Hamiltonian H: 

<-=-£• < 18 > 

9q 

hence, 

3H 

qp = -q-j^ = qF , (19) 
where F is the force acting on the particle. Let us average over time t: 

(...) = - f\..dt . (20) 
t Jo 

Integrating by parts we get: 

(qp) = -<qp> = (qF) . (21) 
For non-relativistic (NR) particles 

qp = 2E kin = p 2 /m - (22) 
For extremely relativistic (ER) particles 

qp = Ekin = c|p| . (23) 

Now for N particles in a gas 

N N 

-2jqipi) = 2j* F ^ < 24 > 

i=1 i=1 

(By the way, — A £j (qiFt) is called the virial.) If the gas is ideal (i.e. non-interacting 

i 

particles), then the force F is non-zero only at the collision of a particle with the 
wall, and the virial reduces to an integral involving pressure: 



N 



L<qp> 

t=i 



Pnq dS 



divqdV = -3PV, (25) 



where n is a unit vector normal to the wall area element dS and the Gauss theo- 
rem is used for transforming the surface integral to the volume one. So, since the 
thermal energy £ (not including mass) is just the total kinetic energy of molecules, 



NR : 28 = 3PV, P = 2£/{3V) , 



(26) 
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ER : £ = 3PV, P = £/{3V) = e/3 . (27) 

The last equality holds since in extremely relativistic case ^> m. We see 
that 3 here is due to div q = 3, i.e. the dimension of our space. 

In a space of n dimensions, following the same lines, we get div q = n, so 
P = e/n and for ER gas we obtain: 

u s = c/V^ . (28) 

Here we should use n-volume V n instead of V = V3 and postulate the first law of 
thermodynamics for adiabatic processes to be d£ + PdV n = 0, so pressure would 
be the force per unit V n _j - the boundary of V n . 

The same equation of state follows from consideration of the stress tensor 
Ttic of ultrarelativistic plasma, which is diagonal and traceless in a rest frame of 
plasma: T o = e, T u = P = e/n. 

In order to find v c in the case n ^ 3 we need an n+ 1 -dimensional spherically- 
symmetric static generalization of the 3 + 1 -dimensional Schwarzschild metric 
which was found by Tangherlini [[J. (See refs.[|lO|] for the details of aspherical 
and time-dependent black holes metrics in higher dimensional spacetimes). The 
line element of the n + 1 -dimensional Schwarzschild metric is 

ds 2 = (l " ^) & ^ - r2dn - 1 ' (29) 

where dO( n _i ) is the line element on the unit (n— 1 )-sphere and the gravitational 
radius r gn is related to the black hole mass M: 

2 167tG n M 

Here A n _i denotes the area of a unit n — 1 sphere, which is 2n^ /F(y ) (for n = 3: 
T(3/2) = \/tc/2, A2 = 4tc). We consider the spaces with n > 3. The factor in the 
definition of r gn is taken from refs.[[l0|]. The form of the metric ( p9| ) is very easy 
to guess. In weak field limit, when goo — > 1 + 2cp we have 

y.n—2 

<P = ~ 2r™- 2 f ° r T °° ' ^ 

This leads in a natural way to the gravitational acceleration g with the radial 
component 

3cp (n-2)r^ 2 
9n = = t — for r -> 00 , (32) 



which implies the constant flux of the acceleration g equal to 

A n _ 1 r n - 1 g n =87t^G n M (33) 
n — 1 

through a sphere of area A n _ir n_1 at large r. It is not hard to verify that the 
Ricci tensor is zero for the metric (p9|), that is the metric (p9|) satisfies Einstein 
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equations in vacuum and describes a spherically symmetric spacetime outside a 
spherical gravitating body. 

One should remember that the dimension of G n depends on n. It is clear 
in the weak field limit from eqs. ( ^l| ) and ([301), since the dimension of [cp] is the 
square of velocity, i.e. zero for c = 1, and hence 

[G n ] = ^ , (34) 

or [GJ =L n - 1 ,if [M] = L _1 . 

For the coordinate velocity of a radially falling particle we get from eq. @ 
for weak gravitational field: 

v 2 ^ + (I^r 2 n-3v^ (35) 

instead of eq. (fTT|). We see that in the case of n-dimensional space the expression 
for v c remains the same, v c = c/ V3- Number "3" here is not due to the dimension 
of space, it is simply due to cubic polynomial in (||). 



4 Conclusions and Acknowledgements 

The speed of sound in relativistic, radiation dominated plasma depends on the 
dimension of space, while the critical velocity v c = c/ y/3 in the Schwarzschild 
metric is the same for any dimension. 
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Finding Center Vortices 

J. Greensite* 

Physics and Astronomy Department, San Francisco State University San Francisco, CA 



Abstract. I report on recent progress in locating center vortex configurations on thermal- 
ized lattices, generated by lattice Monte Carlo simulations of SU(2) gauge theory. A very 
promising method, which appears to have some important advantages over previous tech- 
niques, is center projection in the direct Laplacian center gauge. 

1 Introduction 

It is a great pleasure for me to contribute this article in celebration of Holger 
Bech Nielsen's 60th birthday. As it happens, my contribution to these proceed- 
ings is quite closely related to ideas developed some twenty years ago by Holger, 
Poul Olesen and other members of the Niels Bohr Institute Holger and his 
co-workers put forward a model of the QCD ground state known as the "Copen- 
hagen Vacuum," which is an explicit realization of the center vortex theory of 
confinement, propounded at about the same time by 't Hooft, Mack, and oth- 
ers The Copenhagen Vacuum is based on a calculation of the QCD effective 
action by techniques of continuum perturbation theory. The finding of Holger 
Nielsen and Poul Olesen, based on an analysis of this effective action, was that 
the QCD vacuum is dominated by a tangle (or "spaghetti") of center vortices. 
Unfortunately, this finding did not settle the issue of confinement, since there are 
always doubts about the validity of perturbation theory in the context of infrared 
physics. After a brief rush of interest, very little work was done on either the cen- 
ter vortex theory in general, or the Copenhagen vacuum in particular, after the 
early 1980's. 

Interest in the center vortex theory revived almost twenty years later, in 1997, 
when methods were developed for locating vortex positions in lattice gauge con- 
figurations, generated numerically by the Monte Carlo technique This ad- 
vance led to another burst of activity, this time mainly numerical, which accumu- 
lated a great deal of evidence in favor of the vortex theory. However, the methods 
which are used to locate vortices in lattice configurations have recently come un- 
der critical scrutiny. In particular, what are these methods really doing, and can 
they fail? A number of authors have pointed out certain weaknesses, and these 
are related to the problem of Gribov copies. In this article I would like to describe 
the original method for vortex location and the difficulties that have been encoun- 
tered, and to discuss an improved procedure which overcomes those difficulties. 

* E-mail: greensit@quark.sfsu.edu 



94117 USA 



Finding Center Vortices 123 



The work I will report was done in collaboration with Manfried Faber and Stefan 
Olejnik. 

2 Thin Center Vortices 

A center vortex is an object whose field strength is concentrated in a point-like re- 
gion in D=2 dimensions, a line-like region in D=3 dimensions, and a surface-like 
region in D=4 dimensions. When a vortex is created on an arbitrary background, 
such that it is topologicalry linked to some closed curve C, then a Wilson loop 
around that curve is changed by a center element of the gauge group: 



PexpIioA^x™] — » zPexp[ioA H dx 



(1) 



where z e Zn is an element of the center subgoup of the gauge group SU(N). 

To create a thin center vortex on a two-dimensional lattice, make the follow- 
ing transformation on the links 

U (x)^zUoM x =0, *i >0 (2) 

as illustrated in Fig. [|. The action changes only locally, at the shaded plaquette P, 
but Wilson loops U(C) change globally by a center element, if the loop C encloses 
the plaquette P. The create a thin vortex in higher dimensions, make the trans- 
formation shown on every xq — xi plane. The stack of shaded plaquettes forms a 
line-like object, in the xz direction, in D=3 dimensions, and a surface-like object 
in D=4 dimensions. 



W(C) 

Fig. 1. Creation of a center vortex (shaded plaquette) in a plane. 



Since creation of a thin vortex changes plaquettes by 

TrtUUU^t] -» zTr[UUU t U t ] (3) 

on the vortex sheet (D=4), the action of a thin vortex is singular in the continuum 
limit. If center vortices are found in the QCD vacuum, then we expect these to 
be extended objects, in which the action of of the vortex sheet is spread out into 
a surface of finite thickness in physical units. It is not hard to show that vortices 
percolating through the lattice will disorder Wilson loops, leading to an area law. 
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3 Direct Maximal Center Gauge 

Although it is easy to describe how to create center vortices on an arbitrary back- 
ground, it is not obvious how to identify such objects in any given configuration. 
The method which our group proposed in ref . [f|] was to fix the lattice to a certain 
gauge, known as "direct maximal center gauge," and then map the links of the 
gauge-fixed lattice onto center elements, a procedure known as "center projec- 
tion." An important insight into what this gauge is doing — and how it can fail — 
is due to Engelhardt and Reinhardt [||] . 

A typical thermalized lattice U H (x), if printed out, looks like a set of random 
numbers. But it is not random, and in fact locally, at (3 ^> 1, the link variables 
approximate a (classical) vacuum configuration 

U H (x)«g(x)gt(x+fi) (4) 

Let us ask the question: What is the best fit, to a given lattice U H (x), by a pure 
gauge g(x)g^(x + ft)? We work with the SU(2) gauge group for simplicity. The 
best fit will minimize 

d ' = 4^Z Tr [( U ^ x) ~ 9Wgt(x+W ) x ( h ' c ) 

x,n 

= ^^Tv[l-gHx)U^x)g(x + iL)] (5) 

where V is the lattice volume. Define 

9iyx) = gt(x)iyx)g(x + fl) (6) 
Minimizing dp is the same as maximizing 

R L = ^Tr[ 9 U^x)] (7) 

which is the same as the lattice Landau gauge condition. What this demonstrates 
is that lattice Landau gauge is is equivalent to a best fit to a pure gauge. 

Next, consider trying to make a best fit to a thin vortex configuration on a 
classical vacuum background 

V^(x) = g(x)Z^x)gt(x + fl) (8) 

where Z H (x) = ±1 . This can be done in two steps. First, since the adjoint represen- 
tation is blind to Z H , we start by finding the best fit in the adjoint representation, 
minimizing 



d A = 4^X Tr A [(U H (x)-g(x)g t (x + W) x (h.c.)~ 



4V 



^2Tr A [l-gt(x)lV(x)g(x + M] 



(9) 
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which is the same as maximizing 

Rdmc = ^Tr A [ s u^{x)] ( 10 ) 



x,M- 



where Tr a indicates the trace in the adjoint representation. This is the direct max- 
imal center gauge-fixing condition, which is obviously just lattice Landau gauge 
in the adjoint representation. Having found g(x) maximizing Rdmc/ we then find 



the Zp_(x] which minimizes 



d 2 = — ^Tr[(lVx)-g(x)Z^x)gt(x + £)) x (h.c.)] (11) 

in the fundamental representation. For fixed g(x), this minimization is achieved 
by setting 

Z^(x)=signTr[eiyx)] (12) 

This step is is known as "center projection." It maps an SU(2) lattice to a Z2 lattice. 
The original and gauge-transformed lattices can be expressed, respectively 

9 U (l (x)=Z K (x)e i6A »' (x) (13) 

with Zp.M the vortex background, and 8A^{x) the fluctuation around the back- 
ground. Direct maximal center gauge finds the optimal Z M .(x) minimizing 6A M .(x). 

Unfortunately, as also pointed out in ref. there are going to be problems 
with this approach to vortex finding, particularly in the continuum limit. It is 
clear that the action density of a thin vortex is singular at "P-plaquettes", which 
are plaquettes p where the product of Z^(x) around the plaquette is negative, i.e. 
Z(p) = —1. P-plaquettes are dual to the thin vortex surface. Since a thermalized 
lattice configuration is locally close to a pure gauge at large (3, i.e. ATr[U(p)] w l,it 
follows that the fit to a thin vortex is very bad at P-plaquette locations. As |3 — > 00, 
the bad fit near P-plaquettes may overwhelm the good fit in the exterior region, 
and the best fit in the continuum limit may be no thin vortices at all. 

Computer simulations tend to support this analysis. In practice it is impossi- 
ble to find the global maximum of Rl or Rdmc; instead there are various methods 
of finding local maxima, known as Gribov copies. Using a sophisticated gauge- 
fixing technique which employs simulated annealing, Bornyakov et al. [g] were 
able to generate Gribov copies in direct maximal center gauge with consistently 
higher values of Rdmc than those achieved by previous approaches based on 
the over-relaxation method. However, these "better" Gribov copies had much 
worse center dominance properties in the scaling regime. Center dominance is the 
agreement between string tensions on the full and center projected lattices. Using 
the improved gauge fixing techniques, Bornyakov et al. found that the projected 
string tension is lower than the full asymptotic string tension by about 30%. 

4 Direct Laplacian Center Gauge 

In confronting the difficulties raised in refs. [||@], we would like to retain the strat- 
egy of finding the best fit to a vortex configuration, but avoid over-emphasizing 
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the inevitable bad fit at P-plaquettes, which seems to be the source of those diffi- 
culties (cf. ref. We therefore adopt a strategy inspired in part by related work 
of de Forcrand et al. jUJ, which in turn draws on the earlier Laplacian Landau 
gauge fixing of Vink and Weise [|l0|]. The idea is to replace the SO(3) gauge trans- 
formation g(x) in (10) by a real 3 x 3 matrix-valued field M(x), chosen to maximize 

Rm = Y_ Tr[M T (x)U A ^(x)M(x + fx)] (14) 

(Ua is the link in the adjoint representation) with the constraint that M(x) is or- 
thogonal "on average." 

l^M T (x)M(x) =1 (15) 

X 

This essentially introduces an adjustable weighting; M(x) is allowed to be small 
where the fit to a thin vortex is bad. 

Rm is maximized (uniquely!) by finding the three lowest eigenvalues, and 
corresponding eigenfunctions, of the lattice Laplacian eigenvalue equation (no 
sum over a = 1 , 2, 3) 

^V M {x,y)^{y)=\ a f^(x) (16) 
v 

where 

£>ij(x,y) = 

- Y_ (PWxJhj^.x+ii + [U At ,(x)]j i 6 v , x _f 1 - 26 xv 6 ab ) . (17) 

There are standard numerical algorithms which solve this equation, and the ma- 
trix field formed by 

M ab (x) = f*(x) (18) 

where b = 1,2,3 denotes the three eigenvectors with lowest eigenvalues, is the 
matrix field maximizing Rm- 

We then map M a b (x) to the nearest SO(3) matrix 

[g A (x)] tJ =fl(x) (19) 

which also satisfies a Laplacian equation 

^XMx.yJfffo) =A ac (x)f?(x) (20) 
v 

This is equivalent to finding the SO(3)-valued matrix field closest to M a b(x) 
which is also a local maximum of Rdmc- The mapping is accomplished in two 
steps: 

1. Polar decomposition 

M(x) =±0(x)P(x) (21) 
to find the SO(3) matrix D which is closest to M. 
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2. Relaxation Q. — > gA to the nearest local maximum of 

Rdmc =^Tr[g A ( x )U A ^Mg A (x + ft)] (22) 

Finally, transform U H (x) — > 9 U |X (x). This prescription for mapping an arbitrary 
gauge field configuration to an (in principle unique) point on the gauge orbit is 
called direct Laplacian center gauge. We apply center projection 

Z H (x) = signTrPiyx)] (23) 

as before, to locate the P-plaquettes forming thin vortices (known as "P-vortices") 
on the projected lattice. 

4.1 Center Dominance 

The most important property to check with the new procedure is center domi- 
nance, since this is where problems were found numerically, in refs. [^] and [pT|, 
with the previous methods. 

The center vortex contribution to the asymptotic string tension is extracted 
from 

W cp (C) = (ZZ...Z) (24) 

in the center-projected configurations, obtained after direct Laplacian gauge fix- 
ing. Creutz ratios Xcp [R> R] are plotted in Fig. || The straight line is the asymptotic 
freedom behavior, with \/a/A — 58 

Note that the Creutz ratios at each |3 are nearly the same (with the exception 
of the 1 x 1 loop), and this is known as "precocious linearity." This feature can 
also be clearly seen in Fig. [| Here we also see the property of center dominance: 
The center-projected Creutz ratios Xcp [R, R] a t all R > 2 agree quite well with the 
full asymptotic string tension cr, extracted from unprojected Wilson loops, both at 
strong and weak couplings 

An alternative way of displaying our data is to plot the ratio 

Xcp [R, R] _ Xpht;s[R,R] ^ 

as a function of R in physical units, where o"Lat is the asymptotic string tension in 
lattice units, while the subscript "phys" refers to physical units. In Fig. [| we see 
that the ratios are not far from unity (center dominance) even down to relatively 
small distance scales (precocious linearity). 

4.2 Other Tests 

To check that P-vortices are really associated with center vortices in the full con- 
figuration, and are not just artifacts of the projection, it is necessary to verify their 
correlation with gauge invariant observables. Such tests were invented in refs. 
|^,U] and Jl3|], and it is important to repeat these with the new method. Here I 
will only touch on the results. 
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Fig. 2. Projected Creutz ratios of various sizes as a function of |3. 



We define the "vortex-limited" Wilson loops W n ( C) to be loops evaluated on 
a subensemble of configurations, selected such that precisely n P-vortices, in the 
corresponding center-projection, pierce the loop. Center projection is used only to 
select the data set. The loops themselves are evaluated using the full, unprojected 
link variables. What we find numerically is that in the limit of large loops, 



W n (C) 

Wb(C) 



(26) 



This is what one expects, if P-vortices locate thick center vortices in the unpro- 
jected SU(2) lattice. 

It is also possible to measure the density of P-vortices as a function of (3. If 
N vor is the number of P-vortex plaquettes, and N T is the total number of plaque- 
ttes on the lattice, then the density of P-vortices in lattice units is given by 



,2 



N T N T a 4 
Total Vortex Area 
6 x Total Volume 



Q 2 = 



(27) 



where p is the center vortex density (area per unit volume) in physical units. 
According to asymptotic freedom, if the P-vortices are locating physical objects, 
we should have 
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Fig. 3. Center projected Creutz ratios for various couplings and lattice sizes. Straight lines 
show the values of the asymptotic string tensions and corresponding error bars, in the 
unprojected configurations, quoted in ref. [Q. 
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But also p is related to the center-projected one-plaquette loop 

W cp (1,1) = (1-p)x(+1)+px(-1) = 1-2p (29) 

which allows us to calculate p from center-projected plaquette values. It is found 
that the the data for p is quite consistent with the asymptotic freedom result (^8|) . 

Finally, as suggested by de Forcrand and D'Elia [[13]], one can remove center 
vortices from the gauge-fixed configuration by simply multiplying by the pro- 
jected configuration 

9 U,(x) -> 9 U' M = Z,(x)9U,(x) (30) 

What we find, using the direct Laplacian gauge procedure, is that confinement 
disappears when vortices are removed in this way. 



5 Conclusions 



The conclusion we draw from all this is that center vortices in the Yang-Mills 
vacuum are correctly located by the new procedure, and these vortices are re- 
sponsible for the confining force. Direct Laplacian center gauge combined with 
center projection is found to be an improvement over the original approach to lo- 
cating center vortices in lattice configurations, and presumably works well even 
in the continuum limit. It should be a useful tool for the further study of the cen- 
ter vortex "spaghetti vacuum;" a picture of the QCD vacuum very similar to that 
proposed by Holger Nielsen and Poul Olesen back in 1978. 



Acknowledgements 

Much of this work, and work that preceded it, was carried out in Copenhagen 
at the Niels Bohr Institute, whose hospitality I am happy to acknowledge. This 
research is supported in part by the U.S. Department of Energy under Grant No. 
DE-FG03-92ER40711. 



References 



H. B. Nielsen and P. Olesen, Nucl. Phys. B160 (1979) 380 

G. 't Hooft, Nucl. Phys. B138 (1978) 1; G. Mack, in .Recent Deveiopmenfs in Gauge 
Theories, edited by G. 't Hooft et al. (Plenum, New York, 1980). 

L. Del Deb bio, M. Faber, J. Greensite, and S. Olejnik, Phys. Rev. D55 (1997) 2298, h ep- 
lat/9610005. 

L. Del Debbio, M. Faber , J. Giedt, J. Greensite, and S. Olejnik, Phys. Rev. D58 (1998) 



094501, hep-lat/980102/ 



M. Engelhardt and H. Reinhardt, Nucl. Phys. B567 (2000) 249, |hep-th/9907139 . 

V. Bornyakov, D. Komarov, and M. Polikarpov, Phys. Lett. B497 (2001) 151, hep- 

lat/0009035. 



7. M. Faber, J. Greensite, and S. Olejmk, hep-lat/0103030 



Finding Center Vortices 131 



C. Alexandres M . D'Elia, and Ph. de Forcrand, Nucl. Phys. Proc. Suppl. 83 (2000) 437, 



hep-lat/9907028 . 



9. Ph. de Forcrand and M. Pepe, Nucl. Phys. B598 ( 2001) 557,|hep-lat/0008016 . 

10. J. Vink and U. Weise, Phys. Lett. B289 (1992) 122, frep-lat/9206006[ 

11. T. Kovacs and E. Tomboulis, Phys. Lett. B463 (1999) 104. 

12. C. Michael and M. Teper, Phys. Lett. B199 (1987) 95; 

G. Bali, C. Schlichter, and K. Schilling, Phys. Rev. D51 (1995) 5 165,|hep-lat/9409005 . 

13. Ph. de Forcrand and M. D'Elia, Phys. Rev. Lett. 82 (1999) 4582, |hep-lat/9901024 



Blejske Delavnice Iz Fizike, Letnik 2, ST. 2, ISSN 1580-4992 
Bled Workshops in Physics, Vol. 2, No. 2 

Zbornik delavnic 'What comes beyond the Standard model', 2000 in 2001 

Zvezek 1: Zbornik ob 60. rojstnem dnevu Holgerja Bech Nielsena 

Proceedings to the workshops 'What comes beyond the Standard model', 2000 
and 2001 

Volume 1: Festschrift dedicated to the 60th birthday of Holger Bech Nielsen 

Uredili in oblikovali Norma Mankoc Borstnik, Colin D. Froggat in Dragan 
Lukman 

Publikacijo sofinancira Ministrstvo za solstvo, znanost in sport 
Tehnicni urednik Vladimir Bensa 

Zalozilo: DMFA - zaloznistvo, Jadranska 19, 1000 Ljubljana, Slovenija 
Natisnila Tiskarna MIGRAF v nakladi 100 izvodov 



Publikacija DMFA stevilka 1479 



